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1.0  SUMMARY 


This  final  report  on  research  included  within  the  Tunable  High  Brightness  Semiconductor  Sources  work  unit 
[202D1100]  includes  several  technology  advancements.  First,  theoretical  advances  in  mid- wave  infrared 
(MWIR)  type-I  quantum  well  laser  efficiency  improvement  are  presented  utilizing  deep  understanding  of 
Auger  recombination  effects  and  how  to  change  them  through  proper  laser  design.  Experimental  results  are 
presented  to  confirm  this  effect  at  a  wavelength  of  approximately  2  fim.  Future  directions  are  presented, 
including  metamorphic  buffer  layer  grown  material.  Second,  a  surface-emitting  distributed  feedback  type-II 
quantum  well  laser  is  introduced  and  experimental  advances  are  presented.  Graphene-based  nonlinear  optical 
materials  are  introduced  for  their  application  to  broad-band  laser  technologies,  as  well  as  associated 
functions  such  as  sensor  protection.  Advances  in  solder  technology  for  optoelectronics  (including,  but  not 
limited  to  mid-IR  lasers)  is  presented  with  some  novel  high  density  indium  with  extraordinary  surface  and 
crystal  features.  Finally,  a  unique  intracavtiy  difference-frequency  generation  (DFG)  semiconductor  laser  is 
introduced  and  demonstrated. 
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2.0  INTRODUCTION 


There  are  few  light-weight,  efficient,  optical  sources  available  in  the  ~2-4  fim  band.  Presently,  the  two 
biggest  competitors  in  this  arena  are  Cr-doped  solid-state  lasers,  and  non-linear  frequency  converted  sources 
(both  of  which  are  being  pursued  within  Sensors  Directorate  (AFRL/RY)).  The  former,  while  reasonably 
efficient  (~15%  wall  plug  efficiency,  not  including  cooling),  is  bulky,  and  is  most  effective  in  the  2.3-3  fi m 
band.  The  latter  necessitates  the  complexity  of  a  multi-element  system,  though  does  allow  broad  spectral 
tuning  utilizing  a  single  source.  There  are  other  solid-state  approaches  that  may  be  suitable  for  direct  longer 
wavelength  generation,  they  are  complicated  by  cryogenic  cooling  which  adds  to  the  package  complexity  and 
reduces  overall  efficiency.  While  semiconductor  lasers  have  addressed  some  of  this  spectral  region,  with  a 
few  notable  exceptions,  they  require  cryogenic  cooling,  affecting  the  necessary  power  budget  for  the 
component. 

There  are  important  infrared  countermeasures  (IRCM)+  bands  I,  II  and  IV  within  this  region  that  are 
presently  addressed  only  through  nonlinear  frequency  conversion.  Sources  operating  in  the  MWIR  region 
also  are  important  for  laser  radar  (LADAR)  applications  because  the  wavelength  range  between  3.5  and 
3.9  fin l  is  where  diffraction,  scattering,  and  divergence  are  optimally  balanced.  For  similar  reasons,  this  same 
band  would  be  critical  to  ground-based  secure  free-space  communication  systems;  these  applications  remain 
critical  to  attain  Department  of  Defense  (DoD)  needs.  If  these  applications  are  enabled  with  compact, 
efficient,  in-band  mid-IR  sources,  it  is  possible  to  proliferate  sensor  protection  (layered  sensing:  assured 
global  access),  search  and  target  identification  (layered  sensing:  engagement  quality  information),  as  well  as 
the  jamming  (layered  sensing:  spectrum  dominance  and  control)  in  these  same  bands,  not  only  to  the  highest 
value  assets  within  the  theater  of  operations,  but  to  lower  value  assets  as  well,  assuring  complete  sensor 
coverage  of  the  battlespace.  Furthermore,  with  such  sensor  abundance,  sensor  networking  and  multi-modal 
sensing  become  an  added  benefit.  This  may  be  embodied  by  bi-static  in-band  illumination/sensing,  chemical 
sensing,  etc. 

short-wave  infrared  (SWIR)  and  MWIR  lasers  are  proving  to  have  many  applications  in  the  sensing 
community,  including  remote  chemical  sensing  of  numerous  organic  compounds  such  as  methane[l],  an 
ideal  spectral  regime  for  low-altitude  laser  radar[2]  and,  if  sufficiently  bright,  infrared  countermeasures. 
Although  there  are  several  potentially  successful  approaches  to  developing  semiconductor  sources  in  the 
short-wave  and  mid- wave  infrared,  as  one  approaches  the  2-4  fim  band,  challenges  exist  that  must  be 
overcome  for  successful  adoption  as  a  product.  Type-I  quantum  well  lasers  are  the  mainstay  of  the 
semiconductor  laser  industry,  having  dominated  early  on  and  have  been  used  to  achieve  a  broad  wavelength 
regime.  However,  performance  of  type-I  quantum  wells  suffer  at  longer  wavelengths,  due  primarily  to  an 
increased  Auger  recombination  rate.  We  chose  to  approach  the  problem  by  addressing  the  Auger 
recombination  mechanisms  directly  and  demonstrate  this  capability  in  the  region  of  A  =2  fim. 

In  Section  5.0,  we  introduce  graphene  as  a  nonlinear  optical  element.  We  present  the  exploration  of  nonlinear 
behavior  of  graphene  composites  in  Section  5.1,  which  has  applications  to  sensor  protection.  For  additional 
laser  functionality,  the  “open  cavity”  nature  of  the  vertical  external  cavity  surface  emitting  laser  (VECSEL) 
provides  an  extremely  attractive  platform  as  tunable,  high  power  laser  sources  with  nearly  diffraction-limited 
beams.  [3,  4]  The  VECSEL  is  considered  to  be  a  hybrid  between  diode  pumped  solid  state  laser  and  a 
semiconductor  laser.  Essentially  it  is  comprised  of  a  vertical  cavity  surface  emitting  laser  (VCSEL)  with  the 
top  distributed  Bragg  reflector  (DBR)  removed  and  replaced  by  an  external  curved  mirror.  The  subsequent 
insertion  of  an  appropriate  saturable  absorber  with  fast  response  time  leads  to  the  production  of  ultrafast, 
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high  power  pulses.  Section  Section  5.2  presents  design  and  experimental  demonstration  of  graphene-based 
absorbers  to  mode-lock  VECSELs.  Graphene  saturable  absorbers  have  previously  been  used  to  mode-lock 
both  solid  state  and  fiber  lasers.  However,  unlike  a  VECSEL,  these  are  able  to  withstand  reasonable 
absorption  without  any  detrimental  effects.  VECSELs  have  been  more  challenging  to  integrate  with  graphene 
due  to  their  high  sensitivity  to  losses.  In  an  attempt  to  solve  this  issue,  we  report  the  combine  the  saturable 
absorption  demonstrated  and  studied  by  layer(s)  of  graphene  with  the  antiresonant  saturable  absorber.  [5,  6] 
This  results  in  a  low  loss  cavity  with  a  high  saturation  fluences  as  well  as  ultrashort  pulses. 

Such  a  laser  system  is  therefore  quite  attractive  for  several  Air  Force  related  applications  such  as  LADAR, 
IRCM,  and  chemical  sensing  and  identification.  Air  Force  applications  necessitate  compact  and  efficient 
sources  that  address  multiple  functions  since  air  platforms  limit  the  accommodation  of  many  individual 
components  for  each  mission.  Thus,  VECSELs  with  tunable  pulse  widths,  wavelengths,  and  repetition  rates 
provide  the  necessary  functionality  to  carry  out  several  missions  in  one  compact  device. 

Because  VECSELs  can  access  a  wide  range  of  wavelengths,  their  use  in  LADAR  systems  for  acquiring, 
tracking,  and  classifying  various  targets  is  possible.  By  using  pulsed  transmission,  which  may  help  minimize 
thermal  and  environmental  propagation  effects,  lasers  operating  in  the  ultraviolet  (UV)  to  infrared  can  detect 
any  scattered  radiation  and  process  them  in  a  laser  radar  receiver  to  remotely  determine  target  characteristics. 
Further,  proactive  IRCM  using  LADAR  systems  can  allow  for  protection  from  weapon  tracking  systems. 
Various  chemical  and  biological  agents  possess  characteristic  resonances  which  span  over  a  wide  wavelength 
range.  Using  a  mode-locked  laser  in  a  cavity  ring  down  spectroscopy,  where  the  center  wavelength  may  be 
tuned  to  the  particular  chemical  of  interest,  can  be  of  great  advantage  to  rapid  chemical  analysis  systems. 
This  would  allow  one  to  probe  necessary  vibrational  modes  of  the  individual  constituents  of  the  molecular 
species  in  question.  VECSELs  designed  for  the  mid-to-far  infrared  regime  are  especially  attractive  since 
many  of  the  chemical  and  biological  absorption  characteristics  exist  in  that  region. 

Additionally  in  Section  6.0,  we  report  related  work  that,  while  may  not  have  been  associated  with  the 
physical  principals  proposed  to  develop  tunable  high-brightness  semiconductor  sources,  do  affect  the  general 
performance  of  lasers,  or  can  be  used  in  concert  with  semiconductor  sources  for  the  purposes  of  enhancing 
desired  effects  for  specific  applications.  First,  in  Section  6.1  we  present  a  unique  method  to  produce  metallic 
nano-particles  that  on  a  semiconductor  surface  that  has  unique  polarization  characteristics.  We  then  introduce 
related  work  we’ve  pursued  in  order  to  obtain  very  low-void,  high  quality  solder  layers  in  Section  6.2.  In 
Section  6.3  we  present  the  work  of  an  AFRL/RY  funded  entrepreneurial  research  project  on  intracavity 
nonlinear  down-conversion  in  which  we  demonstrated  tunable  MWIR  lasers  with  a  cavity-enhanced 
nonlinear  conversion  efficiency.  Finally,  we  make  some  concluding  remarks  in  Section  7.0. 
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3.0  MID-WAVE  TYPE-I  QUANTUM  WELLS 


Short-wave  and  mid-wave  infrared  lasers  are  proving  to  have  many  applications  in  the  sensing  community, 
including  remote  chemical  sensing  of  numerous  organic  compounds  such  as  methane[l],  an  ideal  source  for 
low-altitude  laser  radar[2]  and,  if  sufficiently  bright,  infrared  countermeasures.  Although  there  are  several 
potentially  successful  approaches  to  developing  semiconductor  sources  in  the  short-wave  and  mid- wave 
infrared,  as  one  approaches  the  2-4  fi m  band,  challenges  exist  that  must  be  overcome  for  successful  adoption 
as  a  product. 

Within  the  last  15  years,  The  quantum  cascade  laser  (QCL)  has  proven  to  be  a  very  capable  semiconductor 
laser  platform,  owing  to  the  ability  to  rely  on  an  engineered  approach  to  intersubband  transitions,  making 
them  very  attractive  for  lower  energy  lasers.  However,  these  lasers  are  limited  by  band  offset  in  the 
conduction  band[7,  8],  and  the  desire  for  higher  energy  sources  is  progressing  toward  more  exotic  material 
system[9,  10].  Equally  as  promising  are  type-II  quantum  well  lasers  which  rely  on  electrons  and  holes  being 
confined  in  adjacent  layers [1 1].  When  these  layers  are  sufficiently  thin,  their  electronic  wavefunctions  have 
significant  overlap,  and  are  well  suited  for  wavelengths  above  ^2.8  fi m,  although  somewhat  restricted  at  the 
longer  wavelengths,  compared  to  QCLs  due  to  band  filling[12].  In  the  short-wave  and  mid- wave  infrared 
wavelengths,  their  low  Tq  necessitates  low  temperatures  to  achieve  appreciable  powers,  even  with  10  or  more 
quantum  wells  [11]. 

Type-I  quantum  well  lasers  can  address  the  2-4  /im  wavelength  band  with  the  gallium  antimonide  (GaSb) 
material  system.  Although  immature  relative  to  III-P  and  III-As  systems,  GaSb  provides  a  reasonably  flexible 
platform  for  the  (AlGaIn)(AsSb)  material  system,  and  allows  for  various  bandgaps,  refractive  index-,  and 
strain-profiles  within  the  active  region.  Ternary  InGaAs  quantum  wells[13]  are  suitable  for  A  <2.1  /im,  but 
are  generally  limited  to  low  compositions  of  indium  and  thin  layers  to  avoid  excessive  strain.  At  longer 
wavelengths,  it  becomes  necessary  to  utilize  the  InGaAsSb  quaternary  quantum  well  to  control  strain  and 
achieve  higher  indium  concentration [14].  Utilizing  a  “quinternary”  barrier  material  of  InAlGaAsSb,  room 
temperature  cw  operation  well  over  3  /im  has  been  demonstrated[15,  16].  Performance  of  this  material 
system  suffers  at  longer  wavelengths,  due  primarily  to  increased  Auger  recombination. 

We  chose  to  approach  the  problem  by  addressing  the  Auger  recombination  mechanisms  directly  and 
demonstrate  this  capability  in  the  region  of  A  =2  /im.  In  Section  3.1,  we  address  the  principles  of  Auger 
recombination  mitigation,  and  outline  our  approach.  In  Section  3.3,  we  describe  the  growth  and  fabrication 
of  broad  area  lasers  used  to  measure  the  Auger  carrier  losses,  followed  by  our  results  and  discussion  in 
Section  3.5. 


3.1  Type-I  Method  Auger  Engineering  Principals 

Type-I  quantum  well  lasers  are  attractive  at  wavelengths  below  about  2  /am,  but  at  longer  wavelengths,  they 
tend  to  suffer  from  increased  Auger  losses,  as  the  photon  energies  become  resonant  with  higher  subbands. 
This  effect  is  nominally  proportional  to  N^h,  where  Nth  is  the  threshold  carrier  density.  A  reduction  in  Nth 
has  been  the  most  common  approach  to  achieving  efficient  lasers.  Even  then,  it  has  been  shown  that  the 
proportionality  constant  (know  as  the  4  Auger  coefficient”,  and  typically  denoted  as  C)  dramatically  increases 
above  2  //m[17]. 

Unfortunately,  an  increased  differential  gain  (thus  practically  reducing  Nth)  is  often  coupled  with  increased 
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Auger  coefficient,  diminishing  the  returns  of  the  lower  operating  carrier  density.  The  carrier  density 
dependence  of  the  Auger  coefficient  is  not  commonly  considered,  which  results  in  incomplete,  and 
sometimes  incorrectly  designed  structures.  In  our  calculations,  we  avoid  this  empirical  relation,  and  rather 
calculate  the  Auger  processes  with  Simulase®,  a  product  of  Nonlinear  Control  Strategies,  Inc.  (NLCSTR), 
which  utilizes  many-body  theory,  determining  a  carrier-dependent  lifetime  (r)  which  we  use  to  evaluate 
quantum  well  design.  The  typical  approximations  used  for  Auger  scattering  equations  which  utilize 
simplified  Coulomb  coupling  matrix  elements  typically  predict  values  for  Auger  recombination  which 
deviate  from  actual  values  by  an  order  of  magnitude. 

For  the  purposes  here,  we  use  semiconductor  Bloch  equations,  quantum  Boltzmann  scattering  to  determine 
dephasing  of  P,  the  macroscopic  optical  polarization,  calculated  from  the  microscopic  polarizability,  as  well 
as  the  lineshape  of  absorption/gain[18].  Spontaneous  emission  is  determined  with  semiconductor 
luminescence  equations  with  higher  order  Coulomb  correlation  source  terms.  These  rigors  avoid  issues  with 
lineshape  errors  in  N  and  T  as  well  as  spectral  positions. 

Critical  for  narrow  bandgap  structures,  Auger  recombination  calculations  utilize  scattering  equations  in 
second  Born-Markov  approximation  with  fully  resolved  k-,  spin-,  and  subband- summations  as  well  as 
z-Integrals.  This  technique  is  more  thoroughly  discussed  in  Reference  19.  This  gain  model  requires  only 
input  defining  structural  layout  (layer  widths,  compositions),  and  basic  bulk  bandstructure  parameters  (e.g. 
Luttinger  parameters,  band  gaps,  and  band  offset[20]).  When  comparing  to  experiments,  the  only  free 
parameters  are  therefore  inhomogeneous  broadening,  Shockley-Read-Hall  (SRH)  recombination,  internal 
efficiency  (wrapped  into  pump  absorption)  and  internal  modal  optical  loss.  These  “free  parameters”  may  be 
measured  independently,  and  will  vary  depending  on  e.g.  interface  quality  and  compositional  control. 


We  then  apply  these  to  the  steady-state  values  of  a  simple  set  of  rate  equations.  Above  threshold,  these  are 
approximately: 


dN 

dt 

dNp 

dt 


+vg(Yg-a)  Np 


(1) 

(2) 


In  these  equations,  N  is  the  2-D  carrier  density,  rg  is  dominated  by  the  optical  absorption,  and  is  analogous  to 
internal  quantum  efficiency.  The  carrier  excitation  by  the  pump  is  denoted  as  J  ( J  =  <&/TkjJp),  and  M  is  the 
number  of  quantum  wells  in  the  structure,  assuming  generated  carriers  are  distributed  evenly  throughout  the 
wells.  The  SRH,  spontaneous,  and  Auger  lifetimes  are  given  by  ra^c,  respectively.  The  quantum  well  gain  is 
g ,  Np  is  the  2-D  photon  density,  rp  is  the  geometric  photon  lifetime,  T  is  the  3-D  mode  overlap  with  all  of  the 
quantum  wells,  and  a  is  the  waveguide  loss  (composed  of  everything  except  the  active  region  absorption).  In 
arriving  at  Equations  1  and  2,  we  ignore  carrier  exchange  between  the  barriers  and  wells  due  to  the  low 
barrier  carrier  density,  accounting  for  these  carriers  as  if  they  were  generated  directly  in  the  well.  We  also 
ignore  the  contribution  of  spontaneous  emission  to  photon  density.  Our  detection  scheme  did  not  have 
sufficient  sensitivity  to  discern  this  contribution  near  and  below  threshold,  and  spontaneous  emission  is 
negligible  above  threshold. 


There  are  many  types  of  Auger  recombination  that  can  occur  in  these  structures.  In  this  wavelength  regime, 
Auger  losses  are  dominated  by  processes  in  which  an  electron  and  hole  recombine  and  the  excess  energy  is 
used  to  excite  another  electron  high  into  a  conduction  subband  (CHCC  process)[21],  and  depend  strongly  on 
details  of  the  band  structure. 
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Figure  1.  Carrier  lifetime  dependence  on  conduction  band  offset. 

a)  Contributions  to  room  temperature  carrier  lifetime  as  a  function  of  AEC.  Spontaneous 
and  Auger  lifetimes  are  given  for  gth=600  cm -1.  Total  carrier  lifetimes  (including  an  SRH 
lifetime  of  4.5  ns)  are  given  for  gth  of 200,  400,  600,  800,  and  1000  cm~l.  The  total  carrier 
lifetime  associated  with  the  spontaneous  and  Auger  lifetimes  plotted  is  emphasized  with  a  shade 
underneath,  b)  transition  energies  between  the  lowest  conduction  band  state  and  excited  states. 
The  drop  of  Auger  lifetime  in  the  region  that  A  Ec  ~  Eg  indicates  the  increased  availability  of 
states  available  for  Auger  transitions. 


3.2  Engineering  Approaches  to  Auger  Processes 

Generally,  when  varying  the  barrier  potential,  A Ec,  the  dominant  factor  is  the  excess  photon  energy  (  e9) 

with  respect  to  the  barrier  level.  Commonly,  the  approach  is  to  design  for  a  large  A Ec,  improving  carrier 
capture  and  reducing  carrier  density.  The  reduced  carrier  density  is  intended  to  improve  the  Auger  lifetime 
through  the  empirical  relation: 

C=(rcN?h)~\  (3) 

where  C  is  a  proportionality  constant  known  as  the  “Auger  coefficient”.  The  excited  carrier  transitions  from  a 
localized  conduction  band  above  the  potential  barrier  to  a  bulk  state.However,  for  a  given  wavelength,  the 
Auger  lifetime  lifetime  can  vary  dramatically,  primarily  depending  on  the  E(J  relation  to  A Ec. 

This  relation  is  summarized  in  Figure  la,  where  we  compute  the  lifetimes  at  300  K  holding  the  gain  constant, 
adjusting  carrier  density  to  compensate  varied  differential  gains.  Carrier  lifetime  contributors  (spontaneous 
and  Auger  lifetimes)  are  plotted  as  a  function  of  barrier  height  in  empty-  and  filled-symbols  assuming  a 
threshold  of  600  cm-1,  as  will  be  seen  to  be  a  reasonable  approximation  for  two  of  our  devices.  Assuming  a 
SRH  lifetime  of  4.5  ns  (also  determined  through  fitting  to  experimental  data,  as  shown  in  Section  3.5),  the 
total  carrier  lifetime  is  also  shown  by  solid  black  lines.  These  are  calculated  for  200  cm-1- 1000  cm-1,  with 
the  data  corresponding  to  600  cm-1  specifically  identified  through  the  grey  fill  beneath  it.  Because  the  Auger 
recombination  is  sensitive  to  subband  positioning,  and  we  are  dominated  by  CHCC  processes,  conduction 
subband  positions  are  plotted  relative  to  the  lowest  electronic  subband  in  Figure  lb.  Finally,  for  reference,  we 
plot  the  photon  energy  in  Figure  lb,  corresponding  to  the  value  of  Eg  on  this  same  scale  to  show  where 
resonant  transitions  are  found,  causing  increased  Auger  recombination  wherein  excess  energy  created  from 
an  electron/hole  recombination  are  easily  transferred  from  E\  to  an  excited  electron  state  in  the  conduction 
band  Ex. 
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In  this  wavelength  regime,  all  practical  structures  are  found  in  the  region  where  Eg  »  A Ec.  At  low  values 
for  A Ec,  the  spontaneous  emission  lifetime  (roughly  proportional  to  inverse  threshold  carrier  density),  is 
long.  Moreover,  for  Auger  recombination  in  this  region,  the  final  electron  is  excited  into  states  that  are 
delocalized  in  the  barrier,  reasonably  approximated  by  bulk-states,  and  a  transition  from  a  confined  subband 
to  this  delocalized  band  requires  a  large  transfer  of  momentum[22].  Because  the  Coulomb  matrix  elements 
that  mediate  the  Auger  transitions  strongly  decrease  with  a  change  in  momentum,  the  Auger  transition  rate 
decreases  (lifetime  increases)  with  increasing  excess  energy  above  the  barrier  level. 

At  higher  barriers  (A Ec  >  0.35  eV),  Auger  recombination  dominates  the  carrier  lifetime  as  the  Eg  ~  A Ec. 
This  is  due,  in  part  to  carrier  excitations  to  more  localized  wavefunctions  resulting  in  reduced  momentum 
transfer.  Within  this  region,  coupling  is  very  dependent  on  band-structure.  Of  note  is  the  region  where 
Eq  —  Ei  and  E 7  —  E\  are  close  to  Eg,  which  promotes  resonant  transitions;  the  transition  probability 
between  these  excited  electronic  states  increases  dramatically,  causing  the  Auger  lifetime  to  drop  to  around 
10  ps.  The  effects  can  also  be  seen  where  A Ec  >1.2  eV  in  Figure  1.  In  this  region,  E$  —  E\  is  sufficiently 
resonant  with  Eg,  reducing  rc  by  a  factor  of  2. 

For  values  x  >0.45  in  Al^Gai-^Sb,  the  X  valley  drops  below  the  T  valley,  causing  the  material  to  become 
an  indirect  bandgap[20].  This  condition  is  generally  avoided  near  the  quantum- well,  especially  in  electrically 
injected  lasers,  and  where  optically  pumped  lasers  absorb  significant  pump  in  the  barrier  regime,  typically  to 
increase  net  absorption.  However,  we  may  computationally  explore  the  region  of  strong  electronic 
confinement,  where  the  CHCC  coupling  is  between  two  quantized  states  (where  Eg  <  A Ec),  also  addressed 
in  Reference  22.  In  this  case,  the  Auger  lifetime  is  smallest  when  a  subband  is  resonant  with  the  lowest 
electronic  subband,  with  a  energy  difference  of  ~  Eg,  e.g.  without  requiring  a  transfer  of  in-plane 
momentum.  As  such,  this  relation  is  strongly  dependent  on  the  specifics  of  the  structure,  and  generalities  are 
significantly  more  difficult  to  make. 

For  comparison  to  the  common  C N3  relation  for  the  Auger  recombination,  Figure  2  shows  the  Auger 
lifetime  (repeated  from  Figure  1)  and  how  that  relates  to  the  analytical  Auger  coefficient  C.  The  inverse 
relation  dominates  the  general  trend  because  Nt}-L  is  relatively  constant  for  most  values  of  A Ec.  At  very  low 
values  for  A Ec,  Nth  increases  dramatically  due  to  the  poor  carrier  confinement,  and  this  relation  deviates 
significantly.  We  also  show  the  approximation  from  Reference  22  for  the  CHCC  coefficient  for  the  case 
where  Eg  >>  A Ec.  It  should  be  noted  that  the  term  Vc  is  used  in  place  of  A Ec  in  Reference  22. 

The  three  structures  all  represent  a  region  where  Eg  >  A Ec,  this  ratio  being  1.32,  1.62,  and  3.35  for  HB, 
MB,  and  LB,  respectively.  Only  the  LB  structure  is  well  approximated  by  the  coupling  from  the  discrete  to 
the  continuous  spectrum  in  the  conduction  band[22]. 

3.3  Growth  and  Fabrication 

3.4  Method  and  Characterization 

An  enabling  tool  for  producing  mid-infrared  laser  heterostructures  as  complicated  as 
In^Gai-^Asi-ySby/GaSb  semiconductor  lasers  is  molecular  beam  epitaxy  (MBE),  whereby  monolayers  of 
single  crystals  can  be  synthesized  with  atomic  layer  precision.  A  Gen-II  molecular  beam  epitaxy  reactor  (in 
support  of  this  research  effort)  was  successfully  transferred  from  Wright-Patterson  Air  Force  Base  (WPAFB) 
to  University  of  Missouri,  Columbia  (UM).  It  was  employed  for  the  production  of  As-based  quantum  cascade 
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Figure  2.  Auger  lifetime  and  coefficients  calculated  at  room  temperature. 

Comparison  of  calculated  Auger  lifetime  v.v.  Auger  coefficient,  where  the  latter  (C)  is  taken  to 
be  equivalent  to  (rcNfh)  \  C\for  the  CHCC process  as  given  in  Reference  22. 
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Figure  3.  Energy  diagrams  of  the  three  structures  HB,  MB,  and  LB  structures. 

Conduction  band  is  marked  “CB”,  the  heavy-hole  valence  band  is  labeled  “VB(HH)”,  while 
the  light-hole  valence  band  is  marked  “VB(LH)”.  The  calculated  electron  energy  bands  (Ecn) 
and  hole  energy  bands  ( Evn)  are  also  plotted  for  each  structure. 
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Figure  4.  RHEED  schematic 

Staib  Instruments  RHEED  system  with  a  12  bit  digital  camera  and  K- space  400  analytical 
system. 


lasers  and  type-I  vertical  external  cavity  surface  emitting  laser  (VECSEL)  devices.  Device  structures  were 
synthesized  via  molecular  beam  epitaxy  (MBE)  using  detailed  band  structure  engineering  for  the  purpose  of 
addressing  optical  power  limitations  in  mid-infrared  devices.  Various  VCSEL  structures  with  various  barrier 
heights  were  simulated,  fabricated  and  analyzed. 

In^Gai-a; Asi_ySby/GaSb  structures  are  produced  using  MBE,  where  several  ovens  of  ultrapure  materials 
were  heated  to  produce  an  evaporated  beam  of  molecules.  The  MBE  process  conditions  necessary  for 
producing  high-quality  devices  include  an  ultra-high  vacuum  environment  (10-10  Torr)  yielding  material 
with  impurity  levels  below  ten  parts  per  billion,  precise  control  and  reproducibility  of  composition  and 
doping,  and  high-quality  layers  with  very  abrupt  interfaces.  At  the  facilities  of  one  of  our  collaborators 
(University  of  Missouri),  the  MBE  system  consists  of  three  separate,  but  connected  vacuum  chambers  (i.e.a 
growth  chamber,  a  buffer  chamber  and  an  intro  chamber)  with  ion-  and  cryo-pumps  to  obtain  requisite 
vacuum  levels.  The  growth  chamber  also  included  a  titanium  sublimation  pump  (TSP),  which  was  used  to 
further  enhance  the  vacuum  environment. 

The  intro  chamber  is  primarily  for  housing  MBE-grown  samples  and  for  loading  new  substrates,  while 
vacuum  was  maintained  in  the  other  two  chambers  via  manual  gate  valves.  A  continuous  supply  of  liquid 
nitrogen  throughout  the  cryoshroud  and  substrate  heater  [along  with  university-supplied  chilled  water  around 
the  Knudsen  cells  (k-cells)  and  through  bulk  evaporator  and  cracker  containing  cells]  are  used  to  provide 
cooling,  minimize  outgassing  of  the  vacuum  chamber  and  minimize  thermal  contributions  of  adjacent  cells. 
To  maintain  a  specific  process  condition  for  each  experiment,  a  cluster  of  proportional-integral-derivative 
(PID)  closed  loop  control  systems  that  monitored  temperature  zones  within  each  oven  are  integrated  along 
with  computer  controlled  (pneumatic)  shutters  in  front  of  each  cell  for  precise  control  over  the  deposition 
rate.  The  MBE  system  contains  arsenic,  antimony,  gallium,  indium,  aluminum,  silicon  (doping),  GaTe 
(doping),  and  beryllium  (doping). 

For  real-time,  in-situ  process  monitoring,  a  Staib  Instruments  reflection  high-energy  electron  diffraction 
(RHEED)  system  with  a  6-in  phosphor  screen  is  used  to  analyze  all  device  structures  during  growth.  The  use 
of  RHEED  to  monitor  the  surface  during  changes  in  substrate  temperatures,  growth  rates,  and  flux  ratios  is 
vital  to  the  success  of  achieving  good  surface  morphology.  Apart  from  surface  morphology  characteristics, 
RHEED  is  also  used  in  determining  the  substrate  surface  conditions.  Figure  4  shows  a  typical  setup  and 
interaction  of  RHEED  with  the  substrate  surface  and  recorded  using  a  camera. 
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After  loading,  samples  are  transferred  into  the  growth  chamber  and  ready  for  growth  after  removal  of  the 
protective  oxide  layer.  The  laser  structures  in  this  research  effort  are  grown  on  epi-ready  (100)  GaSb 
substrates  from  WaferTech.  Beam  equivalent  pressures  at  possible  growth  temperatures  for  all  k-cells  cells 
are  periodically  monitoring  using  an  ion  gauge  lowered  directly  in  the  path  of  where  the  substrate  would  be 
positioned  during  growth.  These  measurements  are  also  used  to  calculate  the  activation  energy  of  each  k-cell. 
The  RHEED  patterns  are  used  to  monitor  the  surface  morphology  and  growth  kinetics  at  specific 
temperatures. 

In^Gai-^Asi-ySby  remains  an  important  semiconductor  material  for  mid-infrared  photonic  devices. 
Conduction  and  valence  band  alignment  of  In^Gai-^Asi-ySb^  quantum  wells  with  other  zinc-blende 
compounds  (aluminum  antimonide  (AlSb),  GaSb,  indium  arsenide  (InAs))  provide  advantageous 
performance  enhancements  to  higher  efficiency  gain  mediums  in  ^3  /iir 1  sources.  Unfortunately,  the  ability 
to  synthesize  In^Gai-^ Asi_ySby  pseudomorphically  across  the  compositional  range  has  been  limited  by 
substrate  availability.  Lattice  mismatch  with  GaSb  has  prohibited  the  exploitation  of  the  entire  wavelength 
tuning  range  that  InxGai_xAsi_ySby  could  potentially  offer.  There  has  been  considerable  effort  to  exploit 
high-indium  content  In^Gai-^ Asi_ySby  for  longer  wavelength  operation,  yet  high  misfit  dislocation 
densities  are  inevitable  and  the  miscibility  gap  is  a  formidable  barrier  (see  Figure  5b).  To  date,  the  ability  to 
cover  the  compositional  range  of  In^Gai-^Asl  —  ySby  as  well  as  its  range  of  properties  remains  unresolved. 
In  this  effort,  we  explore  In^Gai-^Asi-ySby/GaSb  structures,  where  the  measured  indium  mole  fraction  ( x ) 
varies  from  x=0  to  x>0.50.  In  addition  to  atomically  smooth  structures  with  growth  two-dimensional  growth 
fronts  Figure  5a,  we  observed  three-dimensional  networks  of  quantum  dashes  and  results  reveal  a 
self-organized  composition  modulation.  Some  physical  features  of  the  quantum  dashes  include  near 
one-micron  lengths,  90-deg  flip  in  orientation,  and  uniformity  across  a  20x20  /m l  area.  Samples  of 
MBE-grown  mid-infrared  VCSEL  structures  under  various  process  conditions  are  presented  in  Figure  6  to 
demonstrate  the  influence  on  the  process  conditions  on  wavelength  tunability.  While  these  devices 
demonstrate  lower  indium  concentration  than  planned  for  the  MWIR  structures,  the  principals  are  identical 
and  we’re  using  the  concepts  transitioning  to  MWIR  as  well. 


Figure  5.  Planview  AFM  of  InGaAsSb  surface  with  varying  concentrations 

a)  Ino.sGao^Sb,  b)  In^Ga^iSb,  c)  In^^Ga^^As^^Sb^^ 


10 

Approved  for  public  release;  distribution  unlimited. 


Figure  6.  Reflectivity  and  photoluminescence  of  several  samples 

a)  InGaSb  VECSEL  structure  (with  DBR),  b)  InGaSb  VECSEL  structure  with  additional  DBRs 
and  modified  growth  conditions ,  and  c)  InGaSb  VECSEL  structure  ( modified  growth  conditions). 

3.4.1  Test  Structures 

Prior  to  laser  growth  and  to  ensure  precise  alloy  composition  and  thickness  for  each  laser  design  presented, 
deposition  rates  are  varied,  measured  in  situ  by  RHEED,  and  confirmed  post-growth  using  HR-XRD. 
Samples  are  also  analyzed  by  atomic  force  microscopy  (AFM),  which  demonstrates  that  these  structures 
were  indeed  produced  via  step  flow  growth. 

In^Gai-^Sb  structures  are  also  produced  to  calibrate  indium  composition  in  the  quantum  well  and  to  ensure 
2  pm  emission,  which  is  confirmed  with  photoluminescence.  For  each  laser  design,  the  MBE  process 
conditions  (flux  ratio,  substrate  temperature,  and  deposition  rate)  are  selected  to  ensure  quality  and 
composition  of  each  binary  and  ternary  layer.  The  shutter  sequence  and  antimony  valve-positioner  are 
computer  controlled,  thus  maintaining  run-to-run  consistency.  Each  quantum  well/barrier  combination  is 
embedded  in  an  AlGaSb  waveguide  that  is  adjusted  to  maintain  a  pump  absorption  of  approximately  85%, 
and  a  transverse  mode  overlap  of  about  5%.  In  order  to  accomplish  this,  HB  and  MB  structures  have  3 
repeats  of  the  quantum  wells,  each  separated  by  100  nm  to  assure  electronic  isolation.  To  achieve  the  same 
overlap  and  absorption,  while  maintaining  single  transverse  mode  operation,  the  LB  structure  necessitates 
only  two  quantum  wells. 

Strained  In^Gai-^Sb/Al^Gai-ySb  quantum  well  lasers  are  produced  by  solid-source  MBE  system. 
Semi-insulating  GaSb  substrates  are  mounted  in  indium-free  molybdenum  blocks  with  diffuser  plates 
(located  between  the  substrate  and  the  substrate  heater)  to  achieve  improved  temperature  uniformity  across 
the  wafer.  During  substrate  outgassing,  deoxidation  temperatures  are  monitored  using  a  10  keV  RHEED 
system  and  are  confirmed  by  a  “1x3”  reconstruction  pattern.  The  growth  temperatures  are  relative  to  the 
GaSb  deoxidation  temperature  and  also  referenced  by  the  “1x5”  phase  transition  while  under  antimony 
over-pressure. 

Each  quantum  well/barrier  combination  is  embedded  in  an  AlGaSb  waveguide  that  is  adjusted  to  maintain  a 
pump  absorption  of  approximately  85  %,  and  a  transverse  mode  overlap  of  about  5  %.  In  order  to  accomplish 
this,  HB  and  MB  structures  have  3  repeats  of  the  quantum  wells,  each  separated  by  100  nm  to  assure 
electronic  isolation.  To  achieve  the  same  overlap  and  absorption,  while  maintaining  single  transverse  mode 
operation,  the  LB  structure  necessitates  only  two  quantum  wells. 
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0  (mrad) 

Figure  7.  Example  of  HR-XRD  £2-2©  rocking  curves. 

This  data  represents  AlSb/GaSb  superlattices  along  with  calculated  diffraction  pattern,  rep¬ 
resenting  faithful  reproduction  material  composition  and  layer  thickness.  Inset  is  as  example 
surface  profile  measurement  by  AFM. 


As  the  gain-guided  structures  are  defined  by  the  optical  pump,  no  lithographic  steps  are  necessary.  The  bare 
wafers  are  thinned  to  nominally  80  /im  using  mechanical  lapping  and  then  cleaved  into  bars  of  various 
lengths  (375  /im  to  1500  /im).  In  this  report,  only  1500  /im  length  cavities  are  used. 

3.5  Type-I  Results  and  Discussion 

Laser  bars  are  mounted  epitaxial- side  up  using  thermal  paste  to  a  low-oxygen  copper  heat  sink  temperature 
controlled  using  a  two-stage  thermoelectric  cooler,  allowing  various  temperatures  at,  and  around  room 
temperature  (from  about  265  K  to  nominally  325  K).  We  maintain  a  slow  dry  nitrogen  purge  to  avoid 
condensation  at  temperatures  below  the  dew-point.  The  lasers  are  pumped  using  a  980  nm  semiconductor 
laser  bar  with  a  spot  125  p m  wide  and  3  mm  long.  Pulses  are  2  /is  long  with  a  1  kHz  pulse  repetition 
frequency,  avoiding  thermal  loading. 

A  fraction  of  the  optical  power  is  collected  using  a  fiber-coupled  extended-range  InGaAs  detector.  Both  the 
pump  laser  and  device  powers  are  time-gated  to  remove  pulse  transients.  The  pump  signal  is  calibrated  and 
compared  to  CW  operation,  integrating  over  multiple  cycles,  to  determine  the  appropriate  power  vs.  detector 
voltage.  While  the  device  power  remains  uncalibrated,  lasing  threshold  is  the  most  important  metric  to 
evaluate  carrier  losses  and  therefore  device  powers  remain  in  arbitrary  units,  although  the  noise  floor  of  this 
detector  is  ^0.01  on  this  arbitrary  scale.  Finally,  we  also  measure  spectra  at  each  temperature  and  pump  level 
to  confirm  the  device  is  behaving  as  we  expect. 

We  measure  the  power  characteristics  from  approximately  350  lasers  spanning  the  three  growth  designs 
mentioned  in  Section  3.1,  and  report  these  in  Figure  8.  The  mean  data  from  each  of  the  three  structures  are 
indicated  by  the  data  points,  while  the  error  bars  signify  standard  deviation  from  the  collected  data. 
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Figure  8.  Laser  power  characteristics  comparing  devices  with  various  barrier  heights. 

Laser  power  characteristics  given  for  the  three  designs  outlined  in  Section  3.1,  as  a  function 
of  equivalent  pump  current  density.  Data  is  shown  on  a  logarithmic  plot  to  emphasize  the 
threshold. 

Although  we  pulse  the  pump  to  avoid  thermal  loading,  we  ignore  transients  allowing  us  to  set  the  left  sides  of 
Equations  1  and  2  to  zero  to  find  the  steady  state.  In  this  case,  we  find  the  threshold  condition  of  the  laser: 

(4) 

M  T  Ta  n  Tc 

where  M  is  the  number  of  wells,  and  r  represents  the  combination  of  all  three  carrier  loss  terms,  and  is 
commonly  referred  to  as  the  “carrier  lifetime.”  Figure  8  shows  the  power  characteristics  of  this  laser.  The 
absorbed  pump  fluence  per  well  is  equivalent  to  rjiJhujp/M  in  Equation  1,  and  is,  of  course,  the  driving  term 
of  the  laser.  It  is  immediately  evident  that  the  HB  structure  has  the  highest  threshold  per  well,  a  result  that  is 
not  unpredictable  due  to  the  large  Auger  recombination  as  discussed  in  Section  3.1  and  calculated  in  Figure  1. 

From  Figure  8,  we  see  the  threshold  per  well  is  480,  283,  and  255  A/cm2/well,  which  corresponds  to 
equivalent  pump  fluences  of  606,  358  and  255  W/cm2/well  for  the  HB,  MB,  and  LB  structures,  respectively. 
Because  rn  is  used  as  a  fit  parameter,  we  only  find  relative  changes  to  the  carrier  lifetimes.  Experimental 
results  compare  favorably  to  theory  in  the  MB  and  HB  structures  for  an  internal  absorption  oli  =10  cm-1  and 
a  SRH  lifetime  ra  =4.5  ns.  In  contrast,  the  LB  structure  had  a  larger  internal  loss  of  (cti)  =20  cm-1,  but  the 
SRH  lifetime  was  larger  as  well  (20  ns),  offsetting  some  of  this  additional  loss.  This  is  not  altogether 
surprising  because  the  LB  structure,  with  only  two  wells,  contains  fewer  interfaces  within  the  high  carrier 
density  region,  where  SRH  recombination  typically  presents  itself.  However,  the  transverse  mode  profile 
contains  longer  tails  and  therefore  could  incur  higher  losses  due  to  optical  scattering. 

From  these  parameters,  we  can  determine  the  threshold  carrier  densities  (Nth)  are  1.15x  1012  cm-2, 

1.67  x  1012  cm  2,  and  2.12x  1012  cm  2  for  the  HB,  MB,  and  LB  structures,  respectively.  From  these,  we  can 
compute  the  carrier  lifetimes  for  the  MB  and  LB  structures  as  2.46  x  and  6.56  x  that  of  the  HB  “standard” 
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temperature  (K)  temperature  (K)  temperature  (K) 

Figure  9.  Measured  and  calculated  threshold  for  various  temperatures. 

Threshold  (in  equivalent  current )  and  T0  and  contributions  from  the  three  carrier  loss  contribu¬ 
tors  Ja,b,cfor  a)  the  HB  structure ,  b)  the  MB  structure ,  and  c)  the  LB  structure. 


structure.  This  represents  a  significant  (>6x)  improvement  of  the  carrier  lifetime,  resulting  in  a  significantly 
lower  threshold  density. 

A  more  definitive  test  of  the  model  and  the  device  improvement  is  the  temperature  dependance  of  the 
threshold,  and  therefore  loss  mechanisms.  This  dependence  is  commonly  denoted  as  Tq  and  empirically 
behaves  as  an  exponential  growth[21].  Although  this  relation  is  not  completely  accurate[19],  it  is  often  a 
suitable  fit  for  a  wide  variety  of  temperatures. 

Figure  9  shows  the  temperature-dependance  of  the  HB  (Figure  9a),  MB  (Figure  9b),  and  LB  (Figure  9c) 
structures.  These  curves,  while  including  the  rji  term,  are  not  normalized  by  the  number  of  wells  as  in 
Figure  8.  We  also  show  in  Figure  9  the  contributing  terms  to  the  total  threshold  fluence.  In  each  case,  the 
relative  contribution  of  the  SRH,  spontaneous,  and  Auger  are  denoted  by  Ja,b,c,  respectively,  and  are 
computed  as  a  function  of  temperature.  The  same  terms  for  (c^)  and  ra  are  used  for  this  graph  as  those  used 
to  compare  in  Figure  8. 

We  can  see  in  each  of  these  examples,  the  measured  Tq  is  closely  predicted  by  the  model.  In  the  best  case, 
these  deviated  by  3  K  for  the  HB  structure.  We  expect  the  kink  at  about  300  K  should  have  been  more  evident. 
This  kink  is  due  to  increased  Coulomb  matrix  elements  as  Eg  drops  at  higher  temperatures  approaching  A Ec, 
where  less  momentum  is  required  for  the  transition.  In  the  worst  case,  this  deviation  was  as  high  as  15  K  for 
the  LB  structure,  although  this  effect  could  be  attributed  to  carrier  transport  (e.g.thermionic  emission). 
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For  each  of  these  cases,  we  can  compare  the  fraction  of  Auger  recombination  to  that  of  the  standard  HB 
structure.  Indeed,  despite  the  increased  threshold  carrier  density,  Jc  is  reduced  from  1.13  kA/cm2  to 
0.25  kA/cm2,  a  4-fold  decrease  in  carrier  losses  due  to  Auger  recombination.  However,  due  to  the  threshold 
carrier  density  disparity,  this  is  actually  a  5.7  x  increase  in  Auger  recombination  lifetime,  reflecting  a 
significant  improvement  to  Auger  losses. 

As  stated  previously,  the  metric  is  the  quasi-empirical  Auger  coefficient  (C),  related  to  rc  by  Equation  3. 
Although  this  relation  is  not  strictly  held,  we  can  calculate  the  value  of  C  relative  to  the  HB  structure  for 
reference;  this  value  is  found  to  be  6.4  x  lower  for  the  MB  structure,  and  a  significant  19.2x  lower  for  the  LB 
structure.  We  show  then  that  C,  most  often  considered  constant  for  a  given  bandgap,  to  be  easily  varied  by 
20  x,  and  we  can  widen  this  more  dramatically  with  stronger  electron  confinement  (as  high  as  ^0.55  eV  with 
barriers  composed  of  Alo.45Gao.55Sb),  although  this  would  not  improve  the  device  any  more,  but  rather  make 
the  lower  performing  devices  less  suitable. 

We  therefore  show  this  concept  for  directly  engineering  the  Auger  recombination  rate  is  supported  by 
experimental  results.  This  approach  and  methodology  is  not  only  appropriate  at  2  /4m,  but  can  be  expanded 
to  longer  wavelengths.  Such  a  numerical  study  is  outlined  in  Section  3.7. 


3.6  IMF-grown  Materials 

One  of  the  promising  techniques  for  MWIR  lasers  was  the  use  of  interfacial  misfit  (IMF)  to  grow  Ill-Sb 
low-energy  active  regions  on  III- As  high  thermal  conductive  DBRs.  Fig.  1(a)  shows  a  schematic  illustration 
of  the  fabricated  VECSEL  structure  on  a  gallium  arsenide  (GaAs)  substrate.  The  bottom  DBR  is  comprised 
of  25  pairs  of  Alo.95Gao.05As/GaAs.  The  active-region  consists  of  9  In0.2Ga0.sSb  QW  separated  by  pump 
absorbing  Alo.25Gao.75Sb  barriers.  An  AlSb  confinement  layer  on  each  side  of  the  active  region  ensures 
carrier  confinement  in  the  structure,  and  a  thin  GaSb  capping  layer  is  added  on  top  to  prevent  oxidation.  A 
low-resolution  TEM  image  of  the  structure  is  shown  in  Figure  10a,  the  structure  matches  the  epitaxially 
grown  structure  (only  a  few  DBR-pairs  are  captured  in  image).  Figure  10b  also  shows  the  ^ -field  profile  for 
the  lattice-matched  VECSEL.  The  key  feature  of  the  antimonide  VECSEL  based  on  the  GaAs/AlGaAs  DBRs 
is  the  ability  to  achieve  lasing  despite  the  7.78%  mismatched  interface  between  the  antimonide  active  region 
and  the  arsenide  based  DBR.  The  first  quantum  well  is  located  within  500  nm  of  this  highly  mismatched 
interface.  Thus  for  such  a  VECSEL  to  work  the  growth  has  to  be  optimized  to  achieve  minimal  threading 
dislocations  in  the  active  region  and  to  achieve  a  reproducible  relaxation  of  the  metamorphic  buffer. 


3.6.1  Limitations  of  the  IMF  for  SWIR  active  regions 

Unfortunately,  what  was  found  was  that  when  compared  to  the  same  active  region  on  lattice-matched  GaSb 
DBR,  the  structure  grown  on  IMF  had  a  much  higher  threshold  power  density  than  that  of  the  lattice-matched 
VECSEL  (nominally  13  x  higher)  and  the  slope  efficiency  is  approximately  1/4  that  of  the  structure  grown  on 
Ill-Sb.  It  can  be  speculated  that  the  deterioration  in  the  performance  of  the  VECSEL  when  grown 
mismatched  on  GaAs  compared  to  when  it  is  grown  lattice  matched  on  GaSb  can  be  attributed  to  the 
presence  of  threading  dislocations  in  the  active  region  when  it  is  grown  on  the  GaAs/AlGaAs  DBR.  In 
Equation  1,  the  second  term  on  the  right-hand  side,  taken  alone,  represents  the  resulting  decay  response  when 
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Figure  10.  Layer  structure  and  image  of  IMF  structure 

a)  Schematic ,  b)  cross-section  TEM  image  and  plot  of  \E\2 -field  for  lattice-mismatched  2  pm 
Ill-Sb  VECSEL  grown  on  GaAs/AlGaAs  DBR. 


carriers  are  excited.  If  the  lifetimes  (ra^c)  are  constant,  an  exponential  results1.  Through  the  spontaneous 
emission  term  of  the  form  /3spT 75  (where  f3sp  is  the  rate  at  which  spontaneous  photons  are  coupled  to  the 
optical  cavity  modes,  and  T  is  the  modal  overlap),  this  carrier  decay  is  imprinted  on  the  photon  generation. 
As  is  custom  in  such  laser  systems,  this  term  is  neglected  in  Equation  2,  an  equation  which  is  appropriate 
only  above  threshold  where  the  spontaneous  emission  contributes  insignificantly  to  the  photon  generation. 

It  is  a  reasonable  to  approximate  the  carrier  decay  terms  by  observing  the  time-resolved  photoluminescence 
signal.  Figure  11a  shows  several  examples  of  the  spectrally  integrated  PL  signal  as  a  function  of  time,  after 
excitation  time  (occurring  at  approximately  2  ns).  The  two  sets  of  data  are  displayed  for  the  average 
excitation  power  densities  of  30  and  3000  W/cm2.  The  measured  transient  photolumiscence  (PL)  intensities 
were  generally  well  fitted  with  the  double  exponential  function  of  I (t)  =  I  dark  +  A\e  t^Tl  +  A2C  t/T 2, 
where  Idark  is  the  dark  counts  without  the  pulsed  laser  excitation,  A\  and  A2  are  the  amplitudes  of  each  term, 
and  t\  and  72  are  the  decay  time  constants.  These  parameters  were  determined  by  the  fittings  shown  by  the 
solid  lines  in  Figure  11a  and  are  summarized  in  Table  1. 

To  identify  the  origin  of  the  two  decay  time  constants  of  t\  and  72,  we  have  studied  the  time  decay  by 
spectral  integrations  in  limited  ranges.  For  example,  the  PL  decay  measured  on  the  quantum  well  with  IMF 
was  spectrally  separated  into  the  two  regions  of  1367-1405  nm  and  1406-1444  nm  and  integrated  in  each 
region.  This  spectral  separation  resulted  in  single  exponential  decays  and  the  decay  time  constants 
determined  in  the  lower  and  higher  photon  energy  regions  very  nicely  reproduced  the  t\  and  72  values  given 

1  As  discussed  previously,  these  carrier  lifetime  contributions  are  considered  constant  only  in  the  small  carrier  density 
modulation  scheme  -  that  is,  they,  themselves,  are  carrier  density  dependent,  causing  a  “multi- exponential”  decay.  This 
complex  decay  can  aid  in  the  distinction  between  multiple  non-radiative  effects. 
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Figure  11.  Time-resolved  photoluminescence  and  STEM  images  of  structure 

a)  Time-resolved  photoluminescence  measured  on  QW  grown  with  IMF  with  the  average 
excitation  poiwers  of30W/cm2  and  3000W/cm2,  b)  STEM  image  of  lattice  mismatched  Ill-Sb 
VECSEL  on  GaAs/AlGaAs  showing  the  source  of  the  threading  dislocations  at  the  IMF  array 
and  subsequent  reduction  in  the  density  as  a  function  of  distance  from  the  mismatched  interface 


Table  1.  IMF  PL  Decay  Characterization 

Time-resolved  photoluminescence  characterization 


Excitation 

W/cm2 

QW  with  IMF 

M/M  n  (ps)  r2  (ps) 

QW  without  IMF 

M/M  T\  (ps)  r2  (ps) 

30 

4.47 

1006 

222 

16.78 

1560 

250 

150 

3.46 

1000 

260 

6.44 

1530 

243 

750 

3.55 

1050 

310 

6.23 

1450 

248 

1550 

2.98 

1004 

330 

6.79 

1550 

258 

3000 

2.69 

1019 

351 

6.58 

1590 

257 
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in  Table  1,  respectively.  Therefore  r\  and  72  are  the  decay  time  constants  related  to  the  quantum  well  ground 
states  and  the  excited  states,  respectively.  The  comparison  of  the  two  sets  of  data  in  Table  1  on  the  quantum 
wells  with  and  without  IMF  shows  a  clear  difference  of  the  excitation  power  dependence  of  the  r 2  value. 
With  IMF  the  72  value  and  the  amplitude  A2  increased  with  the  excitation  power,  while  they  were 
independent  of  the  excitation  power  without  IMF.  After  the  pulsed  laser  excitation  of  the  quantum  well 
samples,  the  excited  electrons  and  holes  experience  energy  relaxation.  Therefore  the  72  value  of  the  excited 
states  includes  this  energy  relaxation  process  and  is  much  shorter  than  the  t\  value  of  the  ground  states.  The 
increase  of  the  excitation  power  generally  enhances  the  band  filling  of  the  quantum  well  ground  states  with 
the  photo-generated  carriers,  and  this  reduces  the  energy  relaxation  rate  in  both  the  conduction  and  valence 
bands  due  to  the  filling  of  the  energy  states  above  the  band  edges.  This  will  result  in  the  increase  of  the  r 2 
value  for  the  quantum  wells  with  IMF.  The  enhanced  steep  decay  in  Figure  1  la  for  the  higher  excitation 
corresponds  to  the  increase  of  the  amplitude  A2.  The  fact  that  the  quantum  wells  without  IMF  did  not  show 
any  excitation  power  dependence  indicates  the  presence  of  residual  defects  that  prevent  the  band  filling  effect 
in  the  quantum  well  ground  states. 

This  band  filling  effect  is  consistent  with  threading  dislocations  originating  from  the  IMF.  Figure  lib  shows 
an  transmission  electron  microscopy  (TEM)  micrograph  of  the  Ill-Sb  VECSEL  that  demonstrates  a 
significant  density  of  dislocations,  which  is  reduced  farther  away  from  the  IMF  interface  (at  the  GaSb/GaAs 
junction).  Carrier  lifetimes  due  to  defect  recombination  can  be  approximated  via  techniques  described  more 
in  depth  in  References  23  and  24: 


1  _  iikT  l 

T/  ’  ^ 

T dislocation  Q  v 

where  /i  is  the  carrier  mobility,  q  is  the  electron  charge,  l  is  the  projected  dislocation  length  (estimate),  V  is 
the  volume  being  analyzed.  In  this  way,  we  estimate  the  dislocation  density  Np  as  a  function  of  distance 
from  the  GaSb/GaAs  interface  (this  information  is  presented  in  Figure  lib).  The  dislocation  density  reduces 
from  ^108  cm-2  to  about  107  cm-2  at  a  distance  of  3  /im  from  the  interface.  Thus  it  is  evident  from  the 
relationship  between  carrier  lifetime  and  dislocation  density  that  the  higher  threading  dislocation  density 
found  in  the  Ill-Sb  VECSELs  grown  lattice  mismatched  on  GaAs/AlGaAs  DBRs  is  detrimental  to  the 
VECSEL  performance.  Furthermore  this  high  threading  dislocation  density  in  the  Ill-Sb  on  GaAs/AlGaAs 
can  also  be  linked  to  the  high  threshold  density,  lower  efficiency  and  reduced  maximum  output  power 
compared  to  a  lattice  matched  Ill-Sb  VECSEL. 

3.7  Longer  wavelengths 

Following  on  our  success  with  devices  operating  at  A  ^2  /jl m[25,  26],  we  have  proceeded  to  structures 
designed  to  operate  at  3.5  /im,  which  is  much  longer  than  type-I  lasers  are  typically  shown  to  operate.  In 
order  to  do  this,  we  must  rely  on  quaternery  material  (InGaAsSb  quantum  wells),  and  compensate  the  strain 
with  quaternary  barriers  (AlGaAsSb).  Figure  2  shows  a  device  that  is  nominally  similar  to  Reference  14, 
where  the  14  nm  wide  quantum  well  is  made  of  Ino.5Gao.5Aso.19Sbo.8i,  flanked  by  Al^Gai-^As^Sbi-y 
barriers.  In  Figure  2a,  the  aluminum  concentration  in  the  barrier  material  is  altered  to  adjust  the  A Ec,  while 
the  arsenic  is  adjusted  to  maintain  a  constant  strain  throughout  the  simulated  structures. 

All  practical  devices  for  such  a  structure  are  found  for  0.4  <  A Ec  <0.82  eV.  The  lower  bounds  (^0.4  eV)  is 
where  barrier  material  is  GaSb.  At  the  upper  bounds  (for  values  x  >0.45  in  Al^Gai-^AsSb),  the  X  valley 
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Figure  12.  Simulated  Auger  and  spontaneous  lifetimes  for  ~3.5  /xm  designed  quantum 
wells. 

a)  Varied  A Ec  and  constant  strain ,  and  b)  at  the  largest  Auger  lifetime  of  (a),  varying  strain 
through  arsenic  content. 


drops  below  the  Y  valley,  causing  the  material  to  become  an  indirect  bandgap[20].  This  condition  is  generally 
avoided  near  the  quantum-well,  especially  in  electrically  injected  lasers,  and  where  optically  pumped  lasers 
absorb  significant  pump  in  the  barrier  regime,  typically  to  increase  net  absorption.  There  is  a  significant 
inefficiency  in  this  regime  due  (known  as  the  “Richardson  Factor”)  to  electron  momentum  mismatch  upon 
transition  from  the  barrier  to  the  well  region. 

The  devices  numerically  investigated  and  reported  in  Figure  12  is  in  a  distinctly  different  regime  than  that 
previously  reported[25],  in  that  the  final  state  in  the  CHCC  Auger  process  is  confined  within  the  quantum 
well.  In  this  region,  Within  this  region,  coupling  is  very  dependent  on  band-structure.  Of  note  is  the  region 
where  £3  —  E\  is  close  to  Eg,  which  promotes  resonant  transitions;  the  transition  probability  between  these 
excited  electronic  states  increases  dramatically,  causing  the  Auger  lifetime  to  drop,  thus  reducing  rc  by  a 
factor  of  2. 

In  Figure  2b,  we  hold  the  bandgap  approximately  constant  (Alo.35Gao.65AsySbi_y  bariers),  and  adjust  the 
arsenic  concentration  of  the  barriers,  which  has  the  largest  on  strain.  We  can  adjust  the  strain  in  the  barrier 
from  positive  to  negative  values  by  varying  the  amount  of  As.  We  see  that  influence  of  the  band  energies 
within  the  quantum  well  remains  small,  and  this  reflects  little  change  in  the  Auger  and  spontaneous  lifetimes. 
We  conclude  these  structures  are  not  dramatically  strain  driven. 

We  intend  to  grow  several  structures  at  University  of  Missouri  (much  like  the  previous  investigation  at  2  /x m), 
to  confirm  these  predictions  in  the  confined  excited  state  structures.  This  will  also  hopefully  provide  a 
suitable  platform  for  demonstrating  high  power  devices,  at  an  unusual  wavelength  for  type-I  quantum  well 
lasers. 
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3.7.1  Type-I  Method:  Strain  Engineering 


Instead  of  being  bound  to  preset  strains  based  on  the  quantum  well  chosen  and  the  lattice  constant,  relative  to 
a  binary  compound  (e.g.  GaAs,  GaSb,  indium  phosphide  (InP),  etc.),  design-space  would  be  dramatically 
opened  up  if  the  strain  were  considered  another  free  parameter.  While  there  are  problems,  which  will  be 
identified,  it  is  possible  to  accomplish  this  through  use  of  metamorphic  buffer  layers,  which  allow  for  nearly 
full  relaxation  of  a  ternary  crystal  on  a  binary  substrate.  This  concept  is  based  on  Reference  27,  in  which  a 
peak  power  of  several  hundred  mW  has  been  demonstrated,  relying  on  metamorphic  buffer  layers  to  adjust 
the  semiconductor  crystal  lattice  constant  to  allow  for  different  material  composition  and  strain  profiles  to  be 
grown  on  it. 

NLCSTR  performed  a  systematic  optimization  strategy  and  found  regions  which  should  improve  the 
performance  by  more  than  about  a  factor  of  20  when  compared  to  Reference  27  and  should  allow  for 
Watt-level  room  temperature  cw  performance.  NLCSTR  used  its  fully  microscopic  many -body  models  to 
calculate  the  intrinsic  material  properties  of  various  quantum  well  structures  for  target  wavelengths  of  about 
2.9-3. 2  fim.  All  intrinsic  properties  relevant  for  laser  performance  such  as  gain  and  absorption  as  well  as 
carrier  losses  due  to  radiaitive  and  Auger  processes  were  calculated.  The  wells  are  assumed  to  be  grown  on 
metamorphic  buffers  that  lead  to  an  unstrained  Alo.64Ino.36Sb  cladding. 

Calculations  were  performed  for  target  wavelengths  of  2.9  /xm  (that  of  Reference  27),  and  3.2  /xm,  showing 
scaling  to  longer  wavelengths.  The  well  widths  were  varied  between  6  nm  and  13  nm.  Five  different  barrier 
compositions  were  assumed.  Data  were  calculated  for  temperatures  ranging  from  250  K  to  400  K  and 
intrinsic  carrier  densities  from  the  low  absorptive  to  the  high  gain  regime  -  in  total  27  different  structures 
were  analyzed. 

Type  I  quantum  well  structures  for  this  wavelength  regime  are  limited  by  strong  carrier  losses  due  to  Auger 
processes.  Thus,  NLCSTR  focused  on  looking  for  structural  variations  that  lead  to  reduced  Auger  losses. 
Figure  13a  shows  the  Auger  loss  time  as  function  of  the  intrinsic  carrier  density  at  300  K  for  the  wells  that 
were  found  to  yield  record  performance  in  Reference  27  (dashed  line),  as  well  as  for  structures  of  various 
level  of  optimization.  It  was  found  that  by  varying  the  well  width  and  barrier  height,  the  Auger  loss  time  can 
be  reduced  by  almost  three  orders  of  magnitude  for  a  given  carrier  density.  This  dramatic  dependence  is 
rather  surprising  and  will  be  explored  experimentally  to  confirm  these  results.  Even  if  the  reasoning  is  found 
to  be  sound,  total  non-radiative  carrier  lifetimes  will  be  shorter,  driven  by  terms  such  as  SRH  losses. 

Figure  13b  compares  the  threshold  current  density  for  the  optimized  structures  from  Figure  13a  to  that  of  the 
previously  grown  device[27].  The  Auger  losses  are  reduced  with  increasing  optimization.  However  as 
discussed  above,  this  also  increases  the  threshold  carrier  density  and  changes  the  radiative  loss  current.  Thus, 
despite  a  reduced  Auger  loss  in  “Optimization  1”  as  compared  to  the  well  in  the  previous  study  the  threshold 
current  is  not  reduced  in  that  well. 

This  also  leads  to  the  fact  that  the  threshold  current  reduction  is  not  as  strong  as  the  reduction  of  the  Auger 
losses.  However,  although  the  reduction  of  the  threshold  is  about  one  order  of  magnitude  smaller  than  the 
reduction  of  the  Auger  losses,  the  threshold  reduction  is  still  very  significant  -  up  to  about  a  factor  25  for 
“Optimization  4.”  Further  study  seems  to  indicate  the  carrier  loss  is  dominated  by  spontaneous  emission;  this 
situation  is  similar  to  the  one  typically  observed  in  lasers  for  wavelengths  of  1.3  /i  m  and  shorter.  This  leads 
to  high  Tq  values  and  should  allow  for  multi-watt  level  operation  at  room  temperature. 

University  of  New  Mexico,  Center  for  High  Technology  Materials  (UNM)  studied  the  growth  of 
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Figure  13.  Equivalent  thresholds  and  simulated  Auger  lifetimes  for  A  ^2.9  pm  devices 

a)Simulated  Auger  lifetimes  as  a  function  of  sheet  carrier  density  for  structures  grown  lattice- 
matched  to  Al0mIn0.36Sb.  These  include  the  previously  grown  structure  ( black  dash)  operating 
at  2.9  pm,  and  perturbations  ( colored )  which  amount  to  simply  lower  barrier  heights.  The  dots 
mark  the  threshold  carrier  density  for  typical  intrinsic  optical  losses,  b)  ( Equivalent )  threshold 
current  densities  for  optimized  structures. 


metamorphic  buffers  on  GaAs  and  GaSb  substrates.  It  was  found  that  substrates  based  on  GaSb  show  less 
cross-hatching.  This  indicates  a  better  growth  quality  with  less  defects  due  to  strain  mismatch  between  the 
top  buffer  layer  and  the  substrate.  UNM  continued  this  study  investigating  improved  growth  conditions  for 
growing  on  GaAs.  Lower  substrate  temperature  and  higher  Sb-flux  were  found  to  yield  improved  growth 
quality. 

Growth  results  were  also  compared  to  results  obtained  by  using  an  IMF  layer  to  relax  the  strain  between 
substrate  and  buffer  instead  of  a  direct  buffer  growth,  a  technique  introducted  in  Section  3.6.  Characterization 
of  the  IMF  result  as  well  as  the  direct  buffer  growth  on  GaAs  and  GaSb  using  AFM  and  Nomarski  showed 
virtually  no  change  in  the  surface  quality  between  the  three  techniques.  Both  characterizations  showed 
extremely  rough  surfaces  (greater  than  20  nm  root-mean- squared  (RMS)  via  AFM).  The  X-ray  diffraction 
begins  to  differentiate  the  techniques,  with  the  IMF  interface  technique  clearly  showing  the  GaSb  peak  next 
to  the  broad  metamorphic  peak  and  the  straight  metamorphic  buffer  approach  showing  a  single  broad 
metamorphic  peak.  The  major  consideration  when  moving  an  antimonide-based  structure  to  GaAs  is  the 
optical  properties.  Using  cryogenic  PL  to  measure  the  intensity  of  a  pair  of  InGaSb  quantum  wells  grown  on 
the  different  metamorphic  buffers  showed  that  the  technique  that  yields  the  best  PL  is  the  straight  buffer 
growth  on  GaSb  substrate  in  Figure  14. 

In  general,  the  concept  relies  on  the  fundamental  premise  that  lower  barrier  structures,  which  do  little  to 
confine  the  carriers,  also  have  a  exceedingly  weak  involvement  from  Auger  recombination  due  to  the 
difficulty  in  meeting  the  necessary  momentum  change  in  order  to  make  the  appropriate  conduction  band 
transition,  introduced  in  Section  3.2.  In  such  structures  the  carriers  are  not  strongly  confined  in  the  wells  and 
significant  numbers  of  holes  and/or  electrons  can  spill  from  the  well  into  the  barrier,  resulting  in  a  leakage 
current,  thus  reducing  the  quantum  well  efficiency. 
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Figure  14.  Measured  PL  from  77  K  to  300  K  for  devices  operating  at  approximately 
3  /mi. 

a)  300K  PL  at  different  pump  powers,  b)  77K  PL  at  different  pump  powers  for  Ino^Gao^Sb 
QWs  on  Alo  QffnO'ZQSb/GaSb  metamorphic  buffer  layer  on  GaSb. 


Figure  15.  XRD  for  metamorphic  buffer  layer  grown  material  designed  to  operate  at 
approximately  3  pm. 

From  The  structure  is  converted  from  GaSb  to  AlSb  with  layers  slowly  introducing  increasing 
amounts  of  indium  into  the  material  The  lattice-constant  is  also  given  for  each  layer  step. 


22 

Approved  for  public  release;  distribution  unlimited. 


bandgap  wavelength  (pm) 


c n 
c 

CD 


CD 


CD 

CD 

< 


^  103 


a) 


carrier  density  (cm'2) 


b) 


Alxlni-xGao.4Sb 


Figure  16.  Metamorphic  Buffer  Layer  grown  material  designed  at  approximately  3  pm. 

a)  Auger  lifetime  calculated  as  a  function  of  sheet  carrier  density  which  include  a  bulk  treatment 
of  the  barrier  material  ( solid  lines)  as  well  as  a  quantized  subband  model  of  the  barriers  ( dashed- 
lines).  b)  Calculate  threshold  current  density  as  a  function  of  barrier  aluminum  composition. 


To  simulate  these  structures,  a  thorough  model  similar  to  to  that  introduced  in  Section  3.1  is  used  to  calculate 
gain  and  losses  on  this  structure.  Because  the  calculations  are  built  up  from  basic  principals,  there  is  little 
need  to  alter  the  theoretical  approach  to  calculating  the  gain  on  a  non-binary  substrate.  Besides  Auger  losses, 
free  carrier  absorption  (intraband  absorption)  losses  are  generally  expected  to  limit  the  performance  of 
devices  in  this  wavelength  regime.  Additionally,  for  realistic  applications,  the  far  field  broadening  has  to  be 
limited,  as  previous  designs  would  result  in  impractically  high  broadening [27]. 

Beginning  from  designs  presented  in  Reference  [27],  devices  were  simulated  adjusting  the  barrier  material 
composition,  which  adjusts  the  A Ec  and  A Ev.  The  initial  design  is  similar  to  the  “deep-well”  structure 
presented  in  Section  3.1,  although  at  a  longer  wavelength.  Figure  16a  shows  the  increase  in  Auger  lifetime  as 
a  function  of  carrier  density  . 

Typically,  the  use  of  a  bulk  model  is  used  to  simulate  barrier  states,  thus  easing  the  calculation  requirements. 
However,  it  seemed  that  based  on  results  in  Figure  16a,  there  is  a  significant  difference  in  the 
Auger-calculated  when  including  subbands  in  barriers.  Such  a  discrepancy  between  the  results  for  both 
models  had  not  been  encountered  in  testing  the  codes  for  materials  at  shorter  wavelengths  -  around  1  micron. 
In  the  latter  case,  the  results  of  the  subband  model  where  found  to  depend  somewhat  on  the  barrier  width  and 
converged  toward  the  results  of  the  bulk  model  in  the  limit  of  wide  barriers.  The  reason  for  the  discrepancy 
between  the  results  of  the  two  models  for  the  3.x-micron  materials  is  yet  to  be  determined.  Tests  of  the 
subband  model  for  various  barrier  widths  have  shown  an  insignificant  dependence  and  did  not  indicate  a 
convergence  toward  the  bulk  model  results.  However,  when  comparing  to  literature [14],  it  seems  that  the 
barrier  subband  model  more  faithfully  predicts  the  published  results.  One  should  note  this  is  only  as  good  as 
the  reported  growth  structure  matches  what  is  grown. 
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Table  2.  Wall  plug  efficiency  characterization 

Tabulated  laser  performance  for  typical  published  results  at  around  2  pm  and  3  pm. 
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3.7.2  Practical  Implications 

Within  this  discussion,  it  is  easy  to  lose  sight  of  how  this  reduction  in  parasitic  carrier  losses  affects  the 
device.  Figure  8  shows  that  reducing  Auger  currents  directly  affects  the  lasing  threshold  of  the  device.  For 
one-dimensional  devices,  and  excluding  changes  to  thermal  conditions,  once  the  threshold  gain  has  been 
achieved,  the  slope-efficiency  is  only  dependent  on  the  fraction  of  out-coupling  loss  to  total  loss.  As  such, 
wall  plug  efficiency  (WPE)  may  be  presented  as: 

VWPe(J  >  Jth )  =  jy  —  ,  (6) 

where  t]d  is  the  differential  efficiency  of  the  laser,  and  V  is  the  device  voltage.  Typical  values  for  these 
parameters  are  given  in  Table  2. 

When  analyzing  improvements  to  a  device’s  wall  plug  efficiency,  we’d  like  to  compare  the  efficiency  of  two 
devices  at  the  same  output  power.  Therefore,  it  helps  if  we  restate  Equation  6  as  a  function  of  output  power: 

VWPE(Pout  >  0)  =  (7) 

We  then  plot  a  ratio  of  Equation  7,  for  two  different  values  of  Jth  wherein  the  parasitic  recombination  rates 
change  the  threshold  current.  This  improvement  can  be  dramatic  for  smaller  powers,  and  asymptotically 
reaches  unity  for  large  powers.  The  lasers  defined  in  Table  2  are  plotted  in  Figure  17. 
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4.0  TYPE-II  QUANTUM  WELL  SURFACE-EMITTING  DISTRIBUTED  FEEDBACK  LASER 


To  counter  the  Auger  problem  introduced  in  Section  3.0  and  extend  to  longer  wavelengths,  it  becomes 
advantageous  to  look  to  a  non-homogeneous  quantum-well  structure.  The  type-II  “interband”  cascade  laser 
achieves  this  by  confining  the  holes  in  a  center  “well”  while  confining  the  electrons  in  layers  immediately 
adjacent.  If  the  layers  are  sufficiently  thin,  the  wavefunction  overlap  can  be  considerable,  as  much  as  ^70% 
that  of  a  type-I  quantum-well. 

Carriers  may  be  optically[29,  30]  or  electrically [31]  excited.  Optically  pumped  lasers  have  demonstrated 
outstanding  results  and  are  predicted  to  be  able  to  scale  from  2.5-9  pm[32]  at  cryogenic  temperatures.  At  low 
temperatures,  these  devices  exhibit  multi-Watt  performance,  but  this  drops  off  dramatically  at  higher 
temperatures[30].  The  electrically  injected  lasers  are  most  often  employed  in  a  cascade  structure  (quantum 
wells  excited  serially),  and  this  realization  has  been  termed  the  interband  cascade  laser  (ICL).  ICLs  have  been 
employed  at  temperatures  up  to  340°C[31]  and  have  demonstrated  over  0.5  W,  cw  operation. 

In  this  AFRL/RY  type-II  quantum-well  laser  project,  the  goal  was  to  investigate  the  employment  of  type-II 
quantum  wells  with  a  surface-emitting  distributed  feedback  (SE-DFB)  cavity  embodiment.  Two  of  the  most 
important  steps  to  develop  the  SE-DFB  laser  are  to: 

1.  develop  an  appropriate  exposure  technique  and  process  for  the  grating  features 

2.  develop  an  effective  technique  to  transfer  this  pattern  from  photoresist  into  semiconductor 


4.1  Type-II  Method:  Active  Region  Modeling 

In  discussions  with  NLCSTR,  we  are  pleased  to  discover  that  Simulase®,  the  software  we  are  using  for 
rigorous  many -body  gain  and  carrier  loss  calculations  in  type-I  structures  (outlined  in  Section  3.0)  is  also 
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Figure  17.  Improvements  to  wall-plug  efficiency. 

Analytical  wall-plug  efficiency  improvements  based  on  reasonable  values  for  reduction  in 
parasitic  currents  such  as  Auger  and  spontaneous  recombination.  The  black  curve  represents 
improvements  to  a  2- pm  laser,  and  the  red  curve  is  an  approximate  3- pm  laser,  both  who's 
values  are  given  in  Table  2. 
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Figure  18.  Standard  type-II  quantum  well  laser  band  schematic. 

The  three-layer  structure  creating  the  quantum  well  is  composed  of  a  hole  confinement  layer 
flanked  by  electron-confinement  layers.  When  these  are  sufficiently  small ,  the  wavefunctions 
overlap  significantly,  aiding  in  device  efficiency. 


able  to  adequately  predict  gain  in  type-II  quantum  wells.  The  first  step  is  to  compare  published  active  regions 
in  the  literature  to  those  predicted  by  the  software.  This  effort  is  somewhat  challenged  because  the  reports  of 
ICLs  are  rarely  thorough  and  often  do  not  provide  sufficient  detail  of  the  active  region  to  simulate  faithfully. 

The  primary  concern  in  type-II  quantum- well  lasers  is  wavefunction  overlap.  The  quantum-well  is  a 
three-layer  structure,  flanked  by  barriers.  Hole  confinement  (valence-band  vacancies)  is  typically 
accomplished  in  a  small  central  layer,  and  is  most  often  flanked  by  two  narrow  electron  confinement  layers 
(electrons  confined  in  the  conduction  band).  Figure  18  shows  a  typical  structure  designed  for  electrical 
injection  after  Reference  11.  The  design  often  relies  on  strong  confinement  of  the  carriers  within  the  three 
wells  through  use  of  wide-bandgap  barriers  made  of  aluminum  arsenide  (AlAs).  While  this  is  a  very  effective 
approach  for  electrically  excited  devices,  optically  excited  devices  require  narrower  gaps  in  which  the  pump 
light  undergoes  single-pass  absorption.  Unfortunately,  as  one  reduces  the  bandgap  of  the  barriers  to  absorb 
pump  photons,  both  the  electron-  and  hole-confinement  becomes  more  bulk-like  and  the  overlap  (as  well  as 
probability  of  locating  the  carrier  within  the  well)  drops  dramatically. 


4.1.1  Double-Confined  Structure 

Potentially,  to  utilize  narrow-bandgap  barrier  material,  it  becomes  advantageous  to  implement  a  narrow-band 
absorbing  layer  within  a  secondary  confinement  region  (notated  as  “secondary  confinement”  in  Figure  19) 
which  confines  both  the  holes,  and  to  a  lesser  extent  the  electrons  upon  weak  confinement  of  the  narrow-gap 
“barriers”  of  the  absorption  layers.  Figure  20a-d  show  the  wavefunctions  of  the  lowest  excited  wavefunctions 
of  the  conduction  and  valence  bands  as  a  function  of  the  secondary  confinement  layer  material.  The  spacer 
layer  was  the  traditional  AlAs  material.  In  the  case  of  small  indium  compositions  (Figure  20a, b),  the 
secondary  confinement  structure  is  very  influential  in  defining  the  wavefunction,  where  the  center  hole  well 
is  readily  coupled  to  the  secondary  confinement.  For  large  indium  concentrations,  the  smaller  bandgap 
(Figure  20c, d)  promotes  hole  wavefunctions  which  are  mostly  confined  within  the  center  hole  well,  although 
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Figure  19.  Conceptual  approach  to  the  double-confined  structure 


the  electron  wavefunction  remains  strongly  influenced  by  the  secondary  confinement  structure.  Not 
surprisingly,  as  the  wavefunctions  become  more  localized  (e.g.  Figure  20d)  the  gain  increases.  This  is 
depicted  in  Figure  20e,  which  shows  the  calculated  absorption  of  such  structures  given  in  Figure  20a-d. 

It  should  be  noted  these  structures  are  very  influenced  by  the  width  of  the  secondary  confinement  structure. 
In  order  to  achieve  sufficient  pump  absorption,  the  secondary  confinement  layer  must  be  too  wide 
(Beer-Lambert  law  says  absorptionoc  1  —  e~ad  where  a  is  the  material  absorption  and  d  is  the  material 
thickness).  The  designs  in  Figure  20  have  pump  absorptions  on  the  order  of  1%  because  the  total  thickness  is 
too  small,  even  for  10  repeats  of  the  structure.  We  must,  therefore  utilize  a  more  traditional  structure  utilizing 
extremely  thick  absorber  layers  without  the  double-confined  structure. 

4.1.2  Single-Confined  Structure 

The  weakly  confined  structures  created  with  the  narrow-bandgap  barrier  materials  limit  the  short-wavelength 
limit  because  at  some  point,  the  electrons  cease  to  remain  confined.  Figure  21  shows  result  of  parametrically 
adjusting  the  hole  well  thickness  (reported  in  monolayers)  and  its  influence  on  absorption/gain. 

A  few  notes  as  a  result  of  this.  The  thicker  hole-well  decreases  the  AE^,  therefore  this  results  in  a  longer 
wavelength.  The  model  has  a  systematic  error  (a  structure  we  know  should  emit  at  4.3  fim  is  predicted 
numerically  to  emit  at  6.98  /xm).  Although  the  actual  source  of  this  error  is  still  being  investigated,  we 
suspect  the  exact  subband  positions  are  not  being  calculated  accurately.  While  these  provide  a  less  relevant 
effect  for  type-I  lasers,  type-II  lasers  are  particularly  sensitive  to  these  positions.  At  thinner  hole  well 
thicknesses,  the  holes  are  no  longer  confined  as  they  couple  to  the  IncuGao.eAsSb  material. 


4.2  Device  Design 

We  have  focused  on  two  structures,  one  with  five  wells,  and  one  with  seven.  These  should  provide  sufficient 
pump  absorption  and  gain.  The  designs  result  in  decent  modal  gain  overlap  (~4%),  and  single  transverse 
mode  operation.  The  most  difficult  design  parameter  is  the  modal  grating  overlap  which  will  provide  the 
“mirror”  in  the  distributed  feedback  laser.  This  overlap  would  ideally  be  about  1.5  x  10-4,  although  our 
present  designs  are  limited  to  about  a  Tgrat  of  5  x  10-5.  These  structures  have  been  sent  to  be  grown 


27 

Approved  for  public  release;  distribution  unlimited. 


2000 


15%  In  20%  In 

i  1  i  1  i  1  i  1  i  i  i  1  i  1  i  1  i 


30%  In 


40%  In 


-2000 


e) 


2468  10  2468  10  2468  10  2468  10 

wavelength  (pm) 


Figure  20.  Absorption  and  wavefunction  calculations  for  type-II  quantum  well  active  re¬ 
gion 

Calculations  showing  wavefunction  calculations  along  with  absorption  (negative  gain)  cal¬ 
culations  utilizing  many -body  theory. The  curves  in  (e)  represent  a  different  carrier  density: 
O.l(blue),  1,  2,  3,  4,  5,  and  10(red)xl012cm~2. 
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Figure  21.  Influence  of  hole-well  thickness 

Calculations  showing  absorption  and  influence  of  the  hole-well  thickness.  The  curves  represent 
a  different  carrier  density:  O.l(blue),  1,  2,  3 ,  4,  5 ,  and  10(red)xl012cm~2. 

commercially  and  returned  to  fabricate,  mount,  and  test. 

4.3  Type-II  Device  Modeling  Results  and  Discussion 
4.3.1  Introduction 

When  designing  SE-DFBs,  the  most  critical  influence  over  the  device  are  the  grating-waveguide  parameters 
£,  n,  and  asurf  in  the  coupled-wave  equations.  Their  derivation  begins  with  approximations  for  estimating 
these  parameters.  Rough  rules  of  thumb  for  approximating  their  values  are  derived.  This  is  followed  by  an 
outline  of  a  potential  technique  for  computing  the  parameters  for  general  grating-waveguides  (not  just 
rectangular  gratings). 


4.3.2  Coupled  Differential  Equations 


The  differential  equations  suitable  for  coupled  distributed  feedback  (DFB)s  which  include  surface  emission 
are  given  by  the  seminal  paper  of  Noll,  and  Macomber[33].  First,  the  total  field  is  assumed  to  be  of  the  form: 


E{x,  z)  —  R(z) 


E+(x,z)  +  E  (x,z) 


+  S(z ) 


E+(x,  z)  —  E  (x,  z) 


(8) 


where  E+/~  are  considered  the  right-  and  left-going  waves,  and  R  and  S  are  defined  as  normalized 
amplitude  functions.  With  these  considered,  along  with  the  coupled  wave  equations,  the  coupled  partial 
differential  equation  (PDE)  set  for  R  and  S  are  solved  by  the  same  manuscript  to  be: 


*ffl-[a(z)  +  iS  +  i(}R(z) 
oz 

-^^--Hz)  +  iS  +  i(]S(z) 


inS(z) 

(9) 

inR(z) 

(10) 
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where  S  defines  the  finite  length  of  the  resonant  wave  vector,  a  is  the  waveguide  propagation  loss  (mentioned 
below),  and  (  is  defined  relative  to  the  infinite  length  problem. 


4.3.3  Approximate  Estimates  of  Grating  Parameters 

A  general  approach  to  finding  the  grating-waveguide  parameters  in  a  manner  consistent  with  rigorous 
boundary  conditions  was  described  by  Noll  for  transverse  electric  (TE)  modes  of  rectangular  gratings[33].  A 
later  paper  by  Finger  et  al.  adopted  the  Noll  approach  for  transverse  magnetic  (TM)  modes  of  rectangular 
grating  waveguides [34]  -  a  formalism  that  would  be  necessary  if  deploying  the  SE-DFB  for  QCLs.  Both 
papers  solve  for  the  two  infinite-length  solutions  (i.e.have  exact  periodicity  along  the  z-axis).  The 
infinite-length  gain-frequency  eigenvalues  are  then  used  to  find  the  grating  parameters  for  the  coupled-wave 
equations  that  describe  finite-length  devices. 

Although  no  code  is  available  for  calculating  the  grating  parameters,  we  are  able  to  extrapolate  from 
calculations  that  were  tabulated  in  Noll.  What  we  do  have  available  are  planar-waveguide  codes.  A  planar 
waveguide  code  represents  the  grating  layer  with  a  duty-cycle  weighted  average  of  the  dielectric  constant 
values  (i.e.n2)  on  each  side  of  the  grating  layer.  What  planar-waveguide  codes  can  find  is  the  optical 
confinement  in  the  metal  as  well  as  a  value  of  the  metal  loss. 

An  approximate  estimate  of  the  surface  loss  coefficient  asurf  is  given  by 

ow/~8r#/A  (ii) 

where  Tg  is  the  optical  confinement  fraction  of  the  waveguide  mode  in  the  grating  layer  and  A  is  the  grating 
periodicity.  This  formula  is  a  rough  fit  (within  a  factor  of  two)  to  past  modeling  analysis  results  of  830-nm 
devices.  Extrapolation  to  longer  wavelength  is  included  in  the  grating  periodicity  denominator  ineq:alpahsurf 
as  per  the  wavelength-coordinate  scale  invariance  of  the  wave  equation.  The  effect  of  varied  metal  index  is 
less  clear. 

To  gain  insight  into  the  coupled- wave  equation  parameters  k  and  (  it  is  useful  to  review  their  derivation. 
They  are  found  from  the  mode  gain-frequency  eigenvalues  a  +  i5  of  the  two  solutions  E+  and  E~  from 
idealized  infinite-length  problem  (equivalent  to  exact  periodic  conditions).  These  values  are  related  to  the 
grating  parameters  as  follows 

k  —  ^  —  ioT  +  S~  —  5+)  (12) 

and 

C  —  2  ^ a~  ~  ~  ^+)  (13) 

The  imaginary  part  of  the  eigenvalues  can  be  partitioned  into  metal  loss  and  surface  loss  parts 

=  arnetal  +  asurf  (14) 

whereas  the  E~  solution  has  no  surface  loss 

=  ametal  (15) 

The  two  metal  loss  parameters  cannot  be  calculated  without  a  complete  code  but  their  magnitude  in  Noll  is 
typically  on  the  order  of  the  mode  loss  from  planar  models. 

& planar  2  metal  ^ metal)  (16) 
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The  a+  metal  loss  was  always  found  to  be  greater  than  a  ,  although  the  two  values  tend  to  be  closer  together 
for  weak  gratings,  which  will  be  the  case  for  the  type-II  SE-DFB  designs  created  for  this  effort. 

The  real  part  of  k,  and  ('  are  derived  from  the  S  part  of  the  eigenvalues  and  cannot  be  estimated  from  a  planar 
model.  However,  it  turns  out  that  when  the  grating  duty  cycle  is  around  50  percent,  which  coincidentally 
gives  the  strongest  surface  output,  5~  ~  5+  and  Equations  12  and  13  can  be  adjusted  accordingly. 

Putting  all  these  approximations  together  we  arrive  at  the  rule  of  thumb  for  k: 

hi  ~  —  Qtsurf-  (17) 

The  phase  of  k,  is  dependent  on  the  coordinate  system.  It  rotates  over  2n  radians  when  the  z-coordinate 
origin  shifts  over  Az  =  A/2.  The  exact  coordinate  origin  is  not  important  in  most  analysis  problems, 
therefore  neither  is  the  absolute  phase  of  hi.  The  insignificance  of  the  phase  of  k  is  verified  in  device  models. 
It  only  becomes  important  in  cases  where  devices  with  facet  reflectors  are  modeled. 

The  value  of  Im(()  is  estimated  from  loss  parameters  whereas  Re(()  is  unknown  without  the  5± 
eigenvalues.  However,  Re(()  can  be  ignored  because  it  simply  redefines  the  zero  of  the  S  axis.  Hence,  we 
arrive  at  the  approximation  for  (: 

C  ~  2°^surf  ^ planar •  (18) 

Equations  11,  eq:kappa  and  eq:zeta  provide  a  basis  to  proceed  with  device  calculations.  However,  there  is 
still  uncertainty  in  the  asurf  estimator  in  Equation  1 1  because  it  is  based  on  a  different  metal-index  value 
than  occurs  at  4  /im  wavelength. 

4.3.4  General  Computation  Method 

The  single-layer  rectangular  gratings  models  in  Noll  and  Finger  led  to  difficult  analytical  and  numerical 
problems.  The  root  of  the  difficulty  arises  from  the  classical  theory  of  deep  gratings  where  the  Rayleigh 
method  breaks  down.  The  field  in  the  so-called  “selvedge”  region  must  have  its  own  form  and  must  be 
boundary  matched  to  the  adjacent  layers.  For  example,  a  separate  grating-layer  eigenfunction  method  was 
employed  in  Noll.  But  just  finding  the  eigenvalues  within  the  grating-layer  was  difficult.  The  additional 
problems  of  matching  boundary  values  at  the  edges  of  the  grating  layer  and  finding  the  laser-resonance 
eigenvalues  made  automated  computation  problematical. 

Instead,  we  would  like  to  have  a  method  that  efficiently  and  reliably  models  any  grating-waveguide  structure  - 
not  just  single-layer  rectangular  gratings.  It  would  then  be  possible  to  analyze  structures  with  mixed  layers  of 
dielectric  and  metal  that  are  achievable  with  modern  micro-fab  technology.  With  such  an  analytical  tool  it 
may  be  possible  to  design  practical  grating-waveguides  that  achieve  both  strong  coupling  and  low  metal  loss. 

One  possibility  examined  in  the  past  is  to  put  the  problem  on  a  grid.  Steady-state  solutions  can  be  found  by 
iteration  or  simulation  of  time  evolution.  While  such  techniques  can  work,  they  have  a  problem  in  that  the 
frequency  eigenvalue  should  be  accurate  to  5  ^0.01  cm-1  out  of  nk  ^1  x  104  cm-1  or  one  part  in  106.  But 
errors  from  finite  grid  size  and  iteration  noise  can  be  on  this  order  of  magnitude.  Moreover,  grid  techniques 
are  computationally  intensive  (i.e.very  slow  to  run). 
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Floquet-expansion  methods  such  as  used  in  Noll  and  Finger  can  find  the  resonance  eigenvalues  to  high 
precision.  The  issue  is  instead  how  to  accurately  and  reliably  solve  the  equations  for  each  order  function  with 
arbitrary  grating  shapes.  The  order  functions  are  defined  by  the  set  of  coupled  differential  equations 


d?  °° 

-j^Em(x)  -  K2MEm(x)  +  k2n2(z)Em(x)  =  k2  ^  eq(z)Em-q (x)  (19) 

q=— 00, /0 

The  grating  structure  is  contained  in  the  dielectric  function  e(x,  z)  the  Fourier  transform  of  which  eq(x) 
which  appears  in  the  previous  equation.  More  formally,  the  problem  is  to  solve  the  coupled  set  of  differential 
Equation  19  subject  to  the  boundary  conditions  that  there  are  outgoing  waves  only  -  whether  radiating  or 
evanescent.  All  of  the  non-grating  layers  have  the  general  form  of  solution 


Em{x)  =  Aeipx  +  Be~ipx 


(20) 


so  that  the  outgoing-wave  boundary  conditions  can  be  expressed  as 

C  A  /D^PtTI  ry  ry 

rp  /  \  _  I  x  ^  max 

j  TD  -ipmX  ry  ^  ry  . 

^  \  ^ min 

where  xmin  and  xmax  are  the  interfaces  to  the  outer-most  layers  and  where 

Pm  =  +y/nW  -  Kg. 


(21) 


(22) 


The  eigenvalue  problem  of  second-order,  ordinary  differential  equation  (ODE)  with  boundary  conditions  on 
both  sides  of  the  field  generally  have  well-defined  solutions. 


4.3.5  Shooting-Method  Numerical  Technique 

In  the  fully  self-consistent  solution,  each  equation  for  Em(x)  in  Equation  19  constitutes  a  non-homogenous 
ODE.  The  zero-order  equation  (m=0,-2)  is  just  the  standard  planar-  waveguide  equation  when  the  other 
Em(z)=  0.  Only  specific  values  of  (3  =  nk  allow  the  zero-order  equation  to  satisfy  the  outward-going 
boundary  conditions  on  both  sides.  These  are  the  usual  waveguide  modes. 

Planar  modes  are  found  by  testing  trial  /3-values  (also  be  expressed  as  5  =  5/3  value)  to  see  if  they  satisfy  the 
boundary  conditions.  The  classical  method  is  to  evaluate  the  determinant  of  the  boundary-matching  equation 
matrix.  This  works  for  simple  layers  of  uniform  index.  Alternately,  a  “shooting”  method  can  integrate  the 
differential  equation  to  the  opposite  side  of  the  structure  given  two  starting  points  at  one  edge.  These  two 
points  can  easily  be  generated  because  the  form  of  solution  is  analytic,  the  initial  amplitude  is  arbitrary,  and 
the  p-value  is  known  for  any  given  trial-/3  value.  In  general,  the  solution  of  the  far  side  has  the  form  given  in 
Equation  21  and  -  depending  on  which  side  you  start  -  the  value  of  the  A  or  B  coefficient  should  vanish 
when  a  mode  is  found.  Thus  one  can  construct  a  test  function  such  as 


\B\ 2 


(23) 


which  blows  up  when  A  —^0.  The  steepness  of  the  test  function  makes  it  easy  to  zero  in  on  precise 
eigenvalues.  The  small  parameter  7  avoids  division-by-zero  errors.  The  normalization  of  Eq{x )  is  arbitrary 
but  it  is  conveniently  chosen  to  be  power  integrated  to  unity. 
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The  shooting  method  for  finding  Eq(z)  can  still  be  applied  when  the  Em(z)  functions  are  non-zero.  The 
gain- frequency  values  will  just  shift  a  little.  In  fact,  the  gain-frequency  eigenvalues  of  the  zero-order  equation 
when  self-consistent  order  functions  are  found  are  the  very  values  of  arf  that  we  seek  for  computing  k  and  (. 

The  above  discussion  outlines  how  the  zero-order  problem  is  handled.  That  still  leaves  the  other  orders 
m/0,-2.  But  the  problem  is  still  the  same-finding  solutions  to  the  non-homogeneous  ODE  subject  to 
out-going- wave  boundary  conditions.  The  difference  is  that  the  amplitude  of  the  fields  cannot  be  subjected  to 
normalization  but  must  be  consistent  with  the  non-homogenous  term  in  the  ODE.  A  shooting  method  can  be 
used  again  except  in  this  case  there  is  no  eigenvalue  search.  Instead  a  search  can  to  be  done  on  the  real  and 
imaginary  part  of  the  A  or  B  coefficient  launched  on  one  side  of  the  structure  such  that  the  outgoing  wave 
boundary  condition  is  satisfied  on  the  opposite  side  of  the  structure.  This  is  done  for  each  order  and  the 
whole  process  is  iterated  until  it  converges.  We  can  anticipate  that  an  iterative  method  will  converge  quickly 
because  the  grating  layer  is  only  a  perturbation  on  the  main  mode. 

The  above  outlines  the  method.  There  are  many  other  details  such  as  use  of  analytic  solutions  whenever 
possible  to  speed  up  the  integration  and  the  method  of  numerically  integrating  across  the 
metal- semiconductor  interface.  The  TM  problem  gets  a  little  more  complicated  but  can  still  be  expressed  as 
an  ODE  (with  a  new  derivative  term)  and  also  solved  using  a  shooting  method.  While  we  choose  to  apply  this 
type  of  cavity  to  the  type-II  quantum  well,  where  TE  calculations  are  the  most  appropriate,  quantum  cascade 
lasers  are  entirely  TM  polarized,  and  therefore  the  coupled  mode  solution  would  be  required  for  creating  a 
QCL  SE-DFB. 


4.3.6  Optical  Gain  Parameters 


The  material  gain  parameters  including  the  internal-quantum  efficiency  (r/2),  the  current-  gain  coefficient  Aj, 
and  the  gain  offset  b  are  used  in  the  unsaturated  gain  formula: 


g  =  r 


f  Vi  Aj  J 

\^r 


(24) 


where  J  is  the  equivalent  current  density,  and  d  is  the  quantum-well  thickness.  These  and  other  parameters 
can  be  inferred  from  a  publication  on  optical  pumping  of  type-II  semiconductor  lasers[35],  which  described 
an  epitaxial  structure  similar  to  that  in  our  project.  These  devices  had  optical-to-optical  efficiency  of  18%  at 
83  K,  which  corresponds  to  differential  quantum  efficiency  (DQE)  of  0.37  after  correcting  for  the  quantum 
defect  from  1.85  (im  to  3.8  /im.  Other  parameters  given  in  the  paper  include:  the  front-  and  back-facet 
reflectivity  =0.02  and  =0.3,  respectively;  the  mode  loss  a  =3  cm-1  and  the  optical  pumping  “stripe”  area 
W  x  L  ^250x2500  /im.  The  authors  do  not  estimate  the  internal  quantum  efficiency  rji  but  one  can  back 
out  a  value  from  the  DQE  formula  for  front-facet  emission 


DQE  =  rn 


iln  1 


2  L 


U/ 


2L 


RfRb 


T  CM 


(25) 


Equation  25  and  the  data  in  Reference  35  imply  7^=0. 63.  The  threshold  mode-gain  coefficient  can  be 
calculated  from  the  standard  threshold  gain  condition  (gain  equaling  losses): 


»  =  2 Lln 


-2_) 
RfRb  J 


CM 


(26) 
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Table  3.  Calculated  and  assumed  waveguide  parameters  for  two  structures  sent  to  be 
grown. 


KM8_060_WG3_7x 

KM8_060_WG3_5x 

Nw 

7 

5 

neff 

3.45205 

3.28149 

C^metal 

3.0333  cm"1 

6.1513  cm-1 

r 

.005007 

.03233 

r 

L  grat 

4.81  xl0“5 

1.01  xl0“4 

which  predicts  #^=13.2  cm  1.  The  threshold  pump  power  reported  is  Pth  ~2  W,  and  may  be  converted  to 
an  equivalent  threshold  current  density  using  the  following  equation: 

Jequiv  =  ~^X  =  480 A/ cm2  (27) 

where  Xp  is  1.85  /im  used  in  Reference  35.  The  optical  confinement  and  active  thickness  were  not  given  in 
the  paper.  However,  we  have  analysis  of  similar  epitaxial  structures  in  the  early  part  of  this  project  with  seven 
active  zones[25].  That  design  had  and  T  =  9  x  10-3  and  d- 46.2  nm,  referenced  to  66  nm  combined 
thickness  of  the  electron  and  hole  confining  layers  in  each  quantum  region.  The  paper  had  ten  quantum 
regions  but  the  ratio  should  be  about  the  same,  and  therefore  we  will  assume  T= 0.0129  and  d- 66  nm  for  the 
reported  structure.  The  b  parameter  is  unknown  as  we  have  only  one  data  point  to  fit  a  straight  line,  but  we  do 
know  that  b  >0.  We  therefore  guess  that  b  «  300  cm-1  (comparable  to  practical  values  for  GaAs).  Then 
from  Equation  24  we  estimate  Aj  as: 

Aj  =  -T-  (gth  +  Tb) .  (28) 

ViFJth 


4.3.7  Grating  Parameters 

There  are  two  epitaxial  designs  to  be  considered  in  the  device  modeling  runs.  They  will  be  assumed  to  have 
the  same  Aj  and  6,  gain  and  free-carrier  mode-loss  parameters  as  estimated  above.  The  waveguide 
parameters  for  the  two  epitaxial  designs  are  known.  Other  parameters  will  be  assumed  to  have  properties  that 
were  estimated  in  Reference  36.  The  prediction  in  the  device  modeling  runs  that  follows  are  an  attempt  to  be 
more  predictive  of  actual  performance.  Particular  waveguide  parameters  are  given  in  Table  3. 


4.3.8  Device  Designs 

With  estimates  of  the  basic  material  parameters  established,  we  can  now  turn  to  the  simulation  of  specific 
device  designs.  Three  stripe  sizes  were  considered:  300x2000,  200x3000  and  150x4000  fim2  (the  stripe 
area  is  designed  to  be  the  same  in  each  case).  For  each  stripe  design,  both  the  7-  and  5-quantum  region 
materials  were  modeled.  Other  design  details  include  the  unpumped  grating  dimensions,  the  grating 
curvature,  and  the  grating  chirp.  Particular  parameters  of  grating  curvature  will  depend  on  real  waveguide 
designs  as  well  as  thermal  profiles. 
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Table  4.  Numerical  results  of  laser  modeling 


stripe  size 

(P  m2) 

epi 

epi 

M 

9th 

(cm-1) 

DQE 

Jth 

(A/cm2) 

0-0  eff 

@lat 

(deg) 

@long 

(mrad) 

Strehl 

300x2000 

5 

1.77 

19.4 

0.20 

620 

0.088 

15.6 

1.56 

0.11 

300x2000 

7 

1.77 

12.8 

0.21 

666 

0.088 

15.8 

2.50 

0.05 

200x3000 

5 

1.91 

18.3 

0.21 

556 

0.093 

8.0 

0.61 

0.61 

200x3000 

7 

1.91 

11.3 

0.22 

591 

0.098 

7.7 

0.55 

0.51 

150x4000 

5 

2.04 

17.1 

0.20 

521 

0.088 

5.2 

0.48 

0.85 

150x4000 

7 

2.04 

10.4 

0.22 

519 

0.093 

5.3 

0.43 

0.81 

The  total  grating  length  was  chosen  to  be  long  enough  to  capture  most  of  the  emitted  radiation.  However,  the 
total  grating  length  and  magnification  are  constrained  in  that  the  in-plane  virtual  sources  should  be  located 
outside  the  grating.  This  can  occur  when  the  unpumped  part  of  the  grating  is  very  long  or  the  magnification 
is  very  high. 

The  grating  curvature  is  adjusted  to  attain  a  nominal  stripe  magnification  of  M  «  2.  This  rough  rule  of 
thumb  gives  strong  enough  curvature  to  suppress  filamentation  (the  amplitude-phase  coupling  prevalent  in 
many  semiconductor  lasers),  yet  is  small  enough  to  make  the  beam  quality  easy  to  achieve  (ultra-high 
curvature  designs  usually  have  poor  beam  quality). The  grating  chirp  was  used  to  optimize  the  beam  quality 
in  each  case. 

The  grating  was  made  wide  enough  to  avoid  reflections  from  the  edges.  It  was  found  that  an  insufficiently 
wide  grating  led  to  filamentation.  A  similar  effect  was  found  in  a  study  of  electrically-pumped  SE-DFB 
devices.  The  softly  extended  edge  of  the  Gaussian  pump  beam  as  compared  to  a  sharply-defined  electrical 
contact  stripe  makes  this  issue  even  more  important  for  optically-pumped  devices. 

It  is  important  to  note  that  facet  reflection  can  destabilize  devices  both  spatially  and  spectrally.  Past 
GaAs-based  devices  used  absorbers  zones  to  suppress  facet  reflections.  An  alternative  method  may  be  where 
the  laser  longitudinal  mode  is  intentionally  misaligned  with  respect  to  the  crystal  cleave  direction  (tilted)  so 
that  end-reflections  are  deflected  off-axis. 


4.3.9  Modeling  Results  and  Discussion 

Table  4  summarizes  the  results  from  the  numerical  modeling  of  the  six  designs.  The  gratings  in  all  designs 
had  to  be  extended  2.5  mm  beyond  the  end  of  the  stripe  in  order  to  reduce  end  losses.  In  comparison,  they 
only  need  to  extend  0.4  mm  for  typical  l-/xmGaAs-based  SE-DFB  devices,  essentially  scaling  with 
wavelength.  The  longer  end-length  results  from  the  smaller  surface-diffraction  coefficient. 

The  300x2000  fim  “stripe”  designs  lased,  but  had  much  lower  beam  quality  and  more  instability  than  the 
other  stripe  designs.  The  grating  in  that  case  had  to  be  very  wide  to  avoid  outer-edge  reflections,  otherwise 
filamentation  would  set  in.  The  narrower  designs  were  better  behaved  and  had  good  beam  quality,  especially 
the  4  mm  long  designs,  which  achieved  a  Strehl  of  0.8,  performing  very  nicely  for  lasers  of  this  type. 

The  field  and  beam  characteristics  from  one  of  the  4  mm  stripe  designs  are  shown  in  Figures  22  and  23. 
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stripe  &  grating 


Gain  field 


R{z) 


S{z) 


Near  field  at  grating 


virtual  waist  plane 


Figure  22.  Grating  and  field  images  from  the  design  with  150x4000  ^m2  stripe  size  and 
7-quantum  regions 

The  overall  field  dimensions  are  770  /am  by  9000  pm.  The  z-axis  points  to  the  right.  The  nearly 
colorless  near  field  indicates  good  beam  quality  ( phase-represented  by  color).  R(z )  and  S(z ) 
are  the  right-  and  left-propagating  field  slowly  varying  amplitude  functions,  respectively, 
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Owing  to  the  longer  wavelength,  the  virtual  waist  size  is  roughly  four  times  larger  than  with  l-fim  devices. 
The  narrow  5 -degree  lateral  far- field  full  width  half  maximum  (FWHM)  for  the  4-mm  length  design  is 
associated  with  its  high  aspect-ratio. 

Thermal  influences  to  the  refractive  index  were  not  included.  Their  inclusion  would  slightly  lower  the  beam 
quality  but  can  be  compensated  by  fine  tuning  higher-order  terms  in  the  grating  phase  expansion.  Such 
refinements  would  be  necessary  after  experimental  devices  and  estimates  of  heat  load  based  on  experiment 
are  available. 

A  stability  test  was  performed  on  the  150x4000  /j,m2  design  with  the  7-quantum  region  material  design.  The 
test  consisted  of  incrementally  increasing  the  antiguiding  until  the  fields  filamented.  The  antiguiding  value 
reached  o^=4.5  at  filamentation,  which  is  considerably  higher  than  typical  values  measured  in  GaAs-based 
devices  {olb  ~2).  This  numerical  test  result  does  not  guarantee  that  the  actual  device  will  not  filament  but  it 
does  at  least  indicate  that  non-filamentation  is  possible.  The  lower  sensitivity  is  likely  due  to  the  larger 
diffraction-limited  spot  (again  due  to  wavelength  scaling).  The  predicted  ~9%  optical-to-optical  efficiency  is 
somewhat  lower  than  that  obtained  in  Reference  35.  Two  item  underlie  the  low  efficiency:  non-radiative 
recombination  (as  parameterized  by  rji)  and  losses  from  absorption,  both  in  the  metal  and  from  free-carriers. 
The  non-radiative  recombination  is  a  materials  issue  and  will  not  be  addressed  here,  although  we  predict 
there  is  some  error  in  our  gain  prediction  that  will  hopefully  be  refuted  by  experiment. 

To  see  how  different  parameters  affect  the  efficiency,  a  few  trial  cases  were  run.  First  the  5  quantum  well 
150x4000  /im2  (stronger  grating)  was  run  using  the  same  5  cm-1  absorption  coefficient  as  in  the  7 
quantum-well  ~  150x4000  /im2  modeling.  This  produced  a  DQE  of  0.28  versus  0.20  with  10  cm-1  loss. 
This  shows  that  grating  strength  is  in  part  limiting  the  efficiency.  To  see  an  extreme  case,  grating  and 
absorption  parameters  for  a  typical  980  nm  SE-DFB  device  were  run  ( aw/  ^50  cm-1).  This  run  produced 
a  DQE  of  0.45  -  a  very  dramatic  improvement.  Yet  the  combined  metal  and  free-carrier  loss  was  6  cm-1, 
actually  larger  than  for  the  7  quantum  well  150x4000  /xm2  case.  Thus,  it  is  clear  that  grating  strength  is  in 
part  limiting  the  efficiency,  not  just  absorption.  Conceptually,  a  stronger  grating  diffracts  light  out  before  it 
has  propagated  very  far  in  the  medium  thereby  avoiding  loss.  Ideally,  we  should  have 
asurf  &metai  +  a/c,  where  Oifc  is  the  free  carrier  loss. 

But  these  observations  are  only  parametric.  The  metal  loss  is  expected  to  increase  dramatically  as  the  grating 
is  made  stronger,  thus  offsetting  the  gain  from  increased  grating  strength.  Ideally,  we  would  have  grating  that 
is  both  stronger  ( aw/  ^50  cm-1),  and  has  a  very  low  loss  (<10  cm-1).  But  given  the  optical  properties  of 
Au  and  Ag  the  simple  etch-relief,  metal-coat  grating  approach  may  never  achieve  this  goal.  This  points  out 
the  need  for  a  grating  layer  structure  that  could  achieve  sufficiently  strong  coupling.  One  thought  is  air-gap 
grating,  which  has  a  large  index  contrast 

4.4  Fabrication  and  Packaging 
4.4.1  Fabrication  Steps 

Figure  24  presents  the  proposed  steps  to  fabricate  the  SE-DFB  outlined  and  described  herein: 

initial  epitaxy  The  epitaxy  is  grown  by  MBE  on  a  GaSb  substrate.  This  epitaxial  structure  is  on  the  order  of 
five  microns  thick,  and  includes  an  optical  waveguide  in  the  growth  direction,  acting  as  a  slab 
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Figure  23.  Simulated  slit  scans  of  the  far  field,  near  field  and  virtual  waist 

Devices  are  the  150x400  fmr  design  with  and  7 -quantum  regions 
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Figure  24.  Proposed  type-II  SE-DFB  processing  steps. 

a,c,e)  show  3-D  renderings  of  structures,  while  b,d,f)  show  cross-sectional  views,  g)  shows  the 
notional  embodiment  of  the  final  device  under  test.  NOT  TO  SCALE. 
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waveguide  (no  guiding  in  the  lateral  direction).  Within  the  optical  waveguide  there  are  several 
(nominally  five)  type-II  quantum  wells.  The  optical  core  absorbs  the  pump  laser  frequency,  and  these 
excited  carriers  relax  into  the  quantum  wells,  recombining  and  giving  off  stimulated  emission. 

patterning  the  grating  Figure  24a, b  shows  the  epitaxial  structure  on  the  GaSb  substrate  (gray)  as  well  as  a 
commensurate  diamond  structure  (blue),  which  will  facilitate  heat  spreading  and  removal  under  the 
laser.  The  first  step  is  to  create  a  surface-grating  structure  into  the  epitaxial  structure.  When  the  period 
is  properly  selected,  this  will  provide  feedback  (from  a  “forward”  to  a  “backward”  traveling  wave),  as 
well  as  some  portion  emitting  from  the  surface,  and  into  the  substrate.  The  “surface-emitting  mode” 
will  ultimately  be  evanescent,  and  all  of  the  power  will  be  emitted  through  the  substrate.  The  mode 
interaction  with  the  grating  teeth  (previously  described  in  Reference  36)  was  previously  calculated  to 
necessitate  grating  teeth  on  the  order  of  0.4  /im  deep  with  a  period  slightly  larger  than  1  fi m  and  a 
nominal  duty  cycle  of  50%.  While  this  feature  size  is  large  for  electron-beam  lithography,  and  small 
for  photolithography,  we  have  determined  that  at  least  the  first  experiments  will  rely  on  electron-beam 
lithography,  owing  to  the  flexibility  of  the  technique,  allowing  for  multiple  periods  and  active  control 
over  the  duty  cycle.  We  chose  to  transfer  the  patterned  e-beam  resist  into  the  semiconductor  through 
dry-etching,  relying  on  BCl2:Cl2  etch  chemistry.  This  process  yields  smooth  surfaces  and  vertical 
side-walls,  and  should  be  compatible  with  capably  achieved  in-house.  While  we  have  some  previous 
unimpressive  results [37],  we  think  the  ending  results  (with  low  CI2  fraction)  indicate  a  direction 
forward  for  patterning  this  grating  structure. 

metallizing  the  structure  To  assure  the  surface-emitting  mode  is  suppressed  (thus  not  wasting  power),  we 
may  metallize  the  surface  of  the  grating.  Because  this  surface  is  interacting  with  the  intra-cavity  mode, 
and  is  distributed  throughout  the  cavity,  this  metal  surface  has  to  be  of  high  quality  and  purity. 
Moreover,  because  the  surface  is  patterned,  the  metal  must  coat  conforming  to  the  surface,  making  it 
within  the  teeth  and  coating  the  side-walls.  We  numerically  investigated  multiple  metals,  but  either 
gold  or  silver  were  deemed  to  be  similar  to  each  other,  and  superior  to  other  metals [36].  Figure  24c, d 
shows  both  the  epitaxial  structure  and  the  diamond  heat-spreader  coated  with  metal.  To  deposit  the 
metal,  we  either  rely  on  sputtering,  or  evaporation.  The  benefit  of  the  sputtering  is  that  it  is  known  as  a 
naturally  conformal  deposition  technique,  although  the  resulting  metal  quality  is  slightly  lower  than 
evaporation.  Evaporation,  on  the  other  hand,  is  directional,  and  will  not  deposit  conformally  without 
some  planetary  motion  of  the  substrate,  varying  the  line-of-site  from  the  deposition  source  (crucible). 
Whether  the  metal  is  gold  or  silver,  a  thin  metallic  blocking-layer  should  be  put  above  the  metal  in 
contact  with  the  semiconductor.  The  purpose  of  this  layer  (likely  nickel)  is  to  avoid  an  intermetallic 
forming  from  the  solder  layer  (e.g.indium  solder  forms  InAu2  which  is  an  optically  lossy  metal). 
Above  this  layer,  we  likely  want  a  gold  layer  which  would  prohibit  oxidation  of  the  Ni  layer,  and 
provide  a  metal  which  indium  will  wet  to. 

The  other  step  in  this  process  is  solder  deposition.  We  have  the  option  to  use  indium  preforms, 
however  we  have  been  developing  thin  solder  layers  for  the  VECSEL  project  (see  Section  6.2). 
Because  this  solder  is  within  the  thermal  path  for  our  lasers,  and  its  thermal  conductivity  is  low, 
thinner  solder  will  improve  device  operation.  A  reasonable  indium  thickness  for  high  power  devices  is 
around  5  /mi,  While  it  is  easier  to  deposit  thinner  layers,  they  don’t  bond  as  well  because  of 
non-planarity  of  the  diamond,  and  to  a  lesser  extent,  the  semiconductor. 

bonding  the  laser  The  laser  is  bonded  near  or  slightly  above  the  solder  reflow  temperature  (in  the  case  of 
indium,  we  often  use  180°C  ,  though  the  reflow  temperature  is  157°C  ).  We  most  often  perform  this 
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on  a  hot  plate  in  our  test  lab  using  our  bonding-control-unit .  If  more  precise  control  of  temperature  or 
pressure  is  necessary,  or  relative  positioning  is  essential,  we  can  use  out  flip-chip  bonder  within  the 
cleanroom. 

final  manifestation  Figure  24g  shows  a  cross-section  of  the  finalized  structure.  The  pump  (A  ym) 
comes  in  at  a  slight  angle  relative  to  the  substrate,  which  is  transparent  to  the  pump  and  the  lasing 
wavelength.  The  pump  light  is  absorbed  in  the  waveguide  core  and  excited  carriers  result  in 
spontaneous  and  stimulated  emission  from  the  type-II  wells.  The  stimulated  emitted  photons  are 
guided  and  amplified  in  the  optical  waveguide,  and  fed  back  within  this  waveguide  using  the  second 
order  of  the  grating.  A  fraction  of  this  light  is  diffracted  normal  to  the  substrate  and  is  the  laser  output. 
The  heat  is  removed  from  the  “surface”  of  the  epitaxial  structure,  which  at  this  point  is  solder-bonded 
to  a  large  thermal-conductivity  diamond  and  subsequently  removed  through  conduction  cooling,  or 
some  other  means. 


4.4.2  Grating  Patterning 


As  outlined  in  previous  updates,  the  grating  pitch  for  a  second-order  grating  necessary  to  provide  reflection 
and  surface-emission  is  A  ~1.1  /mi,  and  the  duty  cycle  should  be  around  50  %.  Features  of  this  size  lay 
between  resolution  typically  achieved  by  the  easy,  and  manufacturable  photolithography,  and  the  more 
difficult,  although  flexible  electron-beam  lithography. 


While  we  discount  this  approach  for  development,  selecting  a  more  flexible  approach  such  as  electron-beam 
lithography,  one  technique  that  may  be  used  for  more  manufacturable  grating  patterning  may  be  proximity 
photolithography.  Photolithography  is  the  technique  of  choice  for  microfabrication,  primarily  due  to  ease, 
variety  of  photoresist,  and  ease  of  manufacturing,  because  wafers  may  be  exposed  en  masse,  instead  of 
individual  patterns  written.  While  there  are  many  different  types  of  photolithography  (contact,  proximity, 
projection,  deep-UV),  to  some  extent,  they  are  all  limited  to  the  resolution  allowed  by  the  incoherent,  UV 
light  utilized  in  the  system.  Mercury  lamps  are  most  typical  and  often  have  a  filter  about  i-line  (A  =365  nm) 
of  mercury.  This,  combined  by  system  optics,  intimacy  of  mechanical  contact  of  mask  with  the  wafer,  limits 
“manufacturing  features”  to  about  0.9  y m  and  “research  features”  to  about  0.5  ym.  For  us,  then,  the  features 
achievable  by  photolithography  are  very  close  to  the  necessary  resolution  to  create  our  patterns,  but  because 
we  have  the  advantage  of  a  quasi-one-dimensional  periodicity,  we  may  rely  on  phenomena  such  as  Talbot 
imaging  to  aid  us  in  transferring  photomask  images  into  the  photoresist.  In  truth,  because  of  photoresist 
edge-bead  effects,  it  is  difficult  to  get  true  contact  between  the  mask  and  the  photoresist-covered  wafer,  so  we 
have  to  be  prepared  to  have  some  amount  of  air-gap  between  the  two.  Practically,  if  we  can  control  the  air 
gap  distance  (our  Karl  Suss  MA6  has  positioning  accuracy  of  about  1  ym),  proximity  exposure  is  preferable 
because  it  minimizes  any  contact  of  the  mask  with  surfaces  that  may  result  in  mask  feature  degradation. 


To  analyze  the  proximity  exposure,  we  use  the  scalar  Rayleigh-Sommerfeld  diffraction  integral[38].  This 
approach  allows  for  finding  even  near-field  solutions  to  a  thin  mask,  which  includes  evanescent  fields.  This  is 
manifested  in  solutions  to  the  following  integral: 


U (x,  y,  z) 


m 


U(x',y',0) 


aperture 


|  (1  -  ikR)  ^ 


dx'dy' 


(29) 


where  U (x,  y ,  0)  is  the  field  immediately  after  the  mask  (incident  field  multiplied  by  the  aperture  function),  k 
is  the  wavenumber,  and  R  is  the  magnitude  of  the  distance  between  any  point  (xf,  y',  0)  in  the  aperture  and  a 
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Figure  25.  Rayleigh-Sommerfeld  diffraction  intensity  patterns  for  coherent  proximity 
lithography 

Computed  Rayleigh-Sommerfeld  diffraction  intensity  patterns  of  a)  amplitude  mask ,  and  b) 
tt- shift  phase  mask  with  a  period  of  of  1.2  pm.  The  computation  window  is  20  pmx5  pm,  and 
the  Talbot  Plane  is  at  7.89  pm  below  the  mask,  defined  by  the  i-line  illumination. 


point  in  the  final  plane  (x,  y ,  z).  Equation  29  can  be  thought  of  as  the  Huy  gen’s  principal,  wherein  the  first 
Rayleigh-Sommerfeld  solution,  which  utilizes  the  Green’s  function  G(T^)  =  exp{ikR)  /  R  —  expfkR')/ R' . 
The  radiation  pattern  associated  with  this  source  is  given  by: 


Hi 


1  d  feikr\ 
2 frdz  \~R  ) 


(30) 


While  this  scalar  solution  is  fairly  straight  forward,  it  still  requires  very  large  amount  of  computation, 
typically  ^1010  inner  loop  computations  which  grow  as  total  number  points  squared. 

We  see  in  Figure  25a,  which  has  a  plane- wave  illuminated  amplitude  grating  (aperture  is  either  0  or  1),  that 
this  mask  has  about  1 -micron  of  depth  accuracy  in  tolerance  that  will  reliably  create  the  mask.  In  Figure  25b, 
the  amplitude  mask  is  replaced  by  an  equivalent  phase  mask,  where  the  aperture  either  transmits  1  or  -1.  The 
latter  example  is  more  tolerant  to  depth  errors,  accuracy  of  about  4  pm  around  3  pm  depth.  This  depth  of 
field  increase  is  mainly  due  to  the  period  doubling  a  result  of  the  requirements  of  phase  masks. 

The  features  fade  near  the  edges  of  the  computation  window.  This  fading,  evident  in  Figure  25,  is  because 
Talbot  theory  relies  on  infinite  periodicity.  However  our  calculations  don’t  assume  reflective  boundary 
conditions,  the  edge  of  the  computation  window  is  also  the  edge  of  the  aperture.  We  will  therefore  see  this  in 
the  actual  mask  patterns,  although  because  they  will  be  so  close  to  the  edge  of  the  laser  cavity  that  we  believe 
this  won’t  result  in  deleterious  effects  within  the  laser  cavity. 
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a)  b) 

Figure  26.  Microscope  image  of  best  exposure  result  of  e-beam  lithography 

Dose  of 300  pC/cm2.  At  this  rate,  the  approximate  exposure  of  complete  devices  is  just  under  1 
hour. 


As  an  alternative,  these  patterns  may  be  written  by  direct  means  like  electron-beam  lithography.  This 
technique  is  best  suited  for  small  patterns  of  “reasonable”  size  (100  pmx  100  /mi),  and  critical  dimensions  in 
the  range  of  10  nm.  We  don’t  question  that  this  technique  can  achieve  features  of  0.6  /im,  but  the  large  laser 
dimensions  may  be  complicated  by  the  following: 

stiching  Because  e-beam  lithography  is  optimized  for  writing  small  features,  the  typical  field  of  view  is 
500  /mix 500  /ma.  Because  our  lasers  are  larger  (about  4  mmx  1  mm),  this  means  we  must  stitch 
together  several  fields  of  view,  keeping  the  grating  in  each  field  in  phase  relative  to  its  nearest 
neighbors.  This  should  be  possible,  although  the  mechanical  move  between  fields  has  a  certain 
amount  of  error  associated  with  it  that  must  be  evaluated. 

write  time  The  large,  high-density  pattern  means  we  run  the  risk  of  very  long  pattern  write  times.  Moreover, 
because  the  grating  pattern  has  curved  lines,  and  the  grating  software  breaks  the  software  into 
Cartesian  boxes,  the  curved  lines  may  take  longer  to  write. 

Figure  26  shows  our  initial  test  pattern  done  to  determine  the  approximate  dose  necessary  to  create  the 
SE-DFB  patterns  we  are  after.  ZEP-520A  e-beam  resist  was  spun  onto  GaSb  substrate  samples  of  25  mm2, 
and  then  test  patterns  with  doses  of  200,  300,  400,  450,  500,  and  550  pC/ cm2  were  exposed.  After 
developing  in  HR  ZEP  developer,  the  patterns  were  inspected  with  optical  microscope.  While  doses  of 
200-450  pC/cm 2  were  all  acceptable,  300  pC/cm 2  had  the  best  fidelity,  while  500  and  550  pC/cm 2  patterns 
had  sufficient  undercutting,  a  result  from  back-scattered  electrons,  which  caused  large  portions  of  the  patterns 
to  be  washed  away.  This  calibration  was  necessary  for  the  approximate  patterns  we  plan  on  using,  but 
because  different  substrates  have  different  electron  back-scatter  patterns,  this  calibration  was  essential. 
Moreover,  we  are  pleased  to  now  know  that  each  pattern  approximate  write  time  is  about  1  hour,  which  is 
actually  significantly  less  than  feared. 

The  next  step  is  to  transfer  these  patterns  into  the  wafer  using  inductively-coupled  plasma  reactive  ion  etch 
(ICP).  Progress  on  this  effort  is  reviewed  in  Section  4.4.3. 
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4.4.3  Dry  Etching  GaSb  Materials 


To  remind  the  reader,  for  the  two  designs  found  thus  far,  we  would  like  to  etch  grating  teeth  (each  being 
approximately  600  nm  wide)  down  into  GaSb  about  400  nm.  For  a  controllable  an  repeatable  etch,  we  would 
like  this  etch  to  be  1-3  minutes,  resulting  in  an  etch  rate  of  nominally  0.15  ^m/min.  As  evidenced  by 
Figure  27,  we  have  taken  a  typical  dry  etch  (begun  from  an  existing  recipe  in  our  chamber).  While  the 
etch-rate  is  quite  fast,  the  fractional  roughness  (RMS  surface  roughness,  as  determined  by  AFM,  as  a  fraction 
of  the  etch  depth)  is  very  large.  The  resulting  etch  surface  roughness  is  within  the  optical  cavity  and  will  act 
as  scattering  centers  which  create  a  loss-term,  as  well  as  deteriorate  the  beam  resulting  laser  quality.  We 
perimetrically  changed  the  chamber  pressure,  as  it  also  imparts  characteristic  features  controlling  the  etch 
rate  and  roughness.  At  low  chamber  pressures,  there  are  insufficient  chemicals  to  etch,  so  the  etch  rate  drops. 
At  large  chamber  pressures,  the  reactive  species  mean-free  path  is  too  short,  and  etch  rates  again  suffer. 
Between  these  extremes,  there  is  a  maximum  etch  rate.  Because  below  this  maximum  point,  etching  is  often 
ballistic  (physical),  the  surfaces  can  be  quite  rough,  and  the  fractional  roughness  generally  drops  off  at  larger 
chamber  pressures  because  there  are  sufficient  surface  interaction  times  to  allow  for  chemical  “clean-up”. 
Perhaps  we  see  these  features  in  Figure  27,  but  it  is  difficult  to  discern.  The  etch  rate  peak  may  be  likely 
found  at  40  mT,  but  the  roughness  does  not  drop  as  dramatically  as  one  may  expect. 

Follow-ons  to  this  work  indicate  that  lower  fractions  of  BCI3  gas  (including  no  BCI3)  have  shown  much 
more  promising  results.  While  the  etch-rate  drops,  so  does  the  surface  roughness. 


Etch  rate 


chamber  pressure  (mT) 


Figure  27.  Etch-rates  and  etch  roughnesses  of  GaSb  etched  using  ICP 

Recipe  using  C/3  .'C/2  etch.  This  recipe  has  control  over  the  etch-rate,  although  the  roughness 
variation  is  inadequate  for  a  sufficient  grating  etch.  The  unfilled  data  indicates  approximate 
values. 
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5.0  BROADBAND  NONLINEAR  COMPONENTS 


Currently  VECSELs  employ  semiconductor-based  saturable  absorbers  which  have  a  narrow  tuning  range  and 
require  complex  fabrication  procedures.  By  developing  a  graphene-based  saturable  absorber,  one  can  take 
advantage  of  its  zero  bandgap  structure  and  therefore  its  frequency-independent  absorption  as  well  as  its 
thermal  and  mechanical  stability  to  passively  mode-lock  lasers  over  a  wide  frequency  range  -  potentially 
through  the  terahertz  regime.  Here  we  report  on  recent  developments  in  graphene  mode-locking  of  VECSELs, 
specifically  for  high  power  operation.  Further,  we  present  work  done  on  the  nonlinear  optical  properties  of 
graphene,  specifically  pertaining  to  the  development  of  saturable  absorbers  as  well  as  optical  limiters. 

From  a  practical  standpoint,  semiconductor-based  saturable  absorbers  require  complex  fabrication 
techniques,  careful  band  engineering  and  have  a  very  narrow  bandwidth.  To  overcome  these  issues  graphene 
has  emerged  as  a  promising  saturable  absorber  material  to  mode-lock  lasers.  Although  significant  work  has 
been  done  to  integrate  graphene  into  fiber  and  solid  state  lasers,  there  has  been  limited  headway  in  integrating 
graphene  in  a  VECSEL  due  the  cavity’s  sensitivity  to  losses.  [39-44] 

Single  layer  graphene  consists  of  a  layer  of  hexagonally  packed  carbon  atoms,  0.3  nm  thick  with  a  linear 
absorption  of  2.3  %.[45,  46]  This  single  layer  of  carbon  shows  a  nearly  featureless  absorption  spectrum  from 
the  visible  to  the  Terahertz  due  to  its  zero-bandgap  and  his  linear  dispersion  which  thereby  allows  intraband 
transitions  at  potentially  all  wavelengths  These  carriers  once  excited,  thermalize  extremely  fast  (in  the  order 
of  tens  of  femtoseconds)  to  the  bottom  of  the  band  thus  preventing  additional  absorption.  Hence  when 
irradiated  with  high  fluences  graphene  undergoes  an  optically  induced  transparency  which  aids  in  the 
creation  of  ultrashort  pulses  when  integrated  into  a  laser  cavity.  On  the  other  hand,  functionalized 
graphene-based  materials  display  broadband  optical  limiting  due  to  two-photon  absorption  (TPA),  reverse 
saturable  absorption  (RSA)  and  nonlinear  scattering  (NLS). [47-51]  Suppression  of  the  output  fluence  has 
been  shown  mostly  in  graphene-based  solutions  which  allow  the  possibility  of  a  strong  nonlinear  component 
due  to  scattering.  However  since  such  materials  in  suspension  form  are  not  practical  components  for  to 
integrate  into  a  structure  the  challenge  has  been  to  develop  graphene  films  in  solid  matrices  [52]  which  are 
better  suited  for  a  device  platform. 

Over  the  last  decade,  graphene  has  emerged  as  an  ideal  material  in  optoelectronic  devices  ranging  from 
saturable  absorbers  for  mode-locking  lasers  to  optical  modulators  and  photo-detectors. [46,  53-55]  Further 
graphene-based  materials  exhibiting  tunable  system  parameters  have  been  shown  to  exhibit  enhanced 
electrical  and  thermal  conductivity,  mechanical  stability  and  strong  optical  phenomena.  A  graphene-based 
material  will  hence  be  especially  attractive  as  a  saturable  absorber  in  VECSELs  due  to  its  fast  relaxation 
times,  and  strong  nonlinearity  over  a  wide  wavelength  range  potentially  from  visible  to  the  terahertz  regime. 
However  since  a  single  layer  of  graphene  has  a  significantly  high  absorption  when  taking  into  account  its 
thickness  its  integration  into  a  VECSEL  cavity  which  is  unable  to  sustain  high  loss  has  been  a  challenge. 
Broadband  optical  limiters  on  the  other  hand  would  be  able  to  suppress  high  fluences  at  various  wavelengths 
while  being  optically  transparent.  Optical  limiting  materials  are  extremely  useful  for  systems  vulnerable  to 
laser  damage  such  as  sensors,  eye  protection  etc.  In  this  effort  we  therefore  aim  to  study  various  graphene 
materials  which  exhibit  a  strong  nonlinear  optical  response  together  with  thermal  stability.  The  aim  is  to  try 
and  manipulate  these  remarkable  optical  properties  to  realize  nonlinear  devices,  particularly  saturable 
absorbers  and  optical  limiters. 
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5.1  Graphene  Polymer  Composites 

An  overview  of  the  optimized  processing  techniques  employed  to  obtain  a  functionalized  graphene  solution 
and  nanocomposite  are  described.  In  addition,  we  report  on  preliminary  optical  studies  on  various 
graphene-based  samples.  The  samples  discussed  here  consist  of  a  water-based  graphene  solution, 
graphene-polymer  composite  solution  and  thin  films.  Dependence  of  factors  such  as  concentration, 
deposition  techniques  and  surface  quality  on  linear  and  nonlinear  optical  properties  are  studied.  Furthermore 
we  demonstrate  the  ability  to  tune  and/or  control  linear  and  nonlinear  optical  properties  by  varying 
concentration,  host  matrix  and  path  length. 


5.1.1  Preparation  Techniques 

Our  first  step  was  to  create  a  graphene-polymer  nanocomposite.  This  entailed  optimizing  processing 
techniques  to  ensure  complete  solubility  of  graphene  in  water.  A  schematic  outlining  preparation  steps  is 
shown  in  Figure  28.  The  first  step  was  to  dissolve  a  dispersant  solution  SDS  in  water.  Graphene  was  then 
added  and  the  solution  was  further  sonicated  until  a  visibly  dissolved  solution  was  obtained.  Two 
concentrations  of  graphene  in  SDS-water  are  studied  (a)  1.7  mg/ml  and  (b)  6.4  mg/ml.  Solubility  was 
confirmed  when  the  graphene  did  not  settle  at  the  bottom  of  the  container  for  weeks. 

Simultaneously,  the  polymer  solution  was  prepared  by  dissolving  PVA  in  water  at  a  concentration  of 
24  mg/ml.  We  then  prepared  various  solutions  with  different  concentrations  of  graphene  in  the  PVA  solution. 
As  a  starting  point  three  solutions  were  prepared  where  the  ratio  of  the  graphene-SDS  solution  to  the  PVA 
solution  varied.  This  was  done  to  study  the  effect  of  graphene  concentration  on  the  linear  and  nonlinear 
optical  properties  of  the  composite.  In  addition  to  the  solutions,  spin  coated  and  dropcoated  samples  were 
also  prepared  using  the  graphene-polymer  composite  solutions.  These  were  spin/drop  coated  on  a  substrate 
and  heated  at  100°C  for  ten  minutes  to  evaporate  any  remaining  solvent.  Furthermore  a  graphene-oxide 
solution  was  also  prepared  by  dissolving  graphene  oxide  in  an  SDS-water  solution  at  a  concentration  of 
0.8  mg/ml.  The  graphene-oxide  due  to  its  -O-  bonds  dissolved  in  the  water  fairly  readily[56,  57]. 


(a) 


Water  and  SDS  PVA  and  water 


/j\  /|\  ^  ^  ^ 

Sonication 
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Figure  28.  Schematic  showing  preparation  of  the  graphene-PVA  composites. 

(a)  Two  separate  solutions  of  SDS  in  water  and  PVA  in  water  are  prepared,  (b)  graphene  is 
added  to  the  SDS/water  mixture  and  sonicated  (c)  PVA  is  then  added  to  the  graphene-SDS/water 
mixture  in  various  amounts  to  produce  three  different  graphene-PVA  concentrations. 
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Wavelength  (nm) 


Figure  29.  Optical  transmission  of  graphene-based  solutions 

Graphene- SDS  concentration  of  (a)  1.7  mg/ml  in  5  mm  path  length  cuvettes  and  (b)  6.4  mg/ml 
in  2  mm  path  length  cuvettes.  It  can  be  seen  that  by  increasing  the  concentration  of  graphene  in 
our  solutions  we  can  tune  the  transmission  of  our  samples. 


5.1.2  Graphene  Optical  Characterization  Methods 

Transmission  measurements  are  first  carried  out  on  the  graphene-based  solutions.  The  measurements  were 
carried  out  using  5  mm  and  2  mm  path  length  cuvettes  and  are  shown  in  Figures  29a  and  29b  respectively. 
Due  to  graphene’s  universal  absorption[46,  57,  58]  we  can  see  that  over  a  wide  range  both  the  spectra  of 
graphene  composites  remain  fairly  featureless.  Also,  as  expected  we  observe  that  as  the  amount  of  graphene 
increases,  so  does  the  absorption.Therefore  by  increasing  the  concentration  of  graphene  we  can  manipulate 
the  linear  transmission  of  our  samples. 

Solid  samples  were  also  prepared  by  spin  coating  the  graphene-PVA  composites  on  a  GaAs  substrate.  The 
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Figure  30.  Raman  spectra  of  graphene  composite  films 

Raman  spectra  of  (a)  25%  graphene-SDS  solution  in  PVA  (black)  and  PVA  with  no  graphene 
(red)  and  (b)  50%  graphene-SDS  solution  in  PVA  (black)  is  shown.  For  both  samples  we  can 
deduce  good  surface  quality  due  to  the  low  intensity  of  the  D  peak 
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Figure  31. 1-scan  measurement  technique  schematic 

Method  used  to  conduct  optical  limiting  experiments.  Unlike  conventional  z-scan  method,  here 
the  sample  remains  fixed  and  intensity  is  changed  via  a  variable  attenuator. 


quality  of  graphene  deposition  on  these  were  studied  by  conducting  Raman  measurements  at  488  nm 
excitation  and  50  X  objective.  The  Raman  spectra  is  shown  in  Figures  30a  and  30b.  As  expected  for  carbon 
based  materials  there  exist  three  prominent  peaks  [59,  60].  The  Raman  spectra  of  graphene-based  materials 
includes  a  D  peak  (absent  from  defect-free  graphene),  the  G  peak  located  at  ^  1580  cm-1  and  the  2D  peak 
observed  ~  2700  cm-1.  The  presence  of  a  D  peak  (occurring  due  to  first-order  zone  boundary  phonons) 
indicates  a  disordered  carbon  lattice  due  to  structural  defects,  amorphous  carbon  or  edge  effects.  The  G  band 
results  from  the  in-plane  vibration  of  carbon  atoms  and  is  the  most  prominent  feature  of  graphitic  material. 
The  lineshape  of  the  2D  peak  (second-order  zone  boundary  phonons)  as  well  as  its  intensity  relative  to  the  G 
peak  can  be  used  to  characterize  the  number  of  layers  ’n’  of  graphene  for  n  <  6  [57].  Single  layer  graphene 
is  usually  characterized  by  a  sharp,  symmetric  Lorentzian  peak  and  a  FWHM  <30  cm-1.  As  the  number  of 
layers  increase  the  2D  peak  becomes  broader,  less  symmetric  and  decreases  with  intensity.  Multilayer 
graphene  samples  are  characterized  by  a  multiple  Lorentzian  fit  and  a  FWHM  >50  cm-1  [61]. 

For  these  samples  we  see  that  the  D  peak  is  significantly  lower  than  the  G  peak  thereby  indicating  very  good 
surface  quality.  In  Figure  30a  we  see  the  Raman  spectra  of  a  spin-coated  PVA  film  (red)  which  unlike  the 
graphene-PVA  films  shows  no  prominent  features.  The  FWHM  of  the  2D  peaks  in  both  Figures  30a  and  30b 
are  >50  cm-1  indicating  multiple  layer  graphene.  Although  the  films  visibly  look  fairly  uniform  it  should  be 
noted  that  graphene  is  not  evenly  distributed  across  the  substrate  and  such  a  Raman  spectra  attained  from 
very  localized  regions  with  a  noticeable  cluster  of  graphene  flakes. 

With  the  development  of  well  dissolved  graphene  solution  we  then  carried  out  nonlinear  optical 
measurements  to  study  optical  limiting  mechanisms  in  the  graphene  samples  in  collaboration  with 
AFRL/RXPJ.  An  optical  limiter  ideally  should  attenuate  intense  radiation  while  exhibiting 
high-transmittance  for  low-intensity  light.  Graphene-based  dispersions  have  been  shown  to  exhibit  both 
saturable  absorption  (SA)  and  RSA  .  Each  of  these  behaviors  manifests  itself  depending  on  the  environment, 
material  properties  and  operation  regime.  Carbon  and  graphene  based  samples  have  been  reported  to  show 
strong  nonlinear  optical  limiting  effects  in  the  nanosecond  and  picosecond  regime  [48,  49,  51,  62,  63].  Such 
nonlinear  behavior  has  been  attributed  to  nonlinear  absorption  (NLA),  nonlinear  refraction  (NLR),  thermally 
induced  nonlinear  scattering,  creation  of  localized  excited  states. ..etc. [51,  64].  The  optical  limiting 
measurements  here  were  conducted  on  a  setup  whose  schematic  is  shown  in  Figure  31.  The  setup  is  capable 
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Figure  32.  Plots  showing  optical  limiting  in  the  dispersions  developed  using  1.7  mg/ml 
graphene-SDS  solution. 

The  plots  (a)  showing  optical  limiting  behavior  whereby  output  fluence  decreases  with  an 
increasing  amount  of  graphene  (b)  plot  of  transmission  versus  input  fluence  which  shows  that 
although  the  onset  of  limiting  is  comparable  in  all  three  solutions,  the  input  threshold  (point 
where  transmission  drops  to  have  its  original  value)  is  the  smallest  for  the  solution  containing 
the  highest  concentration  of  graphene. 
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of  measuring  NLA  and  NLR  simultaneously.  The  fluence  out  is  measured  via  an  imaging  system  which  also 
records  the  total  transmitted  energy.  The  fluence-in  and  incident  energy  is  measured  with  another  silicon 
meter  as  shown  in  the  figure.  The  fluence-in  is  directly  proportional  to  the  measured  incident  energy.  These 
nonlinear  measurements  were  carried  out  at  785  nm  wtih  pulse  duration  5  ns.  The  sample  is  kept  fixed  while 
an  attenuator  is  used  to  vary  the  input  fluence.  The  first  step  of  measurements  were  conducted  on  the 
graphene-based  samples  with  concentration  of  1.7  mg/ml  in  a  2  mm  path  length  cuvette.  The  optical  limiting 
curves  are  shown  in  Figures  32a  and  32b. 

In  Figure  32a,  a  plot  of  output  fluence  versus  input  fluence  for  the  graphene  solutions  are  shown.  We  can  see 
that  as  the  percentage  of  the  graphene  solution  is  increased  the  output  fluence  decreases  significantly. 

Figure  32b  shows  a  plot  of  the  transmittance  versus  the  input  fluence.  Here  we  notice  that  the  onset  of 
limiting  (point  where  transmission  starts  reducing)  is  fairly  comparable  in  all  three  solutions.  This  led  us  to 
believe  that  the  onset  of  limiting  is  independent  of  the  concentration.  However  the  input  threshold  (where 
transmission  drops  to  50%  of  its  initial  value,  or  3-dB)  shifts  to  lower  values  with  increasing  graphene 
concentration.  This  is  because  higher  concentrations  of  graphene  are  able  to  block  the  light  more  effectively 
at  higher  incident  fluences  than  those  with  lower  concentration.  For  example,  the  onset  of  limiting  in  all  three 
plots  is  seen  to  be  approximately  0.4  J/cm2.  However  the  graphene-SDS  solution  which  contains  no  PVA  and 
therefore  the  highest  amount  of  graphene  has  the  lowest  threshold  fluence  of  6  J/cm2  for  which  the 
transmittance  drops  to  its  initial  value. 

Attempts  to  tune  the  threshold  value  as  well  mimic  observed  trends,  we  conducted  the  same  experiment  with 
the  higher  graphene-SDS  (6.4  mg/ml)  solution.  Similar  data  as  before  is  shown  in  Figures  33a-c.  As  before 
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we  are  able  to  reduce  the  transmittance  considerably  by  increasing  the  graphene  concentration  shown  in 
Figure  33a.  In  Figures  33b  and  33c  we  plot  the  transmittance  versus  the  input  fluence.  Here  we  notice  that 
with  increasing  concentration  the  input  threshold  increases.  This  result  is  contradictory  to  those  attained 
above  and  therefore  further  investigations  need  to  be  carried  out  to  understand  the  dependence  of 
concentration  on  the  optical  limiting  mechanism.  It  should  be  noticed  in  Figure  33c  that  for  the  higher 
concentrated  solution  the  transmittance  slightly  increases  as  the  input  fluence  increases  up  to  around 
10  J/cm2  which  we  believe  to  be  saturable  absorption.  This  may  play  a  part  in  causing  the  deviation  in  these 
samples  from  expected  results. 

To  further  demonstrate  our  ability  to  tune  the  nonlinear  behavior  in  the  graphene  dispersions  we  conducted 
nonlinear  experiments  on  one  of  the  graphene  dispersions  (25%  graphene  solution  in  PVA)  using  cuvettes 
with  different  path  lengths.  By  decreasing  the  path  length  of  the  cuvettes  we  reduce  the  effective  length  of 
our  samples  and  therefore  increased  the  output  fluence  thereby  shifting  the  input  threshold  value.  In 
Figure  34a  the  linear  transmittance  for  all  three  samples  is  shown.  As  expected,  we  get  a  higher  transmittance 
with  a  shorter  path  length.  The  data  exhibiting  the  output  fluence  and  transmission  versus  input  fluence  is 
shown  Figure  34b  and  34c  respectively.  The  data  shown  Figure  34b  is  as  expected  -  as  the  effective  length  of 
the  sample  decreases  the  output  fluence  increases.  In  Figure  34  (c)  we  see  that  as  before  the  onset  of  limiting 
is  not  significantly  affected  by  path  length  however  the  input  threshold  changes.  For  a  greater  path  length  we 
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Figure  33.  Plots  showing  optical  limiting  in  dispersions 

Power- dependent  transmission  using  the  higher  concentration,  6.4  mg/ml  graphene-SDS  solu¬ 
tion.  The  plots  (a)  showing  optical  limiting  behavior  whereby  output  fluence  increases  with 
decreasing  amount  of  graphene  (b),  (c)  Shows  change  in  the  limiting  threshold  whereby  it 
shifts  to  higher  fluence  s  for  increasing  concentration  but  as  before  the  onset  of  limiting  remains 
fairly  constant  in  all  three  samples.  For  the  higher  concentrations  (25%  (blue)  and  50%(red)) 
graphene-PVA  solution  it  should  be  noted  that  there  is  a  slight  increase  in  transmittance  initially 
with  increasing  input  fluence  which  is  believed  to  be  saturable  absorption. 
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Figure  34.  Evidence  of  variable  optical  limiting  behavior  based  on  path  length 

The  plots  (a)  shows  change  in  optical  limiting  behavior  whereby  with  with  smaller  path-lengths 
we  can  increase  the  output  and  (b)  the  input  fluence  threshold  shifts  to  higher  values  with 
decreasing  path  length. 


observe  a  higher  input  fluence  which  conceptually  was  unexpected.  Therefore  further  studies  are  required  to 
understand  the  unusual  behavior  of  these  systems. 

Although  it  is  clear  from  all  the  plots  discussed  above  that  we  can  tune  the  nonlinear  response  from  the 
graphene-based  samples  the  reasons  for  the  nonlinear  response  are  still  not  understood.  To  better  deduce  the 
reasons  for  the  nonlinear  behavior  we  analyze  the  energy  and  fluence  transmission  in  one  of  our  samples  in 
order  to  analyze  the  contribution  from  NLR.  This  plot  is  shown  in  Figure  35.  NLR  can  cause  the  fluence 
output  to  lower  out  (self  focusing)  or  raise  (defocusing).  In  Figure  35  the  fluence  is  slightly  raised  than  the 
energy  transmission  indicating  some  self  focusing  by  the  sample.  However  since  the  difference  between 
these  is  not  significant  we  can  deduce  that  the  response  is  due  to  NLA  and  any  contribution  from  NLR  is 
small.  Other  physical  mechanisms  can  also  lead  to  limiting  behavior.  For  example  since  our  solvent  choice  if 
water  which  has  a  lower  boiling  point  (relative  to  other  solvents  used  for  such  purposes)  [51]  we  can  attribute 
part  of  the  nonlinear  behavior  to  the  creation  of  microbubbles  and  microplasmas.  These  are  formed  by  rapid 
localized  heating  (due  to  strong  linear  absorption)  giving  rise  to  solvent  bubbles  which  scatter  the  incident 
light  causing  a  reduction  in  transmission.  Other  possibilities  include  the  creation  of  excited  states  (under 
excitation)  possessing  higher  absorption  cross-section  area  than  the  ground  state  absorption  cross-section 
area  thereby  leading  to  RSA. 

The  above  optical  limiting  measurements  were  then  repeated  using  an  excitation  wavelength  of  1064  nm  with 
a  16.7  picosecond  pulse  width.  A  comparison  between  the  NL  response  of  the  785  nm  (5  ns  pulse  width)  and 
1064  nm  (16.7  ps  pulsewidth)  is  shown  in  Figures  36a  and  36b.  Here  we  see  that  optical  limiting  is  sensitive 
to  laser  pulse  width.  In  addition  it  shows  the  broadband  response  of  graphene-based  systems.  Figure  36a 
shows  the  fluence  dependence  effect  with  the  transmission  versus  irradiance  plot.  Here  we  see  that  at  the 
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Figure  35.  The  energy  transmission  and  fluence  out  is  plotted  versus  input  fluence 

The  raised  fluence  out  indicates  positive  NLR  (self focusing)  -however  as  the  difference  between 
the  energy  transmitted  and  fluence  transmission  is  not  drastically  different  we  can  attribute  the 
nonlinear  response  to  NLA  with  very  little  contribution  from  NLR. 
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Figure  36.  The  energy  transmission  for  various  pulse- widths  and  wavelengths 

(a)  normalized  transmission  ofps-scale  pulse  at  1064  nm  compared  to  ns-scale  pulses  at  785  nm 
plotted  as  a  function  of  peak  irradiance.  (b)  same  data  as  (a),  but  plotted  as  a  function  of  pulse 
energy  density  (fluence). 
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Figure  37.  Concentration  dependance  of  linear  and  nonlinear  material  transmittance  in 
films 

(a)  linear  transmittance  of  various  concentrations  of  graphene  composite .  (b)  nonlinear 
transmission  of  same  films  showing  the  concentration  dependance  of  nonlinear  effects. 


same  irradiance  the  longer  pulse  has  a  higher  nonlinearity  indicating  that  our  system  is  fluence  dependent.  It 
is  reported  that  in  the  nanosecond  regime  the  parameters  describing  the  thermal  properties  are  more  dominant 
in  contributing  to  optical  limiting.  Therefore  in  the  nanosecond  regime  the  mechanisms  contributing  to  the 
higher  nonlinear  behavior  are  thought  to  be  due  to  nonlinear  scattering  caused  by  the  formation  of 
microbubbles  and  microplasmas.  However  in  Figure  36b  we  observe  an  irradiance  dependence  indicating 
that  sample  shows  higher  nonlinearity  for  the  shorter  pulse  (higher  irradiance).  This  in  the  picosecond  regime 
is  attributed  to  nonlinear  absorption  (NLA)  mechanisms.  Therefore  the  nonlinear  processes  occurring  here 
are  both  fluence  and  irradiance  driven  and  further  investigations  are  underway  to  understand  this  behavior 
fully  [65]. 

We  now  extend  our  experiments  onto  solid  samples.  Three  dropcoated  samples  were  prepared  using  the 
concentrations  described  above.  Once  the  substrate  was  dropcoated  the  sample  was  annealed  for  10  minutes 
at  100°C  .  Their  limiting  curves  are  shown  in  Figures  37a  and  37b.  As  before,  we  see  that  the  clamping 
levels  are  concentration  dependent.  However  it  seems  that  the  clamping  takes  place  much  more  sharply  than 
the  solution  samples.  The  difference  between  the  dropcoated  and  solution  samples  are  shown  in  Figure  38.  In 
this  plot  both  samples  are  prepared  using  the  25%  graphene  in  PVA  solution.  The  solution  is  tested  using  a 
2  mm  cuvette.  The  dropcoated  sample  thickness  is  estimated  to  be  around  10  microns.  In  Figure  38a  we 
notice  that  the  dropcoated  sample  exhibits  a  higher  absorbance.  This  could  be  due  to  aggregation  of  graphene 
flakes  on  the  substrate  surface.  However  as  seen  in  Figure  38b  we  see  that  the  solid  film  shows  higher  output 
fluence  when  compared  to  the  solution  even  though  it  shows  a  higher  linear  absorbance.  The  output  fluence 
for  the  film  could  be  greater  due  to  its  smaller  path  length  -  as  the  light  interacts  with  the  2  mm  cuvette  over  a 
longer  length.  However  in  Figure  38  (c)  the  solid  film  exhibits  a  lower  input  threshold.  This  leads  us  to 
believe  that  the  other  factors  are  affecting  the  optical  limiting  response  in  the  solid.  This  could  be  due  to 
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Figure  38.  Path  length  dependance  of  linear  and  nonlinear  material  transmittance  (com¬ 
paring  film  to  solution) 

(a)  linear  transmittance  of  various  path  lengths  of  graphene  composite,  (b)  nonlinear  output 
fluence  of  various  path-length  differences ,  and  (c)  normalized  nonlinear  transmission  showing 
the  influence  of  film-  vs  solution. 


graphene  flake  aggregation,  substrate  and/or  annealing  dependent  effects. 

5.1.3  Demonstration  of  Composite  Structures 

From  the  pioneering  work  by  Novoselov  et  al.  ten  years  ago[45]  significant  volume  of  work  has  been  devoted 
to  investigations  of  electronic  properties  of  graphene.  Only  recently  has  focus  shifted  on  realizing  the  scope 
of  graphene  for  optoelectronic  devices.  Much  of  this  attention  has  been  given  to  studying  the  nonlinear 
optical  properties  of  graphene  as  applicable  to  the  development  of  saturable  absorbers  (as  described  above)  as 
well  as  optical  limiting.  Functionalized  graphene-based  suspensions  have  been  used  to  display  broadband 
optical  limiting  as  defined  in  the  seminal  paper  by  Tutt  and  Boggess  [47].  These  optical  limiting  effects  have 
been  shown  to  occur  due  to  mechanisms  such  as  TPA,  RSA  and  NLS.  Although  strong  nonlinear  behavior  (in 
particular  optical  limiting  behavior)  in  graphene  materials  has  been  demonstrated  over  a  wide  range  of 
wavelengths  the  most  promising  results  are  mostly  in  suspension  forms.  [66-68] 

This  is  due  to  the  vulnerability  of  the  solid  host  matrices  to  incur  damage  at  high  fluences.  Further  due  to  the 
high  concentration  of  graphene  in  these  suspensions  optical  limiting  effects  tend  to  dominate  any  saturable 
absorption  due  to  a  strong  nonlinears  scattering  component,  thereby  not  allowing  broadband  tunability  of  the 
nonlinear  optical  response  in  these  samples.  In  this  work  we  expand  on  these  studies  by  developing  low-cost, 
non-toxic  graphene-based  composites  both  in  solution  and  spincoated  form.  Further  we  investigate  the 
nonlinear  optical  properties  of  these  composites  for  various  pulse-widths  (nanosecond  to  femtosecond)  over  a 
wide  range  of  wavelengths. 
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Figure  39.  Demonstration  of  tunable  nonlinear  components  with  concentration 

Shows  the  change  in  transmittance  with  irradiance  for  G-PVA  solutions  in  the  (a)  ps  (1064  nm) 
and  (b)  ns  (785  nm)  regimes  respectively.  The  onset  of  limiting  for  the  G-PVA  composites  is 
lower  in  the  ps  regime  (A)  due  to  greater  TPA  contribution  and  lower  saturable  absorption.  The 
greater  saturable  absorption  observed  for  the  high  concentrations  in  the  ns-  (b)  pulsewidth  is  due 
to  greater  ESA  contribution.  The  nonlinear  coefficients  calculated  via  numerical  fitting  (solid 
lines)  reveal  high  nonlinear  absorption  and  lower  saturation  intensity  for  high  concentration 
and  longer  pulsewidth. 
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Figure  40.  Steps  for  G-DNA  composite  processing 

The  off  white  DNA  CTMA  compound  (a)  turns  to  a  gray  compound  (b)  resulting  in  a  G-DNA 
CTMA  hybrid.  This  compound  is  dissolved  in  butanol  (c)  to  form  two  concentrations  (d)  25 
mg/ml  and  (e)  18  mg/ml  and  AFM  measurements  carried  out  show  lateral  sizes  ranging  from 
1  pm  for  non-aggregated  flakes  to  6  pm  for  clustered  graphene  flakes 


5.1.4  Graphene  DNA  Composites 

Initially  a  graphene-polymer  nanocomposite  is  developed  by  preparing  a  water-based  graphene  solution. 
Together  with  a  surfactant  this  is  dissolved  in  PVA  composite  in  varying  concentrations  and  sonicated.  The 
resulting  stable  solutions  are  tested  using  an  irradiance-scan  (I-scan)  technique  in  both  the  nanosecond  (ns) 
and  picosecond  (ps)  regime.  Such  a  composite  (in  suspension  form)  indicated  a  very  weak  saturable 
absorption  effect  followed  by  strong  optical  limiting  effects  in  the  picosecond  and  nanosecond  regime  as 
shown  in  Figure  39  respectively.  Nonlinear  coefficients  calculated  indicated  higher  nonlinear  absorption 
coefficient  in  the  nanosecond  regime  as  expected,  due  to  strong  scattering  effects.  Once  spincoated  however 
the  PVA  host  was  unable  to  provide  thermal  stability  and  damaged  at  fluences  <  1  J/cm2.  Hence  such  a 
composite  was  deemed  not  suitable  for  device  integration. 

To  try  and  alleviate  these  issues  we  attempted  to  integrated  graphene  into  a  DNA  composite  which  has  been 
shown  to  be  a  promising  material  for  several  optoelectronic  devices  as  well  as  possessing  high  damage 
thresholds.  The  biopolymer  used  is  formed  from  marine-derived  DNA  which  has  been  explored  as  a  bulk 
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Figure  41.  z- scan  characterization  of  G-DNA  films  of  various  concentrations 

(a)  Peak  intensity  dependent  z-scan  measurements  of  18  mg/ml  shows  saturable  absorption 
in  the  form  of  a  peak  in  transmittance  not  seen  for  the  DNA-CTMA  ( inset)  without  graphene. 
With  the  concentration  increased  to  25  mg/ml  shown  in  (b)  higher  pumping  results  in  a  dip  in 
transmission  whereby  optical  limiting  effects  start  dominating  in  the  sample  in  the  form  of  TP  A. 


polymer  material  for  photonic  applications  for  the  past  decade.  [69]  Numerous  devices  and  applications 
containing  the  DNA  biopolymer  have  been  demonstrated;  these  range  from  purely  electronic  applications 
such  as  field  effect  transistors  (FETs)  to  optical  applications  such  as  light-emitting  diodes  (LEDs)  and 
polymer  electro-optic  modulators.  [70-72]  The  chemical  structure  of  DNA  allows  for  attaching  particles,  in 
this  case  graphene  flakes,  through  either  intercalation  or  major/minor  groove  binding,  to  the  DNA  double 
helix  for  enhanced  material  properties.  These  factors  have  therefore  led  this  attempt  to  use  DNA  as  a  host 
matrix  for  graphene.  This  novel  graphene-based  bio-polymer  results  in  combining  the  nonlinear  behavior  of 
graphene  and  the  thermal  stability  and  robust  nature  of  a  DNA  composite.  Graphene  flakes  (obtained 
commercially)  are  dispersed  into  a  DNA- water  solution,  a  surfactant  (CTMA)  is  added  and  the  graphene 
DNA-CTMA  is  precipitated,  filtered  and  dried  out  to  form  a  hybrid  material  (Figure  40  (b)).[69,  72-74]  This 
newly  obtained  hybrid  composite  now  allows  one  to  dissolve  it  in  any  solvent.  Here  we  use  butanol  to  obtain 
uniform  spincoated  films.  Material  characterization  of  the  composite  itself  is  carried  out  and  it  is  seen  that 
the  films  are  fairly  uniform  although  some  evidence  of  aggregation  is  noted  in  the  AFM  data  shown  in  shown 
in  Figure  40  (f)  and  (g)).  z-scan  measurements  on  the  films  with  varying  concentration  (images  of  these  film 
are  shown  in  shown  in  Figure  40  (d)  and  (e))  are  carried  out  at  A=1  pm  in  the  femtosecond  regime.  The 
samples  are  irradiated  with  peak  intensities  ranging  from  226  GW/cm2  to  589  GW/cm2.  The  z-scan 
measurements  are  fit  using  the  same  model  described  by  Equation  32. 

As  shown  in  Figure  41  (a)  and  (b),  as  the  peak  intensity  is  increased  the  saturable  absorption  effect  in  the 
lower  concentration  18  mg/ml  g-DNA-CTMA  film  increases  giving  us  nonlinear  coefficient  of  0.5  cm/GW 
mdlsat=59  GW/cm2.  Unlike  lower  concentrations,  an  increase  in  graphene  results  in  the  appearance  of  a  dip 
in  transmittance  at  comparable  input  intensities.  For  the  25  mg/ml  concentration  at  IQ  =  589  GW/cm22, 
optical  limiting  starts  to  overtake  saturable  absorption  effects  resulting  in  a  higher  NLA  coefficient  of  f3e / /  = 
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12  cm/GW  and  subsequently  higher  saturation  intensities.  The  inset  in  Figure  41  (a)  shows  no  change  in 
transmittance  for  the  DNA-CTMA  film  demonstrating  that  the  nonlinear  response  of  the  g-DNA-CTMA  is 
indeed  due  to  the  embedded  graphene.  Outside  the  tunability  of  the  nonlinear  response,  perhaps  the  most 
important  feature  of  the  graphene-DNA  hybrid  is  how  it  responds  to  thermal  effects  and  is  able  to  potentially 
withstand  high  optical  fluences  without  damage.  Damage  thresholds  of  these  films  are  measured  using  a 
3.63  ns  pulse  at  a  repetition  rate  of  1  Hz.  Initially  no  evidence  of  damage  is  seen  below  178  J/cm2  although 
humidity  exposure  does  affect  this  threshold.  This  value  is  much  higher  than  reported  earlier  for  work  done  in 
testing  damage  thresholds  in  graphene  materials.  [75,  76]  Currently  more  work  is  being  done  to  improve  this 
effect  in  addition  to  reducing  aggregation  of  the  graphene  flakes.  This  will  involve  centrifuging  the  solution 
before  spincoating  to  try  and  separate  out  the  larger  particles.  In  addition,  circular  dichroic  measurements  are 
also  underway  to  study  the  preferred  binding  between  the  graphene  flakes  and  the  DNA  helix. 


5.2  Graphene  Saturable  Absorber  Mirrors 

Throughout  this  program,  a  key  development  effort  has  been  focused  on  graphene  as  a  saturable  absorber  in 
laser  cavities  to  facilitate  short-pulse  operation.  Graphene’s  extraordinary  broad  band  nature  provides  an 
ideal  platform  upon  which  we  may  expand.  The  following  introduces  the  design  and  characterization  of  the 
graphene  and  the  saturable  absorber  mirrors  (SAMs),  as  well  as  a  way-forward  in  continued  characterization 
of  both  graphene-based  and  semiconductor-based  devices. 


5.2.1  Design  of  Saturable  Absorber  Mirrors 

When  the  graphene  is  located  at  a  mirror  and  you  have  bidirectional  propagation  through  the  material,  it  is 
possible  to  utilize  the  standing  wave  to  control  the  light-matter  interaction  with  the  graphene.  As  an  example, 
the  graphene  exhibits  a  relatively  flat  linear  absorption  of  ~2.3%.  However,  when  integrated  into  a  mirror, 
we  can  adjust  the  position  of  the  extremely  thin  graphene  relative  to  the  standing  mode.  The  resulting 
wavefunction  overlap  (0  directly  affects  the  linear  loss.  Figure  42  (a)  shows  how  the  field  is  controlled  in  the 
graphene  layer.  The  confinement  can  be  adjusted  between  nearly  zero  to  4,  where  1  is  equivalent  to 
unidirectional  propagation  through  a  transmissive  graphene  layer.  At  the  center  of  the  Bragg  mirror,  the 
confinement  can  be  altered  by  adjusting  the  spacer  layer.  Figure  42  (b)  shows  how  this  is  controlled  for  a 
given  wavelength  (in  this  case,  A  ~1  /un).  While  the  spacer  layer  affects  the  reflected  mode  and  therefore  the 
standing  wave  slightly,  the  spacer  layer  thickness  reaches  maximum  confinement  when  the  antinode  is 
located  in  the  graphene,  corresponding  to  a  spacer  layer  thickness  of  A/4. 

Figure  42  (c)  shows  the  confinement  factor  of  a  graphene  saturable  absorber  mirror  (G-SAM)  designed  for 
1.04  fi m  again  as  function  of  spacer  thickness,  but  for  the  full  optical  spectral  band.  However,  it  also  presents 
the  spectral  dependance  of  the  confinement.  Figure  42  (b)  represents  a  vertical  scan  at  a  wavelength  of 
1.04  fi m.  The  optical  spectral  bandwidth  of  a  G-SAM  created  using  a  spacer  layer  of  d  thickness  that  would 
function  over  is  much  wider  in  the  anti-resonant  design  (regions  of  blue)  than  in  the  resonant  design  (regions 
of  red). 

The  reflectance  of  the  G-SAM  as  a  function  of  irradiance  (/)  can  be  expressed  as: 


R(I)  —  exp 


Oi  o 


i  +  £/A 


sat 


+  tfeffI)2L 


(31) 
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confinement  factor 


Figure  42.  Calculated  field-confinement  varying  spacer  layer 

(a)  calculated  field  in  a  DBR  structure  with  a  spacer  layer.  Adjusting  the  spacer  layer  controls 
the  field  confinement  (!;)  within  that  layer,  and  therefore  the  light-matter  interaction,  (b)  At  a 
given  wavelength,  the  confinement  can  be  controlled  through  adjustment  of  the  spacer  layer 
thickness,  (c)  the  spectral  confinement  of  the  field  within  the  graphene  as  a  function  of  spacer 
layer  thickness.  Anti-resonant  structures  ( blue  regions)  have  a  much  wider  confinement  width 
than  the  resonant  structures. 


where  I  is  the  on-axis  intensity,  a0  is  the  linear  absorption  coefficient,  /,„/  is  the  saturation  intensity,  j  / /  is 
the  nonlinear  absorption  coefficient,  and  L  is  the  actual  thickness  of  the  nonlinear  material.  The  material 
saturable  absorption  can  be  expressed  as  a  subset  of  this: 

cy 

am  =  TTfij^+^!  <32) 

The  nonlinear  reflectance  of  a  G-SAM  and  dependance  on  the  nonlinear  components  is  shown  in 
Figure  43  (a).  The  linear  loss  is  controlled  through  the  £ aQ  product.  The  modulation  depth  and  unsaturated 
losses  are  controlled  by  Isat  and  /3eff  and  their  relation  to  each  other.  The  device  Isat  is  where  the  nonlinear 
losses  are  cut  in  half,  and  the  roll-over  in  reflectance  at  high  intensity  is  controlled  by  /3e//-  We  can  control 
not  only  the  linear  reflectance  through  confinement,  but  the  nonlinear  terms  (/sat/£  and  ^/3eff  /aQ)  also 
adjust  as  a  function  of  confinement.  Figure  43  (b)  shows  the  calculation  of  the  nonlinear  reflectance  for  three 
difference  confinements.  This  amount  of  variation  for  a  single  absorber  region  is  quite  dramatic  and  reflects  a 
great  deal  of  customization  and  control  possible  through  this  method  of  graphene  confinement. 


5.2.2  Characterization  of  Materials  and  Saturable  Absorber  Mirrors 

Prior  to  use  in  an  antiresonant  G-SAM  structure,  the  graphene  layer  is  characterized  using  a  CARY® 
spectrophotometer,  a  Renishaw®  Raman  system  and  transmission  mode  z-scan.  For  material  analysis,  these 
measurements  are  performed  with  graphene  on  a  quartz  substrate,  which  does  not  interfere  with  the  Raman 
signal  or  the  nonlinear  transmission  measurements.  The  linear  transmission  for  the  graphene  (with  the  quartz 
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reflectivity 


a)  intensity 


Figure  43.  Calculated  nonlinear  reflectance  for  various  GSAM  designs 

(a)  Canonical  intensity-dependent  loss  descriptions,  including  linear  and  nonlinear  losses,  b) 
Intensity -dependent  nonlinear  reflectance  for  single-layer  graphene  with  spacer  layers  resulting 
in  varied  £. 


z-position  (cm) 


Figure  44.  Linear  and  Raman  graphene  measurements  from  the  visible  through  the  near- 
infrared  region 

(a)  Linear  measurements  indicate  a  featureless  linear  transmission  of  graphene  over  a  wide 
wavelength  range  (b)  Raman  measurements  with  a  narrow,  symmetric  2D  peak  and  R  >n/Rc  >1 
indicates  monolayer  graphene  although  some  areas  exhibit  Raman  spectra  characteristic  of 
multilayer  regions  (c)  z-scan  measurements  of  the  single  film  with  a  peak  in  transmittance 
indicates  saturable  absorption  with  Isat  of  3  GW/cm 2. 
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background  removed)  is  shown  in  Figure  44  (a).  Here  we  notice  a  flat  transmission  spectra  over  a  wide 
wavelength  range  indicating  the  broadband  nature  of  graphene  with  some  variability  depending  on  location, 
presumably  due  to  regions  of  multilayer  graphene. 

The  films  are  then  studied  using  Raman  measurements  at  514  nm,  with  a  50  x  objective  and  a  1800  1/mm 
grating.  The  Raman  spectra  of  one  such  transferred  film  is  shown  at  two  different  points  in  Figure  44  (b).  As 
expected  the  Raman  spectra  of  carbon-based  materials  exhibit  three  prominent  peaks.  The  D  peak  located  at 
^1329  cm-1  is  activated  in  the  presence  of  defects.  The  G  peak  (at  1562  cm-1)  corresponds  to  the  doubly 
degenerate  phonon  mode  at  the  Brillouin  zone  center  while  the  2D  peak  at  2684  cm-1  is  a  two-phonon  band 
which  exists  in  the  second  order  Raman  spectra  of  graphene  irrespective  of  defects. [77,  78]  Here  we  observe 
that  the  D  peak  is  suppressed  indicating  a  nearly  defect  free  graphene  surface.  The  small  D  peak  that  is 
observed  is  likely  due  to  wrinkles  or  folding  in  the  graphene  transferred  from  the  copper.  The  2D  peak  is 
symmetric,  undistorted  and  can  be  fit  to  a  single  Lorentzian  with  a  full  width  half  maximum  (FWHM)  of 
33  cm-1.  This,  together  with  the  intensity  ratio  R^dIRg  >1,  indicates  monolayer  graphene.  It  should  be 
mentioned  that  although  the  response  over  a  large  part  of  the  surface  indicated  monolayer  graphene,  certain 
areas  do  exhibit  a  lower,  symmetric  2D  peak  indicative  of  multilayer  regions.  [77-79] 

A  sub-picosecond  z-scan  at  wavelengths  ranging  from  1020  nm  -  1086  nm  is  used  as  a  means  to  probe  the 
graphene  layer  saturable  absorption.  The  z-scan  setup  consists  of  a  tunable  Ti: sapphire  laser  source  with  a 
125  fs  pulse  width  and  a  1  kHz  repetition  rate.  For  the  data  presented  in  this  summary,  the  pulse  energy  was 
0.77  /iJ  at  a  wavelength  of  1086  nm.  The  beam  has  a  Gaussian  profile  with  a  beam  waist  radius  of 
approximately  25  fi m.  The  Rayleigh  length  (zr)  is  calculated  to  be  0.25  cm.  Measurements  are  taken  in 
transmission,  where  the  peak  fluence  is  given  by  IQj  (l  +  -2/4)  .  The  z-scan  measurement  is  shown  in 
Figure  44  (c)  whereby  the  peak  in  transmission  indicates  saturable  absorption,  as  the  sample  is  translated 
through  the  focal  point.  The  transmission  varies  slightly  over  the  surface  of  the  sample  due  to  regions  of 
multilayer  graphene.  To  quantify  this  contribution,  data  is  fit  to  an  intensity-dependent  absorption  model 
consisting  of  both  saturable  absorption  and  optical  limiting  terms  of  the  form  of  Equation  32.  [80-82] 

Based  on  the  relatively  low  fluences  tested,  it  is  sufficient  to  neglect  /3e//,  although  if  a  laser  is  operated  at 
larger  fluence,  effects  such  as  two-photon  absorption  must  be  considered.  The  fit  shown  in  Figure  44  (c) 
assumes  three  layers  of  graphene  and  results  in  a  saturation  intensity  of  3  GW/cm2  which  is  comparable  to 
previous  reports.  [6,  83]  At  other  points,  z-scan  measurements  indicate  single  layer  graphene,  although  this 
data  has  larger  relative  noise.  Additional  measurements  are  carried  out  at  lower  fluences  and  various 
wavelengths  for  thorough  analysis.  Although  the  Isat  does  not  change  dramatically,  the  signal-to-noise  ratio 
over  the  sample  does  vary  with  decreased  input  fluences. 

Based  on  this  material  characterization,  we  calculate  the  optimum  Si02/Si3N4  antiresonant  structure 
demonstrated  here  utilizes  a  spacer  layer  (37  nm)  designed  to  give  a  low  field  enhancement  at  the  graphene 
interface.  Here  the  field  enhancement  at  the  graphene  surface  is  approximately  0.25  which  leads  to  a  G-SAM 
saturation  intensity  of  12  GW/cm2  (I^t^  —  I  sat  /  £  =3/0.25  GW/cm2).  [79]  This  low  confinement 

is  beneficial  for  producing  high  pulse  energies  and  damage  thresholds,  but  can  be  tailored  based  on  the 
application  and  needed  response  of  the  saturable  absorber  mirror.  As  shown  in  Figure  42  (c),  this  also 
corresponds  to  a  very  broad  nonlinear  reflectance  band,  which  can  be  utilized  to  mode-lock  a  broad  spectral 
range  of  lasers,  and  has  been  used  to  do  so  in  Reference  6. 

To  better  understand  the  VECSEL  effects  described  above,  we  are  developing  an  in-house  nonlinear 
reflectance  measurement  system  to  probe  both  the  VECSEL  chip  gain  and  nonlinear  optical  characterization 
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Figure  45.  Nonlinear  reflection  measurement  system 

angle  dependent  nonlinear  reflection  setup  currently  underway  using  a  femtosecond  high  power 
laser. 


of  the  G-SAM.  Currently  our  setup  consists  of  an  angle  dependent  reflectivity  setup  using  a  goniometer  and 
high  speed  detectors  (Setup  shown  in  Figure  45).  The  reflectivity  of  the  SAM  will  be  calculated  as  a  function 
of  incident  pulse  fluence  which  is  changed  using  a  variable  attenuator  consisting  of  a  combination  of  a 
half-wave  plate  and  polarizer.  Beam  profiling  was  carried  out  to  determine  the  fluences  which  will  be  needed 
to  saturate  the  graphene.  The  pulse  fluence  on  the  graphene  SAM  should  be  at  least  10  x  the  saturation 
fluence.  Such  a  measurement  will  allow  us  to  study  and  characterize  each  saturable  absorber  mirror  (with  and 
without  the  graphene)  to  determine  parameters  such  as  saturation  fluence  and  modulation  depth.  These 
measurements  are  extremely  useful  in  predicting  the  most  optimum  design  to  achieve  high  power,  stable 
mode-locking.  Once  the  parameters  needed  to  saturate  the  graphene  are  accurately  calculated  using  the 
nonlinear  reflection  setup  similar  intracavity  fluences  will  need  to  be  induced  in  the  VECSEL  cavity  setup. 

5.2.3  Growth  and  Fabrication 

We  first  introduce  the  gain  component  of  the  cavity  that  we  use  to  experimentally  explore  VECSEL 
mode-locking.  The  chip  is  composed  of  strain  compensated  InGaAs  quantum  wells  and  GaPAs  barriers, 
followed  by  an  AlAs/GaAs  Bragg  mirror  grown  using  metal-organic  vapor  phase  epitaxy  (MOVPE)  by  Prof. 
Wolfgang  Stolz  of  University  of  Marburg.  The  structure  is  bonded  epitaxial- side  down  to  a  diamond 
heat-spreader  and  the  substrate  is  completely  removed,  with  an  anti-reflection  coating  applied  to  the  surface 
to  maximize  the  pump  in-coupling.  The  semiconductor  structure  is  optimized  for  high  power  operation  at 
nominally  1030  nm.2 

2More  information  can  be  found  in  a  similar  structure  -  Heinen,  et  al.,  Electron.  Lett.  48,  516  (2012)  and  is  provided 
by  Prof.  Jerome  Moloney  of  U.  of  Ariz.[84] 
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Figure  46.  V-cavity  schematic  and  power  characteristics 

(a)  Micrograph  of  G-S AM  epoxy-mounted  to  a  copper  die  and  mounted  in  custom  cavity  mount. 

(b)  Shows  a  schematic  of  the  V-cavity  with  the  gain  chip  at  the  fold  and  the  G-S  AM  and 
outcoupler  at  either  end.  Such  a  cavity  with  length  8.5  cm  results  in  (c)  power  curve  for  the 
VECSEL  cavity  whereby  we  were  able  to  attain  an  average  power  of  10  W  with  sub-picosecond 
pulse  widths ,  also  shown  as  a  function  of  pump  power. 


The  G-SAM  integrated  in  the  VECSEL  cavity  consists  of  the  DBR  which  possesses  9  periods  of  alternating 
layers  of  chemical  vapor  deposition  (CVD)  silicon  dioxide  (177  nm)  and  silicon  nitride  (119  nm)  with  a 
37  nm  silicon  dioxide  spacer.  The  graphene  is  transferred  from  a  commercially  purchased  graphene  on 
copper  foil.  The  graphene-Cu  is  spincoated  with  PMMA  at  4000  rpm  for  40  seconds.  The  foil  is  then 
immersed  in  a  1  M  ferric  chloride  solution  in  water  until  the  copper  is  etched.  The  free-standing 
PMMA-graphene  film  is  then  rinsed  with  water  and  methanol  and  lifted  off  using  the  DBR.  The  sample  is 
then  heated  for  20  minutes  at  100°C  to  flatten  out  the  film.  The  PMMA  is  finally  removed  using  acetone. 
[85,  86] 

5.2.4  Graphene  Mode-locked  Laser  Characterization 

We  use  a  standard  V-shaped  cavity  to  obtain  mode-locking.  The  G-SAM  is  placed  at  one  end  while  a  high 
reflectivity  (98%)  outcoupler  with  a  10  cm  focal  length  is  placed  at  the  other,  folded  about  the  VECSEL  chip. 
A  schematic  of  the  cavity  is  shown  in  Figure  46  (a).  The  VECSEL  chip  is  mounted  on  a  heat  sink  maintained 
at  13°C  for  temperature  control,  and  is  pumped  using  a  fiber  coupled  laser  diode  operating  at  808  nm.  The 
total  length  of  the  cavity  is  8.5  cm  with  a  127  pm  spot  size  on  the  G-SAM  which  leads  to  a  maximum 
intensity  of  2.5  GW/cm2  at  the  G-SAM.  This  intensity  does  not  lead  to  any  damage  of  the  graphene  although 
damage  thresholds  lower  than  this  have  been  reported.  [87]  The  high  damage  threshold  here  we  believe  is  the 
result  of  the  low  field  confinement  at  the  graphene  layer.  The  cavity  exhibits  large  mode-locked  output 
powers,  shown  in  Figure  46  (b).  Here  we  report  the  dependence  of  output  power  and  pulsewidth  as  a  function 
of  the  pump  power,  achieving  10.2  W  out  of  a  VECSEL  cavity  without  any  significant  signs  of  roll-over.  An 
increase  in  input  power  also  results  in  a  decrease  in  the  pulsewidth  of  the  system  since  the  higher  pumping 
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results  in  greater  saturation  of  the  graphene. 

The  optical  spectrum  bandwidth  also  changes  with  input  pump  power,  maintaining  nearly  bandwidth  limited 
pulses  throughout  the  entire  region  measured.  At  powers  below  about  6  W,  the  optical  spectrum  is  primarily 
single  lobed,  with  a  single-peak  autocorrelation  function  as  shown  in  Figure  47  (a)  and  (b)  for  an  output 
power  of  2  W.  Here  the  autocorrelation  spectra  and  calculated  function  consists  of  a  single  symmetrical  peak 
fitted  to  give  a  pulsewidth  of  500  fs  which  decreases  with  higher  pumping.  The  optical  spectrum  and 
autocorrelator  spectra  for  this  cavity  are  presented  at  higher  output  power  of  9.3  W  in  Figure  47  (c).  The 
emission  spectra  as  shown  in  47  (c)  displays  an  asymmetric  optical  spectra  with  a  FWHM  of  6.92  nm 
centered  at  ^1044  nm.  The  corresponding  three-peak  autocorrelation  function  pulse  is  shown  in  Figure  47 
(d).  While  the  autocorrelation  function  does  not  provide  a  unique  pulse  shape,  one  possible  temporal  pulse 
envelope  consistent  with  such  a  function  contains  two,  350  fs  FWHM  sech2  pulses  separated  by 
approximately  650  fs.  Although  this  may  not  be  a  unique  solution,  the  spectral  phase  of  such  a  “test”  pulse  is 
given  in  Figure  47  (c).  The  center  autocorrelation  peak  is  fit  with  a  sech2  function  to  show  agreement  with  a 
pulse  width  of  353  fs.  The  function  obtained  via  the  autocorrelated  discrete  Fourier  transform  of  the  optical 
spectra,  including  test  phase  from  Figure  47  (c),  is  shown  in  the  shaded  region  in  Figure  47  (d),  which  is  very 
close  to  transform  limited,  and  is  in  good  agreement  with  the  experimental  autocorrelation  signal.  Finally 
evidence  of  mode-locking  was  verified  via  radio-frequency  (RF)  spectrum  measurements  presented  in 
Figure  47  (e)  with  a  signal-to-noise  of  50  dB  at  1.76  GHz  [<30  dB  in  Figure  47  (e)]  from  higher  noise  floor 
resulting  from  larger  bandwidth  resolution),  and  strong  harmonics  up  to  our  resolved  limit  of  26.5  GHz, 
limited  both  by  our  detector  and  RF  spectrum  analyzer.  The  RF  harmonics  vary  more  than  expected,  and  the 
implications  of  these  results  remain  unclear.  However,  this  RF  spectrum  is  generally  consistent  with  two 
simultaneous  pulses  traversing  the  cavity.  If  this  holds  true,  these  sub-picosecond  pulse  envelopes  each  have 
energies  of  approximately  2.8  nJ.  This  large  average  power  compares  favorably  to  semiconductor  saturable 
absorber  mirror  (SESAM)  mode-locked  VECSELs. 

The  G-SAM  mode-locked  VECSEL  has  a  maximum  operating  internal  cavity  peak  power  of  approximately 
2.5  GW/cm2,  while  the  effective  Isat  is  closer  to  12  GW/cm2,  once  the  optical  field  overlap  is  considered. 
Mode-locking  is  occurring  below  the  saturation  fluence,  in  the  same  operation  mode  previously  seen  in 
SESAM  mode-locked,  as  well  as  graphene  mode-locked  lasers.  [87] 

In  an  effort  to  experimentally  demonstrate  the  influence  of  the  field  confinement,  we  developed  a  G-SAM 
consisting  of  the  silicon  dioxide/silicon  nitride  mirror  and  a  86  nm  SiC>2  spacer  with  a  confinement  of  near 
unity  at  the  graphene  layer.  In  addition,  to  increase  the  fluence  at  the  graphene  layer  we  increased  the  length 
of  our  cavity  and  decreased  the  spot  size  on  our  G-SAM  because  the  losses  (both  linear  and  nonlinear) 
become  larger  and  the  effective  Isat  increases  (indicated  in  Figure  43  (b)).  Although  powers  as  high  as  before 
were  not  achieved  we  have  still  measured  higher  average  power  than  previously  published  and  were  able  to 
get  a  more  stable  and  narrower  linewidth  RF  spectrum  shown  in  Figure  48  (a)  and  single  lobed  optical 
spectrum  shown  in  Figure  48  (b).  Once  we  determine  the  nonlinear  coefficients  and  optimize  the  geometric 
parameters  of  the  various  G-SAMs  developed,  we  intend  to  place  them  in  this  VECSEL  and  demonstrate 
more  stable  mode-locking  at  high  powers. 
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Figure  47.  Optical  spectrum,  autocorrelation,  and  RF  spectrum 

(a)  Shows  the  optical  spectrum  of  the  VECSEL  cavity  operating  at  a  low  power  of2W  with 
a  spectral  bandwidth  of  3  nm  corresponding  to  single  peak  (b)  autocorrelation  trace  (dots). 
At  higher  powers  this  optical  spectrum  becomes  asymmetric  as  shown  in  (c)  with  the  optical 
spectrum  of  the  VECSEL  cavity  operating  at  9.3  W,  with  a  spectral  bandwidth  of  6.92  nm 
corresponding  to  a  nearly  Fourier  transform  limited  pulse  which  corresponds  to  (d)  the  auto¬ 
correlation  trace  (dots)  with  side  lobes.  The  FWHM  of  this  side-lobed  autocorrelation  has  a 
resulting  optical  pulse  of 353 fs  (resulting  from  sech2  fit  to  central  peak).  The  shaded  region 
indicates  the  calculated  autocorrelation  trace  from  the  corresponding  optical  spectrum  using 
the  test  phase  given  in  (c).  (e)  RF  spectrum  showing  the  first  15  harmonics. 
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Figure  48.  Narrowed  RF  spectrum  indicating  stable  mode-locking 

(a)  Narrower  RF  spectra  of  the  first  harmonic  ( <  3kHz)  and  the  single  lobed  optical  spectrum 
in  (b)  indicates  potential  stable  mode-locking  in  the  system  with  the  higher  confinement  G-SAM 
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6.0  ASSOCIATED  RESEARCH 


During  this  period,  there  was  some  research  that  resulted  as  part  of  this  project,  but  the  output  was  not 
directly  related  to  the  end-product  of  the  effort.  They  are  worth  recording  within  this  report. 


6.1  Sb-Nano  wires 

While  growing  Ill-Sb  structures,  we  encountered  a  troubling  result  -  that  is  there  were  surface-defects  after 
finalizing  the  epitaxial  growth,  especially  for  thick  structures  which  are  typically  needed  for  longer 
wavelength  devices.  As  an  example,  this  occurred  on  thick  antimony-based  structures  such  as  DBRs  requires 
an  open  antimony  shutter  to  remain  open  for  the  entire  duration  of  the  growth.  Excess  antimony  condenses 
on  several  cold  surfaces  during  this  process,  including  the  cracked  source  shutter.  At  the  conclusion  of  the 
antimomide  device  growth,  the  substrate  is  cooled  under  an  antimony  overpressure  until  the  temperature  of 
the  wafer  reaches  ^400  °C,  a  process  which  prevents  surface  deterioration  resulting  from  antimony 
desorption.  At  this  point,  the  antimony  shutter  is  closed  and  the  substrate  is  cooled  down  to  ambient  chamber 
temperature.  For  long  antimonide  runs  however,  the  accumulation  of  antimony  on  the  back  of  the  antimony 
shutter  heats  up  during  the  cooling  process,  due  to  its  close  proximity  of  the  antimony  cracking.  As  the 
backside  of  the  shutter  is  vaporizing  trace  amounts  of  antimony,  the  substrate  is  already  near  ambient 
temperature.  The  antimony  hence  crystalizes  on  the  cold  wafer  surface  resulting  in  nano-crystal  formation. 
We  have  thus  generated  DBR  samples  that  have  surface  antimony  and  some  that  are  free  of  surface  antimony, 
achieved  by  using  the  cracker  valve  to  remove  the  Sb  overpressure.  The  primary  reason  for  expecting  the 
presence  of  antimony  on  the  DBR  surface  is  the  significant  change  in  the  reflectivity  of  the  samples  for 
conditions  with  and  without  antimony. 

We  used  a  series  of  analysis  to  determine  the  result  of  this  surface  structure.  First,  atomic  force  microscopy 
was  used  to  notice  there  were  surface  structures  similar  to  nanowires,  all  aligned  with  each  other  and 
orientated  in  nominally  the  same  direction.  An  example  of  this  can  be  seen  in  Figure  49a.  To  determine  the 
composition  of  these  features,  cross-section  transmission  electronic  microscopy  (X-TEM)  is  used  to 
investigate  the  surface  of  the  DBRs,  and  is  imaged  in  Figure  49b,  which  shows  three  distinct  regions  of  the 
nanowires,  each  which  may  be  analyzed  using  cross-sectional  EDS.  The  first  region  (i)  is  the  underlying 
substrate,  which  has  signal  indicating  it  is  pure  AlSb  (in  the  example  shown  in  Figure  49  the  nanowire  was 
grown  on  AlSb).  Region  (ii)  is  an  interfacial  region  where  the  aluminum  is  quenched  and  replaced  by  Sb  to 
form  GaSb,  as  indicated  by  Figure  49c.  After  this  wetting  layer,  the  nanowires  form  through  strain-induced 
growth,  which  causes  pyramidal  structures  of  pure  Sb.  While  it  can  not  be  seen  in  Figure  49b,  each  pyramid 
is  composed  of  two  domains  growing  normal  to  the  pyramid  surfaces. 


Polarization  dependent  transmission  measurements  are  then  carried  out  using  a  UV-VIS-NIR  Cary 
spectrophotometer.  Here  the  measurement  is  taken  in  a  standard  manner  with  polarizers  placed  in  the  path  of 
both  the  sample  and  the  reference  beams.  Transmission  measurements  were  taken  with  both  the  orthogonal 
incident  polarization  directions  (which  we  refer  to  as  Tv  and  T^,  polarized  in  the  v  and  h  directions, 
respectively).  As  expected,  it  was  seen  that  transmission  measurements  taken  at  normal  incidence  on  the 
DBR  without  the  antimony  deposits  (not  shown)  show  no  difference  between  the  v  and  h  polarization 
directions.  However,  polarization  selectivity  is  seen  for  the  DBR  consisting  of  the  antimony  nano-structures 
whereby  the  transmission  intensity  of  the  sample  is  higher  when  incident  with  v  polarized  light  (wherein  2^ 
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Figure  49.  Analysis  of  Sb-nanowires  on  GaSb 

a)  surface  atomic  force  microscopy,  b)  XTEM  of  surface  structure,  c)  cross-sectional  analysis 
EDS  of  layers  within  the  structure,  .  d)Transfer  matrix-based  simulations  are  used  to  fit  the 
modified  reflectance  of  antimony-based  DBR  by  using  a  calculated  antimony-air  composite 
refractive  index,  e)  Degree  of  polarization  shows  a  higher  transmission  for  vertically-polarized 
light.  In  addition,  the  amplitude  of  the  electric  field  at  the  Sb  nano- structure  surface  coincides 
with  the  enhancement  of  polarization.  A  cross-section  schematic  of  the  antimony-nano-structure 
surface  is  inset  indicating  the  direction  of  v  and  h  polarizations. 
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is  perpendicular  to  the  Sb  nanowires  on  the  surface).  A  schematic  of  the  cross-section  of  the  antimony 
surface  and  the  polarization  directions  relative  to  the  antimony  nanowires  is  shown  in  the  inset  of  Figure  49e. 
The  degree  of  polarization  (Tv/Th)  as  well  as  the  electric-field  profile  at  the  surface  of  the  nano-structures  is 
shown  in  Figure  49e.  Here  it  can  be  seen  that  the  DBR  with  the  antimony  nanowires  primarily  favors 
transmission  of  vertically-polarized  light.  The  electric  field  profile  is  calculated  using  the  antimony-air  index 
and  thicknesses  of  the  layers.  The  field  intensity  coincides  with  the  degree  of  polarization  indicating  that  the 
optical  field  is  more  confined  within  the  Sb  nano- structure  layer,  this  allows  for  greater  light-matter 
interaction  thereby  enhancing  the  effects  of  the  nano-structures.  The  DBR  below  allows  for  the  tuning  of 
these  resonances  by  changing  layer  thickness,  indices  as  well  as  the  incident  angle.  This  indicates  that  the 
antimony  nanowires  behave  much  like  a  wire-grid  polarizer  whereby  light  polarized  perpendicularly  to  the 
antimony  nanowires  has  a  higher  transmission.  This  discriminating  of  light  polarization  however  is  not 
significant  due  to  the  possible  coupling  between  the  nanowires  due  to  their  extremely  close  proximity. 


6.2  Laser  Packaging 

Indium  is  a  key  metal  for  making  electrical  as  well  as  thermal  contact  to  a  wide  variety  of  semiconductors. 
Despite  some  disadvantages  such  as  whisker  formation  and  migration  under  lifetime  test[88],  key 
characteristics  such  as  its  mechanical  compliance  and  softness,  combined  with  extremely  high  thermal-  and 
electrical-conductivity  make  it  a  preferred  contact  metal  for  applications  where  the  solder  is  sandwiched  with 
some  force  between  a  device  and  carrier.  For  instance  indium  solder  technology  is  critical  for  bonding 
detector  arrays  to  readout  integrated  circuits  (ROIC)  to  form  focal  plane  arrays  (FPAs)[89]  as  well  as  some 
high  power  devices  such  as  power  amplifiers  or  large  area  lasers  with  heat  spreaders  [90].  In  order  to  achieve 
such  a  layer  of  indium,  several  different  techniques  have  been  employed  including  -  electroplating [91], 
vacuum  deposition  by  sputtering [92]  and  evaporation[93].  The  deposition  techniques  have  offered  better 
control  over  the  thickness  of  the  indium  film  and  are  currently  employed  extensively  in  the  semiconductor 
processing  industry.  Several  significant  challenges  posed  by  deposited  indium  on  various  substrates  include 
void  formation  which  can  be  highly  detrimental  to  thermal  transport[94],  as  well  as  the  presence  of  extensive 
crystalline  columns  on  the  film  surface  which  can  complicate  bonding  to  the  indium  surface.  It  must  also  be 
pointed  out  that  the  softness  of  indium  which  makes  it  such  a  desirable  solder  prevents  the  use  of  any 
mechanical  or  chemical  polishing  techniques  from  improving  the  surface. 

To  overcome  these  indium  thin  film  voids  and  improve  surfaces  more  amenable  to  uniform  bonding,  we  have 
explored  DC  sputtering  of  indium  onto  silicon  substrates  with  using  the  substrate  temperature  as  a  key 
parameter  in  achieving  closely  packed  and  highly  planar  films.  The  planar  surface  morphology  of  a  sputtered 
metal  is  highly  dependent  upon  the  homologous  temperature  at  which  the  process  takes  place.  The 
homologous  temperature  of  the  sample  is  defined  as  the  ratio  of  the  substrate  temperature  during  sputtering 
to  the  melting  point  of  the  metal.  It  has  been  shown  by  Grovenor  et  al.[95]  and  by  Thornton[96]  in  their 
structure  zone  models  that  the  deposition  of  sputtered  metal  films  at  homologous  temperatures  of  <0.2 
results  in  very  small  equiaxial  grains  and  a  planar  surface  morphology  whereas  at  a  much  higher  homologous 
temperature  of  >0.7  the  growth  would  be  dominated  by  very  large  columnar  structures.  The  presence  of 
these  large  columnar  structures  is  considered  to  be  the  dominant  source  of  surface  roughness  when  sputtered 
indium  is  used  as  a  solder[97].  In  most  metals  the  melting  point  is  very  high  compared  to  ambient  room 
temperature,  and  therefore  the  grain  structure  corresponding  to  zone  I  or  zone  T  is  naturally  produced  with  a 
room-temperature  platen  as  the  homologous  temperature  remains  <0.2.  In  the  case  of  indium,  the  low 
melting  point  of  429.76  K  results  in  a  substrate  temperature  of  ^86  K  necessary  to  achieve  a  homologous 
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a)  b)  c) 


Figure  50.  AFM  and  Crossectional  SEM  of  indium  deposited  above  and  below  homolo¬ 
gous  temperature 

SEM  image  of  the  samples  with  indium  sputtered  on  a  silicon  samples  at  77  K  substrate  tem¬ 
perature  and  at  room  temperature.  The  cryogenic  sample  in  part  ( b )  clearly  shows  significant 
reduction  in  the  presence  of  crystalline  indium  grains  compared  to  the  room  temperature  sample 
in  part  (c).  Part  (a)  shows  the  AFM  profile  of  the  planar  region  in  the  cryogenic  sample. 


temperature  of  0.2  and  this  can  be  achieved  by  cooling  the  substrate  with  liquid  nitrogen  during  indium 
sputtering.  This  study  involves  comparing  two  scenarios.  The  first  scenario  represents  indium  deposition 
with  the  substrate  kept  at  room  temperature  (no  active  temperature  control),  where  the  homologous 
temperature  is  ^0.70.  The  second  scenario  involves  depositing  indium  with  a  homologous  temperature  of 
about  0.17  and  is  achieved  using  liquid  nitrogen  to  cool  the  substrate  to  a  nominal  temperature  of  77  K. 
Correspondingly,  the  grain  structures  fall  within  the  zone  T  and  zone  III  regions[95].  The  objective  of  the 
study  is  therefore  to  deposit  and  characterize  sputtered  indium  on  silicon  at  cryogenic  temperatures  to 
achieve  a  highly  planar  and  closely  packed  indium  layer  thus  drastically  improving  the  quality  of  indium 
solder  for  electronic  and  optoelectronic  devices. 

At  room-temperature,  films  were  typically  8.6  pm  thick,  while  low  temperature  sputtering  resulted  in  4.8  pm 
thick  indium  films.  The  main  difference  was  due  to  indium  packing  density.  The  room  temperature  samples 
had  a  very  rough  morphology  while  the  cryogenic  samples  were  near  specular.  The  reason  for  such  an 
improvement  in  the  quality  of  the  samples  can  be  seen  in  Figures  50(b)  and  50(c)  which  show  plan- view 
SEM  images  of  the  two  samples  [Figure  50(b)  at  77  K,  Figure  50(c)  at  300  K].  For  samples  kept  at  77  K 
during  the  sputtering  process  shows  the  presence  of  some  large  surface  features,  which  are  crystalline  indium 
grains,  with  large  areas  of  planar  surfaces.  The  AFM  image  in  Figure  50(a)  shows  a  profile  of  the  flatter 
surface  and  this  measurement  registers  a  root  mean  square  roughness  of  ^20  nm.  In  contrast  the  silicon 
substrate  that  was  kept  at  room  temperature  during  the  sputtering  shows  a  much  higher  density  of  indium 
grains,  with  some  features  as  long  as  50  pm  and  having  facet  formation.  For  room-temperature  depositions, 
the  film  thickness  is  ~7.3%.  Thus  the  growth  of  the  sample  at  cryogenic  temperatures  has  moved  the  growth 
mode  into  a  much  more  planar  regime  compared  to  the  room  temperature  sample.  This  effect  is  consistent 
with  other  metallic  films  structure  zone  model[95].  However  the  issue  of  the  residual  crystalline  structures  in 
the  cryogenic  sample  needs  to  be  understood.  This  can  be  explained  through  the  more  involved  model  by 
Thornton[96]  where  the  observed  surface  is  termed  a  “transitional  surface”  laying  between  a  completely 
planar  zone  and  a  columnar  growth  mode  zone.  This  issue  may  be  resolved  in  the  future  by  further  reducing 
the  homologous  temperature  by  making  use  of  a  helium  based  cryogenic  system  or  by  modifying  the  argon 
plasma  pressure  both  of  which  have  been  effective  in  other  metal  systems[97]. 
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6.3  Intracavity  Down-conversion 


6.3.1  Motivation 

The  lack  of  high-power,  tunable,  efficient  long-wave  infrared  (LWIR)  sources  is  limiting  active  sensing 
technology  and  waveform  agility,  a  critical  need  for  LADAR  type  imaging.  Most  often,  quantum  cascade 
lasers  are  used  for  direct  generation,  which  has  only  achieved  about  100  mW  with  single-digit  efficiencies 
and  little  tunability.  Nonlinear  processes  most  often  uses  rare-earth  solid  state  lasers  and  PPLN,  which  have 
inefficient  generation  at  longer  wavelengths  and  have  relatively  narrow  tuning  bands.  This  bottleneck  is 
removed  by  the  open  cavity  design  of  the  VECSELS  providing  access  to  high  intracavity  powers  required  for 
linewidth  control  and  wavelength  tunabilty  through  intracavity  non-linear  frequency  mixing  processes. 

Based  on  this  concept,  a  blue-green  laser  with  high  power  and  stability  has  been  reported  applying  a  tunable 
intracavity  type  II  sum- frequency  generation  (SFG)  in  a  folded  T-cavity  configuration.  [98]  However  the 
critical  defense  need  of  accessing  wavelengths  in  the  MWIR  and  LWIR  spectral  region  requires  the 
application  of  intracavity  difference  frequency  generation  (DFG)  processes.  Herein  we  demonstrate  the 
efficient  generation  of  MWIR  wavelengths  applying  intracavity  DFG  in  a  T- VECSEL  cavity.  We  also  provide 
a  comparison  of  the  experimentally  obtained  efficiencies  with  theoretically  calculated  values. 


6.3.2  Results  and  Discussion 

The  T- VECSEL  cavity  consists  of  two  laser  cavities  one  of  which  is  folded  around  a  polarization  beam 
splitter  (PBS)  generating  orthogonal  polarizations  of  the  two  lasers  lasing  at  980  nm  and  1 170  nm  in  the 
cavity  as  shown  in  figure  lb  and  lc.  The  regions  between  the  PBS  and  “VECSEL  chip  1”  (s-polarization), 
and  PBS  and  “VECSEL  chip  2”  (p-polarization)  have  orthogonal  polarizations  as  depicted  in  the  schematic 
in  Figure  51.  There  is  a  collinear  region  from  the  PBS  to  a  flat  mirror  via  a  high  reflecting  concave  mirror 
where  these  two  polarizations  are  multiplexed  to  satisfy  a  type  II  phase  matching  condition  through  a  non¬ 
linear  crystal  for  DFG.  Folding  the  collinear  cavities  at  the  common  cout-coupler  allows  for  the  simultaneous 
focusing  of  the  overlapped  lasing  modes  into  a  nonlinear  crystal  thereby  facilitating  efficient  difference 
frequency  generation.  The  crystal  is  made  of  silver  thiogallate  (AGS  or  AgGaS2)  chosen  due  to  its  large 
effective  optical  nonlinearity,  wide  spectral  and  angular  acceptance,  broad  transparency  range  and 
non-critical  requirements  for  temperature  stabilization  and  vibration  control.  The  crystal  offers  a 
transparency  range  between  0.53-12  microns  spanning  over  the  frequency  band  of  the  DFG.  The  dimensions 
of  the  crystal  are  5  x  5  x  5  mm3  cut  at  0=42°  and  0=90°  in  order  to  satisfy  a  type  II  phase  matching 
condition.  [99]  Birefringent  filters  at  Brewster’s  angles  are  used  in  front  of  each  chip  for  independent  tuning 
of  wavelengths  across  its  gain  band.  The  experimentally  generated  DFG  signal  is  depicted  in  Figure  52 
showing  a  range  of  idler  wavelengths  between  5.3  /i m  and  5.5  /am. 

Based  on  the  energy  conservation  law  depicted  in  Figure  51a,  a  small  tunability  in  nanometers  in  the  pump 
wavelength  (e.g.30  nm  for  the  980  nm  chip),  while  keeping  the  signal  constant  at  1 170  nm  provides  a 
tunability  of  several  microns  for  the  idler  or  the  mid-IR  wavelength  as  shown  in  Figure  5 Id.  Similarly  a  small 
change  in  the  phase  matching  angle  in  the  crystal  also  gives  rise  to  a  broad  tunability  of  the  mid-IR 
wavelengths  as  depicted  in  Figure  51e.  The  phase  matching  condition  for  type  II  negative  case  where  type  II 
corresponds  to  the  condition  where  pump  and  signal  wavelengths  have  extraordinary  (p)  and  ordinary  (s) 
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Figure  51.  T-VECSEL  schematic,  experiment  image,  and  theory 

(a)  Optical  down-conversion  process,  (b)  schematic  of  the  T-VECSEL  cavity  depicting  the 
various  optical  components,  (c)  Experimental  set-up.  (d)  Plot  showing  the  tunability  of  the 
idler  wavelength  with  pump  wave  at  constant  signals  1.17  pm  ( red  dashed)  and  1.07  pm  ( red 
solid),  (e)  Phase  matching  curve  depicting  the  tenability  of  the  signal  (green)  and  idler  (red) 
wavelengths  with  phase  matching  angles  in  the  crystal. 
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Figure  52.  Theoretical/experimental  DFG  idler  frequencies  demonstrating  tunability 
from  5.3  and  5.5  pm. 

(a)  Power  tuning  curve  for  pump  and  signal  wavelengths  centered  at  960  nm  and  1170  nm 
respectively,  (b)  experimental  DFG. 


polarizations  respectively  and  negative  refers  to  negative  uniaxial  crystal  (AGS),  is  given  by: [99] 

cos2  9V  sin2  9P  {  n°  nf  \ 2  . 21 

——f  +  7  J  =  (  d  I  xA 

0 nP )  ( nv )  \V  XJ 


(33) 


where  n° ,  and  nes  are  the  refractive  indices  for  the  ordinary  (s)  and  extraordinary  (p)  polarizations  for 

the  pump  and  signal  wavelengths  respectively  and  6P  is  the  angle  that  the  incident  pump  wavevector  makes 
with  the  optical  axis  of  the  crystal. 

The  corresponding  power  conversion  efficiency  is  given  by:  [100] 

=  F  =  32t r2d?efflPah 

Pp  ri0c\s\jnsnpni(l  +  p) 

where  Pi ,  Pp  and  Ps  are  the  intensities  of  the  idler,  pump  and  signal  wavelengths  respectively,  l  is  the  length 
of  the  crystal,  /jl  is  the  ratio  of  the  pump  and  signal  wavelengths,  h  is  a  parameter  related  to  the  beam  waist 
and  length  of  the  crystal,  c  is  the  speed  of  light  in  vacuum,  etaQ  is  the  permittivity  of  free  space,  ns ,  np ,  rii 
are  the  average  refractive  indices  of  the  signal,  pump  and  idler  respectively.  Finally,  def  f  is  the  effective 
non-linear  coefficient  of  the  crystal  which  is  given  as: 


deff  =  c/36  sin  (2 6)  cos  (2 </>)  (35) 

A  power  tuning  curve  based  on  this  formulation  is  depicted  in  Figure  52a  showing  maximum  power 
conversion  achieved  at  around  5  fi m  which  is  in  good  agreement  with  the  experiments  (depicted  in 
Figure  52b). 

In  summary,  we  have  demonstrated  the  generation  of  mid-IR  wavelengths  using  a  unique  method  of 
non-linear  DFG  process  applying  semiconductor  VECSELs  in  a  T  shaped  optical  cavity.  A  type  II  phase 
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matching  condition  has  been  satisfied  to  generate  these  wavelengths.  The  peak  power  conversion  efficiencies 
obtained  experimentally  around  5  fim  agree  well  with  the  theory.  Furthermore  there  is  a  potential  to  obtain 
broad  tunability,  higher  output  powers,  and  extension  to  LWIR  region  with  this  method  and  is  currently  under 
investigation.  Thus  the  semiconductor  VECSEL  technology  in  combination  with  a  non-linear  frequency 
mixing  process  has  significant  potential  to  obtain  high  power/brightness  lasers  with  broad  tunability 
providing  access  to  wavelengths  which  are  difficult  to  achieve  with  conventional  sources  thereby  paving  a 
promising  technology  route  for  future  sensing  and  imaging  applications. 
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7.0  CONCLUSIONS 


There  are  many  summary  comments  that  we  can  say  about  multiple  efforts  within  this  Work  Unit.  In  the 
Type-I  quantum- well  efforts,  we  have  shown  that  we  can  quantum  design  the  Auger  recombination  through 
detailed  band  structure  engineering.  We  experimentally  realize  three  different  laser  structures  with  various 
barrier  heights  and  analyzed  the  room  temperature  power  performance  and  the  and  temperature-dependent 
threshold  variation,  for  sources  operating  at  2  microns.  We  obtain  an  increase  in  the  carrier  lifetime  of  >6x. 
Our  detailed  theoretical  analysis  shows  that  most  of  this  improvement  is  due  to  the  improved  Auger  lifetime, 
with  a  lifetime  as  long  as  ^5.7  x  that  of  the  standard  structure.  Moreover,  a  similar  concept  has  been  utilized 
to  theoretically  show  orders  of  magnitude  improvements  for  lasers  operating  at  wavelengths  above  3  fim. 
Improvements  to  efficiency  of  laser  devices  can  be  very  significant,  although  depends  highly  on  operation 
regime. 

We  have  established  antimony  nano-condensates  as  a  mechanism  for  reduced  reflectivity  from  MBE  grown 
antimonide  DBRs.  We  have  analysed  these  DBRs  using  TEM  and  AFM  and  have  concluded  that  these 
crystalline  structures  consist  of  antimony.  A  simulation  of  the  reflectivity  of  the  DBR  with  the  serrated 
antimony  features  on  the  surface  of  the  DBR  agrees  with  the  observed  changes  in  the  reflectivity  spectrum. 
Finally  an  anneal  process  at  400  °C  restores  the  reflectivity  of  the  DBR  to  its  intended  design  thus  confirming 
the  analysis. 

Because  in  high  power  lasers,  thermal  management  is  of  the  utmost  importance,  we  have  spent  some  effort  to 
create  a  new  type  of  solder  deposition,  where  we  have  deposited  indium  on  silicon  substrates  using  DC 
magnetron  sputtering  comparing  samples  at  liquid  nitrogen  cryogenic  temperature  to  those  deposited  at  room 
temperature.  The  results  of  these  experiments  are  consistent  with  structured  zone  models  for  metal  sputtering 
with  the  lower  homologous  temperature  which  leads  to  a  highly  planar  indium  film  with  a  near  specular 
surface,  while  the  higher  homologous  temperature  results  in  very  rough  morphology  dominated  by  micron 
size  crystalline  grains.  The  indium  film  on  the  cryogenic  sample  still  has  some  remaining  crystalline  features 
which  puts  the  deposition  process  in  a  transitional  zone,  between  completely  planar  and  completely  rough 
morphologies,  in  Thornton  et  al.’s  structure  zone  model.  However  the  large  planar  regions  have  a  very  low 
RMS  roughness  of  ^20  nm.  By  operating  in  a  different  zone  region,  it  is  worth  noting  that  processing 
parameters  such  as  pressure  and  plasma  power  may  be  revisited  to  achieve  faster  film  depositions. 
Furthermore,  we  expect  that  this  technique  may  be  used  to  create  smaller  grain  indium  films  with  other 
deposition  techniques  such  as  evaporation. 
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Acronym 

Description 

AFM 

AFRL/RY 

AlAs 

AlSb 

atomic  force  microscopy.  11,  12,  21,  44,  57 

Sensors  Directorate.  2,  3,  25 
aluminum  arsenide.  26 
aluminum  antimonide.  10,  12 

CVD 

chemical  vapor  deposition.  63 

DBR 

DFB 

DFG 

DoD 

DQE 

distributed  Bragg  reflector.  2,  15,  18,  59,  63 
distributed  feedback.  29 
difference-frequency  generation.  1,  71,  73 

Department  of  Defense.  2 
differential  quantum  efficiency.  33,  37 

FET 

FWHM 

field  effect  transistor.  57 

full  width  half  maximum.  37,  48,  64,  65 

G-DNA 

G-PVA 

G-SAM 

GaAs 

GaSb 

graphene  -  DNA.  iv,  56,  57 

graphene  -  poly  vinyl  alcohol  compopsite.  55 

graphene  saturable  absorber  mirror.  58,  59,  61-64,  66 

gallium  arsenide.  15,  18,  20,  21,  35,  37,  47 

gallium  antimonide.  4,  10-12,  15,  18,  20-22,  37,  40,  43,  44 

HR-XRD 

high-resolution  X-ray  diffraction,  iii,  11,  12 

ICL 

ICP 

IMF 

InAs 

InP 

IRCM 

interband  cascade  laser.  25,  26 
inductively-coupled  plasma  reactive  ion  etch.  43 
interfacial  misfit.  15,  16,  18,  21 
indium  arsenide.  10 
indium  phosphide.  20 
infrared  countermeasures.  2,  3 

LADAR 

LED 

LWIR 

laser  radar.  2,  3 
light-emitting  diode.  57 
long-wave  infrared.  71,  74 

MBE 

MOVPE 

MWIR 

molecular  beam  epitaxy.  9-1 1,  37 
metal-organic  vapor  phase  epitaxy.  62 
mid- wave  infrared.  1-3,  10,  15,  71 

NLA 

NLCSTR 

NLR 

NLS 

nonlinear  absorption.  48,  49,  51,  52,  57 

Nonlinear  Control  Strategies,  Inc..  5,  20,  25 
nonlinear  refraction.  48,  49,  51,  52 
nonlinear  scattering.  45,  54 
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ODE 

PBS 

PDE 

PID 

PL 

PVA 

QCL 

RHEED 

RMS 

RSA 

SA 

SAM 

SDS 

SE-DFB 

SESAM 

SFG 

SRH 

SWIR 

TE 

TEM 

TM 

TPA 

TSP 

UNM 

UV 

VCSEL 

VECSEL 


ordinary  differential  equation.  32,  33 

polarization  beam  splitter.  7 1 
partial  differential  equation.  29 
proportional-integral-derivative.  9 
photolumiscence.  16,  21,  22 
poly  vinyl  alcohol,  iii,  46-50,  53,  56 

quantum  cascade  laser.  4,  30,  33 

reflection  high-energy  electron  diffraction.  9-11 

root-mean-squared.  21,  44 

reverse  saturable  absorption.  45,  48,  51,  54 

saturable  absorption.  48 

saturable  absorber  mirror.  58,  62 

sodium  dodecyl  sulfate,  iii,  46,  47,  49,  50 

surface-emitting  distributed  feedback.  25,  29-31,  33,  35,  37, 43 

semiconductor  saturable  absorber  mirror.  64 

sum-frequency  generation.  7 1 

Shockley-Read-Hall.  5,  6,  13,  14,  20 

short-wave  infrared.  2 

transverse  electric.  30,  33 
transmission  electron  microscopy.  18 
transverse  magnetic.  30,  33 
two-photon  absorption.  45,  54,  55 
titanium  sublimation  pump.  9 

University  of  New  Mexico,  Center  for  High  Technology  Materials.  20,  21 
ultraviolet.  3 

vertical  cavity  surface  emitting  laser.  2 

vertical  external  cavity  surface  emitting  laser.  2,  3,  15,  18,  40,  45,  61-65,  71,  73,  74 
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